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13.  ABSTRACT  (M*MHmmJ00wonSj 

Four  research  areas  were  investigated:  1.  the  propagation  and  coupling  of  hydro- 
magnetic  waves  in  the  magnetospere.  The  coupled  hydromagnetic  equations  in  the 
dipole  model  of  the  magnetosphere  were  solved  numerically.  A  reconstruction  of  the 
long  period  waves  in  an  actual  geomagnetic  storm  was  demonstrated.  2.  non-linear 
wave-particle  interactions  in  the  magnetosphere.  Anomalous  cross-field  diffusion 
of  trapped  energetic  protons  may  result  in  proton  precipitation  in  the  equatorial 
region;  3  the  filamentation  instability  of  large  amplitude  hydromagnetic  waves  in  the 
solar  wind  plasma.  This  work  is  relevant  to  observations  of  Alven  waves  and 
magnetosonic  waves  in  space  plasmas;  4.  parametric  excitation  of  whistler  waves  in 
the  high  latitude  ionoosphere  by  a  high  frequency  heater  transmitting  from  the 
ground.  It  has  been  shown  that  whisler  waves  and  Langmuir  waves  can  be  excited  by 
I  ionospheric  heaters.  The  instability  which  produces  these  waves  offers  a  potential 
|  mechanism  to  generate  large  amplitude  long  period  waves  in  the  ionosphere. 
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I.  Introduction 

This  report  only  discusses  work  that  has  been  carried  out  during  the  funding  period 
from  Feb.l,  1990  to  Jan.31, 1991  supported  under  Grant  No.  AFOSR-88-0127.  This  grant  has 
been  awarded  to  the  Polytechnic  University  since  Feb.  1,1 988.  The  work  accomplished 
during  the  early  funding  period  has  already  been  reported  in  the  previous  two  years' 
progress  reports. 

The  investigation  covered  by  the  Grant  deals  with  the  physical  phenomena 
associated  with  waves  in  the  plasmas  of  the  ionosphere,  magnetosphere,  and  solar  wind. 
The  focus  of  the  study  is  to  explore  the  role  of  the  coupled  hydromagnetic  waves  in  the 
{  geomagnetic  substorms  of  the  magnetosphere. 

During  the  past  year,  our  work  has  achieved  many  interesting  results,  which  have 
all  been  accepted  for  publication.  Basically,  the  topics  of  investigation  are  divided  into  four 
general  categories:  (a)  propagation  and  coupling  of  hydromagnetic  waves  in  the 
f  magnetosphere.  In  this  study,  coupled  hydromagnetic  equations  in  the  dipole  model  of 
the  magnetosphere  are  solved  numerically.  As  an  example  and  practical  application  of 
this  study,  a  reconstruction  of  the  storm  time  Pc  5  waves  of  November  14-15,  1979  by  our 
numerical  solution  is  demonstrated;  (b)  nonlinear  wave-particle  interaction  in  the 
,  magnetosphere.  Energetic  protons  in  the  radiation  belts  are  trapped  by  the  earth's  dipole 
magnetic  field.  In  the  presence  of  large  amplitude  kinetic  Alven  waves,  the  trajectory  of 
these  protons  may  become  chaotic  and  anomalous  cross-field  diffusion  of  protons  may 
result  as  manifested  by  the  proton  precipitation  in  the  equatorial  region;  (c)  filamentation 
instability  of  large  amplitude  hydromagnetic  waves  in  the  solar  wind  plasma.  Large- 
-•  amplitude  Alven  waves  and  magnetosonic  waves  are  frequently  observed  in  the  high 
speed  streams  of  the  solar  wind.  Their  stability  against  a  convective  filamentation 
instability  are  investigated  thoroughly.  The  relevance  of  the  study  with  some  observations 
in  the  space  plasmas  is  also  discussed;  and  (d)  parametric  excitation  of  whistler  waves  in 
£  high  latitude  ionosphere  by  the  HF  heater  transmitted  from  the  ground.  Our  study  shows 
that  whistler  waves  (VLF  waves)  together  with  Langmuir  waves  can  be  excited 
parametrically  by  the  EISCAT's  HF  heater.  The  proposed  instability  process  is  a  potential 
mechanism  to  generate  strong  artificial  PC  waves  in  the  ionosphere. 
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Through  the  support  of  this  grant,  one  graduate  student  has  completed  the 
requirements  of  the  degree  in  August  1989.  His  Ph.D  dissertation  was  included  in  the  last 
year's  progress  report.  A  new  graduate  student  Joe  Huang  had  been  recruited  to  continue 
the  research  work  during  the  past  one  and  half  funding  years. 
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Abstract.  An  instability  that  leads  to  the 
filamentation  of  large-amplitude  Alfvgn  waves  and 
gives  rise  to  purely  growing  density  and  magnetic  field 
fluctuations  is  studied.  The  dispersion  relation  of  the 
instability  is  derived,  from  which  the  threshold  conditions 
and  the  growth  rates  of  the  instability  are  analyzed 
quantitatively  for  applications  to  the  solar  wind  plasma. 
We  have  examined  their  dependence  on  the  filamentation 
spectrum,  the  plasma  fi,  and  the  pump  frequency  and 
intensity  for  both  right-hand  and  left-hand  circularly 
polarized  Alfvfcn  waves.  The  excitation  of 
filamentation  instability  for  certain  cases  of  interest  is 
discussed  and  compared  with  that  of  the  parametric  decay 
and  modulation  instability.  The  relevance  of  the  proposed 
iustability  with  some  observations  is  discussed. 

1.  Introduction 

A  great  deal  of  interest  in  the  stability  of  the  MHD 
system  in  the  presence  of  finite  amplitude  Alfvfcn 
waves  has  arisen  in  recent  years  [Wong  and  Goldstein, 
1986;  Terasawa  et  al,  1986;  Longtin  and  Sonnerup,  1986]. 
It  is  mainly  stimulated  by  the  frequent  observations  of 
large-amplitude  hydromagnetic  fluctuations  in  the  solar 
wind  at  1  AU  [Belcher  and  Davis,  1971],  in  the  high-speed 
streams  of  the  solar  wind  [Abraham-Shrauner  and 
Feldman,  1977],  in  the  upstream  Jovian  bow  shock 
[Goldstein  et  al.,  1985],  and  in  the  interplanetary  shocks 
and  the  terrestrial  foreshock  [Vinas  et  al.,  1984;  Smith  et 
ai,  1985].  These  fluctuations  are  generally  believed  to  be 
circularly  polarized  Alfvgn  waves  propagating 
almost  exactly  field  aligned.  A  class  of  parametric 
instabilities  excited  by  the  circularly  polarized 
Alfv6n  pumps  has  thus  been  investigated.  Two 
different  instability  processes  are  currently  discussed  in 
the  literature.  One  leads  to  the  parametric  decay 
instability  [Galeev  and  Oraevskii,  1963;  Sagdeev  and 
Galeev,  1969;  Cohen  and  Dewar,  1974;  Terasawa  et  al., 
1986],  and  the  other  one  gives  rise  to  the  modulation 
instability  [Lashmore- Davies,  1976;  lonsoh  and  Ohg,  1976; 
Goldstein,  1978,  Derby,  1978;  Longtin  and  Sonnerup,  1986], 
Moreover,  the  modulation  instability  has  also  been 
analyzed  in  the  region  that  describes  the  nonlinear 
evolution  of  the  Alfvbn  waves  propagating  along  a 
static  magnetic  field  [Mio  et  al.  1976a,  b;  Mjilhus, 
1976;  Spangler  and  Sheerin,  1983].  A  derivative  nonlinear 
Schrodinger  equation  possessing  soliton  solutions  is 
derived  to  govern  the  evolution  of  nonlinear  Alfven  waves. 
On  the  contrary,  Ovenden  et  al.  [1983]  have  shown  that 
the  evolution  of  nonlinear  Alfven  waves  is  governed  by  a 
set  of  three  coupled  equations  which  in  turn  are  related 
to  the  nonlinear  Schrodinger  equation  with  known  soliton 
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solutions.  The  possible  applications  of  Alfv6n 
solitons  to  solar  and  astrophysical  plasmas  have  been 
discussed  by  Ovenden  et  aL  [1983]  and  Spangler  and 
Sheerin  [1983]. 

Sakai  and  Sonnerup  [1983]  investigated  the  effects  of 
dispersion  on  the  modulation  instability  etched  by  the 
circularly  polarized  Alfv6n  waves.  Their  analysis 
is  restricted  to  the  situations  in  which  disper-lve  effects 
are  weak  and  the  wave  number  k  of  the  sound  wave  is 
much  smaller  than  the  wave  number  1^  of  the  pump 
wave  (k  «  ko).  These  restrictions  have  been  lifted  in 
the  more  recent  work  by  Wong  and  Goldstein  [1986], 
They  study  the  dispersive  effects  on  both  the  modulation 
and  the  decay  instabilities  in  a  unified  manner  ove~  a 
wide  range  of  physical  parameters.  In  addition,  the 
results  show  that  the  maximum  growth  rate  of  the 
modulation  instability  occurs  as  k  is  comparable  to  1^. 
A  new  but  weaker  instability  existing  in  a  very  narrow 
bandwidth  near  k  to  kc  has  also  been  revealed  by 
their  analysis. 

In  this  paper,  a  filamentation  instability  which  has  not 
been  considered  by  the  previous  workers  is  analyzed.  It 
is  known  that  a  large-amplituoe,  initially  uniform  wave 
propagating  in  a  plasm*,  can  break  up  into  filaments 
because  of  the  filamentation  instability  [Schmidt,  1979; 
Kuo  and  Schmidt;  1983].  This  starts  from  small 
perturbations  in  the  plasma  density,  and  it  results  in  a 
modulation  of  the  plasma  dielectric  constant  anC  wave 
distribution,  which  in  turn  increases  the  density 
perturbations.  The  purpose  of  the  present  paper  is  to 
show  that  the  filamentation  instability  can  be  excited 
together  with  the  parametric  decay  and  modulation 
instabilities  by  the  large-amplitude  circularly  polarized 
AlfvAn  waves  in  the  solar  wind.  The  threshold 
conditions  and  the  growth  rates  of  the  instability  will  be 
determined  and  compared  with  those  of  the  decay 
instability  and  the  modulation  instability.  Some 
observations  will  be  discussed  for  corroborating  the 
predicted  characteristic  of  the  proposed  instability. 

The  organization  of  the  piper  is  as  follows.  In -section 
2  we  derive  the  coupling  equations  for  Alfvbn 
sidebands  and  the  purely  growing  magnetostatic  mode.  A 
dispersion,  relation  is  -obtained  in  section  3  and  analyzed 
for.  the  threshold  fields  and  the  growth  rates  for  various 
cases.  The  numerical  results  sue  also  given  in  section  3. 
Finally  presented  in  section  4  are  a  summary  and  brief 
discussion. 

2.  Coupling  Equations 

We  investigate  the  propagation  and  filamentation  of 
ducted,  large-amplitude,  circularly  polarized1  Alfven 
wives  in  infinite,  spatially  uniform  plasmas  embedded  in  a 

constant  magnetic  field  B0  *=  zBc.  The 
wave  magnetic  fields  are  represented  Ly.  -  '  '  • 
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where  B„  is  the  unperturbed  wave  field  intensity,  assumed 
to  be  constant  and  a  real  quantity  for  simplicity;  the  plus 
and  minus  signs  refer  to  the  right-  and  left-hand  circular 
polarizations,  respectively;  k^,  and  wc  are  the  wave  number 
and  angular  wave  frequency,  satisfying  the  dispersion 
relation  %2  =  (1*-  «V^)koJvA*  wkere*n  va  and  ft  at®  the 
Alfvfin  speed  and  the  ion  cyclotron  frequency, 
respectively. 

The  zeroth-order  velocity  responses  of  electrons  and 
ions  to  the  Alfv6n  waves  can  be  written  as 

V±  =<l±VtyVpi± 

*>  -(5^)(%^)^B5yB0)e,(V''ijt)  +  c.c.  (2) 

where  the  subscripts  e  and  i  refer  to  electrons  and 
ions. 

The  process  under  consideration  is  the  scattering  of  the 
unperturbed  Alfven  waves,  into  sidebands 
(propagating  oblique  to  the  z  direction)  by  the 
simultaneously  excited  density  perturbations-  associated 
with  the  purely  growing  magnetostatic  modes.  Let 

k  =  £k  be  the  filamentation  wave  vector, 
the  density  and  magnetic  field  perturbations  of  the  purely 
growing  raagentostatic  modes  have  the  expression 

n,  *>  n^e^coskx  and  B,  =■  ^B^coskx 

where  n,  and  B,  are  real  amplitudes  and  T  is  the  growth 
rate;  quasi-neutrality  has  beer  assumed.  The  basic 
equations  that  are  linearized  for  analyzing  the  purely 
growing  modes  include  the  continuity  equations,  the 
momentum  equations  for  both  electrons  and  ions,  and  the 
Maxwell  equations; 

9v 

1°.  +  “alT  =  0 

“  i“» +  “oir 

_  _  .*) 

i“v„  +  F«  =  -*(Vvar 
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and  =  (4xn0e/c)(vtty-v,iy)  (5) 

where  n„,  T#j,  and  m(M)  are  the  unperturbed 
plasma  density  and  the  unperturbed  electron  (ion) 
temperature  and  mass,  respectively;  the 
are  the  nonlinear  Lorentz  forces  experienced  by  electrons 
and  ions  tliat  reduce  to  the  ponderomotive  forces  in-  the 
unmagnetiicd  plasma  case.  In  terms  of  Vi±  and  Svtg 
that  represent  the  velocity  responses  of  electrons  and  iors 
to  the  unperturbed  Alfven  waves  and  the 
sidebands,  respectively,,  the  nonlinear  Lorentz  forces  are 
given  by.  [Kuo  and  Lee,  1983]. 

a  +  iMV*b> 

2)3  ' 

HV  ir(M  x  2)]) 15+  fee.  (6) 


where  the  notations  mt~^n(M)  and  ^*sl  are  used.  One 
can  show  in  a  self-consistent  way  that  the  nonlinear 
Lorentz  forces  act  only  in  the  x  direction,  i.e.. 


We  now  subroutitute  the  expressions 
yv.^Je^sinkx 

/V  JV 

”  v«  “  v.bt 

-  yE^e'^ainkx 

and  F^j  *xF^elksinix 
into  (3)-{5).  The  relation 

B£/B0  “  lV/< V  +  k3cJ)]Sf./n0 
is  obtained,  and  tbs  teruitant  equation  can  be  combined 
into  a  single  coupled  mode  equation  for  the  purely 
growing  magnetostatic  mode  : 


(I1  +  IK*vAs/(l  +  kVA/n.ft)l  +  •^c*»}®*/“o  - 

(k/M)(F,  +Fj)  (7) 

This  equation  gives  the  purely  growing  density 
perturbation  as  a  function  of  the  Alfven  waves 
and  sidebands.  Without  the  nonlinear  Lorentz  forces,  (6) 
reduces  to  the  linear  dispersion  relation  of  the 
magnetosonic  eigeamode.  The  magnetostatic  ntode  is, 
however,  a  nonli  nearly  driven  mode. 

The  Alfvfn  sidebands  are  excited  through  the 
beating  current  density  driven  by  the  Alfvfn  pump 
wave  fields  on  the  density  perturbation  of  the 
magnetostatic  mode.  The  coupled  mode  equation  for  these 
sidebands  can  be  derived  from  the  following  fluid 
equations  for  electrons  and  ions  together  with  the 
Maxwell  equations: 


^  +  ^ng^  +  n,^)  -  0 
-  Is  +  V<ng^I±  +  n.'vpii) 

(S) 

5E±  +  zB0  =  -(Ta/nce)V6a+ 

(9) 

-r  -  ^  V”o) 

+  (e/Mc)(^}i  -  x  2B0 

(9*) 

Vx5t±  -  -  Hr 

(10) 

anti 

Vx6B±  -  (4xc0e/c)[(f^ii  -  S^^) 

+.  (n,/n0)C?pj±  -  ^p^.)J 

(11) 

wbere  c,1  *  (Tt  +  Tt)/n^.  - 

-  - 

:  We  have  neglected  the  nonlinear  Lorentz  forces 
(v  *Vv)  in  (9)  and  kept  the  beating  currents  (B,Vp±)  in 
(8}  and  (11)  as  the  driving  sources  of  the  Alfvin 
sidebands;  The  underlying  reason  is;  that  the  purely 
growing  mode  is  more  effective,  in  prodncihg.density 
perturbations-  than  velocity  petturfeticro.  ‘ 


Kuo  et  aL:  Filamentation  Instability  of  A>rv6n  Waves 


9623 


Following  the  functional  dependence  of  the 
magnetostatic  mode,  we  express  the  physical  quantities  of 
the  sidebands  as 

*v«±  -  [(&v««c  4  iySv^+yicoskx 

+  i z5v'ei+ ,s i nkxje^e' ^°*  ^  +  c.c. 

5n+  =  iSn±(sinkx)e1V^0*  ^  +  c.c. 

SB+  =  [(x+ia±y)coskx-iz(k/k0)sinkx]<B+e'’te1^0*  ^ 

+  c.c.  ~  (12) 

where  the  upper  sideband  (k+  =  ke  +  k,  w+  =  wa  +  i-j) 
and  lower  sideband  (k.  =  k0-k,  w  =  wa  -( i -y) *) 
propagate  together  along  the  magnetic  field  and  form  a 
standing  wave  pattern  across  the  magnetic  field,  a± 
are  unknown  parameters  defining  the  polarizations  of  the 
Alfven  sidebands  and  _will  be  determined  shortly; 
the  expressions  of  5E+  can  be  obtained  from 
the  Faraday  law. 

Substituting  (12)  into  (8)-(  11)  and  eliminating 
5n±  from  (8)  and  (9),  we  first  obtain 

=  -(VK>)®±/B0 

~  ■*+(*%/ ko)^B+/B0 

=  KV*V««)(  VnoXBp/Bo) 

-  ( WVkHXvV kW-W  l*Wty 

T  (®B+/B0)l/f 

«^i±y  =  ^^/^((^/^(n./noltBp/Bo) 

+  (V*o)n  +  (V^^/HiX'iiVkoV-l) 

•  (i«±-VA)XC/Bo))/f 

*  =  -(Vk)(A/H.)2(a-^2/n,2)(VA5/M)) 

•  (n,/n0)(Bp/B0) 

+  (VkiXlWiV^/Boil/f  (13) 

where  f  =  1  + 

The  missing  component  £ve±i  can  be  obtained  from  the  z 
component  of  (11),  it  reads 

=  &Ut  4  tt±(kvAVni)(«B±/B0)  ( 

Tliese  relations  of  (13)  and  (14)  are  then  used  in  the 
remaining  two  equations,  i.e.,  the  x  and  y  component  of 
(11),  leading  to  the  determination  of  the  parameter 

“±  -  {(l*Vty[HVk)Vo/koV-!)] 

+  (f  ±  q/M»)(k?0vAJ/^)(l  +  k2/!-,2)} 

/  {(1/f  ±  V^IHVW^V/koV-l)] 

+  (M?i J  [fko2  vA2/  rf  j-(  kg/  k)  *(  uj,*/  kp2c,2- 1)])  ( 15) 

and  -that-  of  'the  coupled  mode  equations  for  the 
Alfvin  sidebands 


{ l-(  Vk)2(lVH))2(  l^iV^JC^/koV-D/f 

=  -(l/f  ±  V^)(ko2vA2/V)(n,/no)(Bp/B0)  (16) 

Equation  (16)  shows  the  sideband  fields  to  be  produced 
by  the  pump  field  and  the  purely  growing  density 
perturbation  of  the  magnetostatic  mode.  Thus,  equations 
(7)  and  (16)  form  a  complete  coupled  set  of  equations 
describing  the  proposed  instability  process. 

3.  Dispersion  Relation 

With  the  aid  of  (2),  (13),  (14),  and  (6),  the  coupling 
term  on  the  right-hand  side  of  (7)  can  be  expressed 
explicitly  as 

(Fe  4  F,)/M  =  ±2k{q/H)(koVA2/q)2 

•  (U^/k^HtykoV) 

•  K4{m/M)(l4a±)(l±tyf\)-(l/f)(k0/k)2 

+‘(fyt^{(142tVty 
-  (Vk)2(£y«^2(iV/^2)/f} 

•  (n./n0)(Bp/B0))(Bp/B0)  (17) 

Substituting  (17)  into  (7)  and  combining  the  resultant 
with  (16),  we  finally  obtain  the  dispersion  relation 

7*  +  k'vAJ/(,+k2vA2/f4q)+kJciI+2kJ(k0vA:</’%)2 

[( l±2V^)-( Vk^qW-D/fKBp/B,,)2 

-  -2k2(k0vAI/M>)2(koV/^2)(llfV’^) 

{ l+^/ko^t^/koVlta^m/K^f  14<«±)(  1±V^) 

-(l/nOco/kJ^fVMj)^  l^iV^lKBp/Bj,)2 

/{Hk0/k)2(q/t^)2(lTa±V^) 

(V/ko2c.2-  D/^ko V/v(A}  (18) 

Equation  (18)  is  a  general  expression  of  the  dispersion 
relation  for  filamentation  instability  of  ducted 
Alfven  waves.  It  will  be  analyzed  in  the 
following  discussion  to  determine  the  threshold  conditions 
and  the  growth  rates  of  the  instability. 

In  principle,  (18)  is  a  quadratic  equation  in  n  and  can 
be  solved  analytically.  However,  since  the  coefficients  of 
(18)  are  complicated  functions  of  the  constant  parameters 
P  =  c*,/\3a  and  K  =  k„VA/q  and  the 

normalized  variable  parameter  K.  =  k/k,  (the  notations 
used  ty  Wong  and  Goldstein  [1986]  are  adopted).  A 
numerical  analysis  of  (18)  with  different  parameters  is 
desirable.  In  doing  so,  w?  vary  K  from  0  to  2  and 
examine  the  dependences  of  the  threshold  intensity  and 
the  growth  rate  on  P  and  K  for  both  right-hand  (R-H) 
and  left-hand  (L-H)  :  circularly  . ;  polarized 
Alfvin  waves.  Presented  in  ‘^Figures  1-4  'are  i-  ;the 
dependence  of  the  normalized  threshold  intensity 
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Fig.  1.  Functional  dependence  of  the  threshold  field 
intensity  on  K  and  P  for  low-  frequency  (<c  =  0.01)  R-H 
circularly  polarized  pump. 


(ij2^)  on  K;  =  (Bp/B,,)^  is  obtained  by 

setting  t  =  0  in  (18).  Figures  1  and  2  correspond  to  the 
R-H  circularly  polarized  waves.  It  is  shown  that  the 
filamentation  instability  has  a  lower  threshold  level  for  a 
wave  with  a  smaller  frequency  (i.e.,  smaller  *).  In 
addition,  the  threshold  of  the  instability  also  varies  with 
p.  For  very  low  frequency  waves  (<c  <  0.05)  the 
threshold  intensity  increases  with  P.  However,  it 
becomes  a  decreasing  function  of  P  for  <c  >  0.1. 
This  trend  has  been  demonstrated  in  these  two  figures, 
wherein  «  «  0.01  and  0.3  are  used  as  the  representative 
parametric  values.  We  have  excluded  the  threshold  values 
in  the  neighborhood  of  X  (marked  on  the  curves)  from  the 
figures.  This  is  because  in  that  wave  number  region  the 
nonlinear  Lorentz  forces  on  electrons  and  cm  ions  are  in 
opposite  directions  with  comparable  magnitude,  so  that  the 

two  terms  proportional  to  (Bp/B0)a  in  (18)  tend  to 
cancel  to  each  otter  and  thus  the  threshold  becomes  very 
high. 

Figures  3  and  4  correspond  to  L-H  circularly  polarized 
waves.  It  is  found  that  the  instability  can  be  excited 
only  by  waves  with  higher  frequencies  (e.g.,  *  > 


Fig.  3.  Dependence  of  the  threshold  field  intensity  on  K 
and  k  for  L-H  circularly  polarized  pump  in  high  P  (» 
1.5)  plasma. 


0.15)  with  reasonable  thresholds.  The  dependence  of  the 
threshold  on  k  for  P  »  1.5  and  that  on  P  for  *  ■  0.3 
are  also  shown  in  Figures  3  and  4,  respectively.  It 
appears  that  the  instability  prefers  to  be  excited  in  the 
region  K  <  0.5.  Otherwise,  the  threshold  increases  very 
fast  with  K.  In  the  region  K  <  0.5,  the  threshold 
decreases  with  both  k  and  P  until  the  instability 
becomes  forbidden  for  the  region  of  very  small  jc  (e.g., 
the  region  below  K  «  X  marked  in  both  Figures  3  and  4 
for  k  m  0.3  and  p  *  1.5). 

We  next  present  the  growth  rates  of  the  instability. 
Shown  in  Figures  5  and  6  are  the  functional  dependence 
of  the  growth  rate  7  on  the  wave  field  intensity  ij* 
of  the  R-H  circularly  polarized  wave  for  low  P  (»  0.5) 
and  high  P  (*  1-5)  cases,  where  k  and  K  are  considered 
to  be  constant  parameters.  In  order  to  obtain  adequate 
information  for  the  dependencies  of  7  in  these  two  cases, 
two  representative  values  for  esch  of  the  two  parameters 
s  and  K  are  chosen  in  the  figure  presentation:  k  *  0.01 
and  0.3  stand  for  the  cases  of  the  low-frequency  and 
high-frequency  waves,  respectively.  We  then  use  K  »  0.1 
and  1  to  characterize  the  regions  of  large-scale  and 
small-scale  filamentation  instability,  respectively.  The 


polarized  pump. 
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*t  -  /r/k0VA  on  the  R-H  circularly  polarized 
pump  intensity  if,  where  K  >  0.1  and  (*,  p)  as 
indicated  are  set  in  the  evaluation. 


results  show  that,  in  general,  the  growth  rate  increases 
with  P  and  K  and  decreases  with  k.  For  the  L-H 
circularly  polarized  wave  case,  the  instability  can  be 
excited  only  by  the  high-frequency  waves  in  the  high  p 
plasma.  We  therefore  choose  P  »  1.5  and  a  «  0.15  and 
0.3  to  evaluate  the  dependence  of  7  on  if.  Again  K 
-  0.1  and  1  are  considered,  and  the  results  are  presented 
in  Figure  7.  In  this  case,  the  growth  rate  increases  with 
a  for  K  *  0.1  and  becomes  a  decreasing  function  of  a 
for  K  -  1. 

4.  Summary  and  Discussion 

We  have  investigated  the  filamentation  instabilities  of 
large-amplitude,  circularly  polarized  Alfv6n  waves 
propagating  along  the  background  magnetic  field.  The 
instabilities  are  excited  via  the  scattering  of  the 
unperturbed  Alfv6n  pumps  (considered  as  a  pump 
wave)  into  sidebands  by  the  density  perturbations  that  are 
associated  with  the  simultaneously  excited  purely  growing 
magnetostatic  modes.  A  four-wave  coupling  process  is 


circularly  polarized  pump  in  a  plasma  with  P  *  1.5. 


then  considered  for  the  analysis  of  the  instabilities.  The 
theory  developed  is  based  on  the  two-fluid  plasma  model 
In  general,  the  fluid  model  is  valid  when  the  condition 
l?c*t/ff  =  Ks*?/y2  «  1  is  satisfied.  In  the  present 
work,  the  maximum  value  of  the  parameter  KJ*?£/2 
for  K  «*  2,  ic  •  0.3,  and  p  •  1.5  is  0.27,  which 
reasonably  Justifies  the  use  of  the  fluid  modeL  For  the 
more  general  consideration  that  includes  the  region 
K*«0/2  >  1,  however,  a  kinetic  plasma  model 
should  be  used. 

The  nonlinear  source  for  the  Alfvtn  sidebands  is 
the  beating  current  driven  by  the  pump  wave  field  on  the 
density  perturbation  of  the  purely  growing  mode,  whereas 
the  nonlinear  Lorentz  force  (which  reduces  _  to  the 
ponderomotive  force  in  the  unmagnetized  case)  introduced 
by  the  spatial  gradient  of  the  resultant  high-frequency 
wave  field  is  the  driving  source  for  the  noooscillatory 
(i.e.,  the  purely  growing)  mode.  These  nonlinear  effects 
result  in  the  coupling  of  the  noooscillatory  mode  with  the 
Alfvin  sidebands  through  the  Alfven  pump  wave. 
A  dispersion  relation  of  the  instabilities  is  thus  derived  by 
combining  the  coupled  mode  equations.  Solving  the 
dispersion  relation,  we  have  determined  the  threshold 
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fields  and  the  growth  rates  for  the  cases  of  the  R-H  and 
L-H  circularly  polarized  pump  waves.  We  have  also 
determined  the  functional  dependencies  of  the  t '•.res  ho  Id 
field  i?th  on  the  wave  number  K  of  the 
nonoscillatory  mode,  the  P  of  the  plasma,  and  the 
frequency  k  of  the  pump  (see  Figures  1-4).  The 
dependence  of  the  growth  rate  7  on  the  pump  inteusity 
if  is  examined  for  several  sets  of  representative 
parameters  (X-AK)  (see  Figures  5-7). 

For  the  geophysical  applications  of  the  concerned 
instabilities,  we  follow  the  examples  discussed  by  Wong 
and  Goldstein  (1986],  The  finite  amplitude  Alfvfin 
waves  that  were  observed  in  the  high-speed  streams  of 
the  solar  wind  and  discussed  by  Abraham-Shrauner  and 
Feldman  [1977]  are  the  left-hand  circularly  polarized 
wave  with  *  *  0.3  and  rf  =  0.05.  As  shown  by  Wong 
and  Goldstein  [1986],  these  waves  are  stable  to  both  the 
modulational  and  the  decay  instabilities  for  (5  >  1.  From 
our  analysis,  however,  it  is  found  that  these  waves  are 
unstable  to  the  filamentation  instability  in  the  small  K 
region.  The  growth  rates  of  the  instability  evaluated  for 
P  =  1.5  are  presented  in  Figure  8;  We  therefore  predict 
the  appearance  of  nonoscillatory,  cross -field  plasma 
density  perturbations  and  magnetic  field  perturbations  in 
the  background  solar  wind  plasma.  However,  these 
predicted  cross-field  density  striations  and  magnetostatic 
structure  probably  have  not  been  observed  by  Abraham- 
Shrauner  and  Feldman  [19771.  The  difficulty  in  the  in 
situ  measurements  of  the  purely  growing  modes  by 
satellites  lies  in  the  large  velocity  difference  between  the 
solar  wind  und  the  satellite.  On  the  other  hand,  it  should 
not  be  difficult  for  the  satellites  to  observe  the  magnetic 
and  density  fluctuations  of  obliquely  propagating  sideband 
AlfVen  waves.  We  have  not  drawn  the  supporting 
evidence  from  the  literature  yet. 

Low-frequency,  left-hand  polarized  finite  amplitude 
Alfvdn  waves  have  also  been  observed  in  the 
upstream  part  of  the  Jovian  bow  shock  [Goldstein  et  aL, 
1985].  It  is  shown  from  our  analysis  that  this  wave  is 
filamentally  stable.  In  fact,  only  the  decay  instability  is 
very  weakly  unstable  for  the  case  k  *  0.02,  if  -  0.05, 
and  p  =  1.5.  Whereas  the  observation  of  right-hand 
polarized  waves  was  not  reported,  it  can  be  speculated 
from  the  linear  dispersion  relations  of  itie  Alfvdn 
waves,  because  in  such  a  low-frequency  region  the 
dispersion  relations  of  right-  and  left-hand-polarized 


Fig.  9r  Dependence  of  7  of  K  for  a  R-H  circularly 
polarized  pump  with  rf  » '0.05  and  k  °  0.02  in  a 
plasma  with  &  >•  1.5. 


Alfvdn  waves  are  almost  identical.  If  this  is 
indeed  so,  the  filamentation  instability  is  found  to  be 
unstable  against  the  right-hand  pumps  in  the  entire  K 
region  of  interest.  The  functional  dependence  of  the 
growth  rate  of  the  instability  on  K  is  presented  in  Figure 
9.  Again,  this  instability  introduces  nonoscillatory,  cross¬ 
field  density  and  magnetic  fluctuations  to  the  background 
plasma  and  thus  can  be  attributed  to  the  observed 
correlations  in  the  magnetic  and  density  fluctuations. 

For  the  other  cases  of  right-hand-polarized  waves  in 
interplanetary  shocks  and  in  the  terrestrial  foreshock 
[Vinas  et  al.,  1984;  Smith  et  al.,  1985],  the  filamentation 
instability  requires  a  threshold  power  higher  than  that 
used  by  Wong  and  Goldstein  [1986]  in  each  case. 
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Combined  Operation  of  Two  Ground  Transmitters 
for  Enhanced  Ionospheric  Heating 
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( Received  Apri*  8.  1988:  Au'epied  June  .'.  1988) 

The  combined  operation  of  an  HF  or  MF  ground  transmitter  and  a  VLF 
transmitter  tor  enhanced  ionospheric  heating  is  discussed.  The  HF  or  MF  transmitter, 
operated  in  a  pulsed  mode,  can  preferentially  produce  short-scale  density  striations 
that  can  render  the  nonlinear  mode  conversion  of  the  subsequently  launched  VLF 
waves  into  lower  h\ brid  waves.  In  addition  to  the  mode  conversion  process,  the 
\  LF  vs  axes,  if  intense  enough,  can  also  excite  meter-scale  density  striations  and 
lower  hybrid  waves  via  parametric  instabilities.  Intensified  density  striations  and 
enhanced  airglow  are  expected,  and  they  can  be  detected  by  incoherent  bachscatter 
radars  and  photometers,  respectively  The  leasibility  and  planning  ol  the  proposed 
expeiimenis  are  addressed. 

! 

I.  Introduction  : 

An  experimental  scheme  is  discussed  for  enhanced  ionospheric  healing  by  the  i 

combined  operation  of  two  powerful  ground  transmitters.  One  is  an  HF  or  MF  ; 

heater  wave,  and  the  other  is  a  VLF  heater  wave.  The  proposed  scenario  is  as  i 

follows. 

The  ionosphere  is  illuminated  by  an  HF  or  MF  wave  first  and.  then,  a  VLF 
wave  subsequently  (see  Fig.  I).  It  has  been  known  both  theoretically  and  experi¬ 
mentally  that  short-scale  (i.e..  less  than  a  few  meters)  ionospheric  irregularities  can 
be  excited  within  seconds  by  HF  heater  waves,  while  it  takes  tens  of  seconds  or 
longer  for  large-scale  (say.  hundreds  of  meters  and  longer)  ionospheric  irregularities 
to  be  generated  (see.  e.g...  Gl  REViCH.  1978:  FE.IER.  1979;  FREY  et  at...  1984:  COSTER 
ci  al..  1985).  As  shown  in  Section  2  of  this  paper,  short-scale  ionospheric  density 
striations  can  effectively  cause  the  nonlinear  scattering  of  VLF  waves  into  lower 
hybrid  waves.  Hence,  we  only  need  to  operate  the  HF  or  MF  heater  in  pulsed  mode 
to  assure  the  excitation  of  meter-scale  ionospheric  irregularities.  During  the  vertical 
ionospheric  healing,  the  HF  heater  wave  frequencies  should  be  less  than  the/oF;  for 

overdense  heating  of  the  ionospheric  F  region.  When  a  MF  heater  is  used,  the  wave  | 

frequency  is  required  to  match  the  local  electron  gyrofrequency  in  the  ionosphere  to 
produce  short-scale  ionospheric  irregularities  via  electron  cyclotron  heating  (Lee  el 
aL  1986). 
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Fig  !  The  propped  experiment  lor  enhanced  ionospheric  heating  Illumination  ol  the  ionosphere  b\ 
a  HF  01  Ml-  ground  tiansmmer  lirst.  and  then  a  \  LF  ground  transmitter  subsequently 


In  addition  to  the  nonlinear  scattering  process,  the  VLF  waves,  if  intense 
enough,  can  also  excite  lower  hybrid  waves  and  ionospheric  irregularities  via 
parametric  instabilities  (LFE  and  Kto.  1984).  The  suggested  diagnoses  for  the 
expected  ionospheric  effects  shall  be  discussed  in  Section  3.  In  principal,  the  two 
different  processes  producing  the  lower  hybrid  waves  can  be  experimentally 
distinguished. 
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2  Heater  Ware- Induced  Ionospheric  Effects 

2  /  Mode  conversion  via  nonlinear  si  altering  process 

Various  plasma  instabilities  can  be  excited  by  HF  heater  waves  during 
overdense  ionospheric  heating  In  seconds,  short-scale  ionospheric  density  striations 
can  be  generated  concomitantly  with  Langmuir  \va\es.  upper  hybrid  waves  etc 
However.  it  requires  tens  of  seconds  tor  the  generation  of  large-scale  ionospheric 
irregularities  by  sell-focusing  instability  or  thermal  filamentation  instability .  Large- 
scale  irregularities  can  gi\e  rise  to  phase  and  amplitude  scintillation  of  beacon 
satellite  signals.  These  wave  interference  patterns  may  not  be  associated  with 
significant  attenuation  of  the  radio  signals.  By  contrast,  the  short-scale  irregularities 
are  able  to  cause  anomalous  absorption  of  the  radio  signals  via  the  nonlinear 
scattering  process  described  below  that  converts  electromagnetic  wave  energy  into 
electrostatic  wave  energy  (eventually.  plasma  kinetic  energy )  m  the  ionosphere 

Lor  simplicity .  a  ducted  w  lustier  ( VLF)  wave  propagating  along  the  f  direction 
is  assumed.  viz  . 

£  ~  ( v  -  n)£  e\p[/(A  c  -  /)].  (I) 

Further,  short-scale  ionospheric  irregularities  induced  by  an  HF  or  MF  heater  wave 
are  represented  by 

on  —  <>/ie\p(/A  1  t.  (2) 

Solving  the  wave  equation,  we  can  show  that  the  ionospheric  irregularities  scatter 
nonlinearly  the  circularly  polarised  (whistler)  wave  into  an  eliiplicaily  polarized 
wave  In  other  words,  a  linearly  polarized  component  of  the  scattered  wave  is 
introduced  by  the  E\1  wave-induced  field-aligned  ionospheric  irregularities 
The  scattered  whistler  (VLF)  wave  field  can  be  expressed  as 

£  =  |(.v  *  i\  )£-  -  t'£-|exp[MAi,r  -  <••>.,/)].  (3) 

where  £,=  £.,(<>/;  ,\„)  representing  the  wave  field  intensity  of  the  circularly  polarized 
component  while  £-  represents  that  of  the  linearly  polarized  component.  A'„  is  the 
unperturbed  plasma  density.  It  is  f  und  that  £.r  and  £>  are  related  by 

£r  =  l[j-\  Bp.  (4) 

where  /=2r.  k  and  /.  -2.7  k,  are  the  scale  size  of  ionospheric  irregularities  and  the 
wavelength  of  the  whistler  wave,  respectively.  If  a  VLF  wave  at  the  frequency  of  10 
kHz  is  injected,  and  the  ionospheric  conditions  are  assumed  to  be  f  (electron 
cyclotron  frequency  j- L. 3  MHz  and  /-(electron  plasma  frequency)=6  MHz.  the 
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whistler  wavelength  in  the  ionosphere  is  lound  to  be  about  500  meters.  According  to 
Eq  (4).  dominates  over  Fr.  provided  that  Hence,  meter-scale  irregularities 
pioduced  by  HF  ol  MF  heater  wave  can  effective!)  render  the  nonlinear  mode 
conversion  of  a  \  LF  (whistler)  wave  into  a  lower  hybrid  wave.  The  parallel  and 
perpendicular  wave  vectors  of  the  lower  hybrid  wave  are  identical  to  the  wave  vector 
of  the  ducted  whistler  wave  and  that  of  the  field-aligned  ionospheric  irregularities, 
respectively 

2.2  I  LF  n  ave-exciteJ parametric  instabilities 

While,  in  principle,  the  nonlinear  scattering  process  occurs  at  any  field  intensity 
of  the  incident  VI. I*'  wave,  parametric  instabilities  can  be  triggered  by  intense 
enough  YLF  waves  as  suggested  in  Ltf-  and  Kio  (1984).  Short-scale  ionospheric 
irregularities  associated  with  zero-frequency  modes  can  be  excited  simultaneously 
with  lower  hybrid  waves  by  the  \  LF  waves.  Therefore,  the  HF  or  MF  heater 
wave-induced  short-scale  ionospheric  irregularities  can  be  further  strengthened  by 
the  \  LF  heater  waves. 

3  Diagnoses  of  Expected  Ionospheric  Effects 

It  is  generally  believed  that  spread  F  echos  in  ionograms  and  scintillation 
phenomena  are  caused  by  large-scale  (say.  hundreds  of  meters)  ionospheric  irregu¬ 
larities  The  generation  of  large-scale  ionospheric  density  irregularities,  therefore, 
can  be  sensed  by  mnosondes  and  scintillation  measurements  of  radio  stars  or  beacon 
satellite  signals.  By  contrast,  the  presence  of  field-aligned  short-scale  (a  few  meters 
and  less)  irregularities  can  be  detected  by  backscatter  radars. 

Backscatter  radar  measurements  can  also  be  made  for  detecting  the  excited 
lower  hybrid  waves  when  the  radar’s  beam  angle  is  90;  with  respect  to  the  earth's 
magnetic  field.  The  recorded  incoherent  backscatter  radar  spectrum  should  look  like 
a  double-humped  ion  spectrum  peaking  at  the  lower  hybrid  wave  frequency.  Since 
lovvei  hybrid  waves  can  accelerate  electrons  effectively  along  the  geomagnetic  field, 
airglow  enhancement  due  to  the  impact  excitation  of  neutrals  by  energetic  electrons 
can  be  expected.  Measurements  of  airglow  at  6300  A  (red).  5577  A  (green),  and  even 
shorter  wavelengths  may  possibly  be  observed.  In  addition,  the  technique  described 
in  C  ARl.sox  el  al.  ( 1982)  for  finding  the  height  distribution  of  enhanced  plasma  lines 
is  capable  of  distinguishing  further  the  electron  acceleration  by  Langmuir  waves 
and  or  upper  hybrid  waves,  and  lower  hybrid  waves  that  are  produced,  respectively, 
by  HF  and  VLF  healer  waves. 

There  are  several  ionospheric  heating  facilities  operated  at  the  frequencies  of  a 
few  MHz  in  the  HF  or  MF  band  in  the  U.S.A..  Europe,  and  the  U..S.S  R.  VLF 
transmitters  have  been  used  at  different  locations  for  the  study  of  wave-particle 
interactions  in  the  magnetosphere.  However,  it  may  not  be  possible  to  use  the 
existing  HF  (or  MF)  and  VLF  transmitters  for  the  proposed  experiments.  Since 
VLF  waves  in  a  narrow  frequency  band  can  be  generated  by  lightning  storms,  the 
experiments  can  be  carried  out  with  available  ionospheric  heating  facilities  during 
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lightning  xtornis  it  is  not  nccexsutx  ts*  nttxc  hghtnine  storms  occuinni*  ttearbx.  li 
lightning  occurs  tit  the  conjugate  location  in  the  opposite  hemisphere,  the  Inzhtnintz- 
induced  \  L.  F  xx  to  es  can  propagate  in  si  ucted  xx  histier  w  ax  s'  lit  sides  and  bounce  hack 
and  forth  lor  sex  end  cxcles  before  thex  die  out. 
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Air  f  orce  Contract  to  Imcoln  1  ah  =  F I962S-MM. -0(i02i  and  the  R  A  DC  contract  through 
Southeastein  Center  tor  Electrical  Engineering  Education  at  the  Massachusetts  Institute  of 
Technologs,  and  h\  the \SF  grant  \TM-is7l32r  and  the  \FOSR  contract  \FOSR  K8-0127 
at  the  Pols  technic  t  nixersitx  Part  ot  this  \sssrk  was  presented  in  the  international 
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Nearly  monochromatic  signals  at  13.6  kHz  ±1  Hz  injected  from  a  ground-based  VLF  transmitter 
can  experience  a  bandwidth  expansion  as  high  as  1%  (-.  100  Hz)  of  the  incident  wave  frequency 
as  they  traverse  the  ionosphere  and  reach  satellite  altitudes  in  the  range  of  600-3800  km  (Bell 
et  a!.,  1983).  The  off-carrier  components,  having  electrostatic  nature,  aie  believed  to  be  induced 
lower  hybrid  wave  modes.  We  investigate  two  different  source  mechanisms  that  can  potentially 
result  in  the  observed  spectral  broadening  of  injected  monochromatic  VLF  waves.  One  is  the 
nonlinear  scattering  of  VLF  signals  by  ionospheric  density  fluctuations  that  renders  the  nonlinear 
mode  conversion  of  VLF  waves  into  lower  hybrid  waves.  These  electrostatic  modes  result  when 
the  injected  VLF1  waves  are  scattered  by  ionospheric  density  fluctuations  with  scale  lengths  less 
than  \/2  it(c/wpt)(Clc /wo)1 13  ~  2  km  in  the  upper  ionosphere,  where  c,  wpc,  flt,  and  uQ  are  the 
speed  of  light  in  vacuum,  the  plasma  frequency,  the  electron  cyclotron  frequency,  and  the  VLF 
wove  frequency,  respectively.  In  the  absence  of  ionospheric  irregularities,  a  second  mechanism  that 
involves  a  parametric  instability  can  excite  the  lower  hybrid  waves  (Lee  and  Kuo,  1984).  This 
process  tends  to  produce  a  spectrally  broadened  transmitted  pulse  with  peaks  at  a  discrete  set  of 
frequencies  on  both  sides  of  the  nominal  carrier  frequency.  By  contrast,  the  nonlinear  scattering 
mechanism  generates  single-peaked  spectra  centered  at  the  carrier  frequency.  Both  types  of  spec¬ 
tra  were  observed  in  the  experiments.  Therefore,  the  two  suggested  source  mechrnisins  contribute 
addiciuely  to  the  observed  spectral  broadening  of  injected  VLF  waves. 


1.  Introduction 

Nearly  monochromatic  signals  at  13.6  kHz  1  1  Hz  in¬ 
jected  from  a  ground-based  VLF  transmitter  were  recently 
observed  to  experience  bandwidth  expansion  as  they  tra¬ 
versed  the  ionosphere  and  reached  the  satellite  altitudes 
in  the  range  of  600-3800  km  [ Bell  et  al.,  1983],  This  ex¬ 
pansion  of  bandwidth,  which  results  in  a  proportional  in¬ 
crease  in  signal-to-noise  ratio,  may  be  as  large  as  1%  (~ 
100  Hz)  of  the  carrier  frequency,  and  the  off-carrier  compo¬ 
nents  are  thought  to  be  electrostatic  in  nature  (/nan  and 
Bell,  1985],  They  occur  only  in  the  presence  of  impul¬ 
sive  VLF  hiss  and/or  a  lower  hybrid  resonance  (LHR)  noise 
band  with  an  irregular  cutoff  frequency,  and  only  for  signals 
whose  frequencies  exceed  the  LHR  frequency  at  the  satel¬ 
lite  location.  Based  on  these  observed  characteristics  of  the 
spectrally  broadened  transmitter  signals,  one  may  suspect 
the  off-carrier  components  to  be  either  electrostatic  whistler 
wave  modes  or  lower  hybrid  wave  modes. 

The  linear  scattering  source  mechanism  producing  quasi- 
eiectrostatic  whistler  wave  modes  was  first  suggested  by 
Bell  ei  al.  [1983],  who  hypothesized  the  creation  of  the 
required  ionospheric  density  fluctuations  by  precipitating, 
low-energy  (<  1  keV)  electrons.  Such  precipitation  events 
have  been  accompanied  by  both  VLF  hiss  and  irregular 
LHR  noise  bands  [McEwen  and  Barrington,  1967;  Laaspere 
et  al.,  1971;  Gumctt  and  Frank,  1972].  Bell  et  al.  [1983] 
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then  speculate  that  the  broadening  of  the  transmitted  pulse 
spectrum  results  from  the  scattering  of  the  initial  signals 
from  the  precipitation-induced  density  fluctuations  and  the 
subsequent  coupling  into  quasi-electrostatic  whistler  wave 
modes  of  short  wavelength.  The  Doppler  shift  associated 
with  these  short-wavelength  modes  is  estimated  to  be  large 
enough  to  produce  the  bandwidth  expansion  of  the  signals 
measured  on  a  moving  satellite. 

However,  we  interpret  the  off-carrier  components  of  the 
spectrally  broadened  transmitter  signals  to  be  lower  hybrid 
wave  modes.  The  main  reason  is  that  the  lower  hybrid 
wave  modes  (eigenmodes  of  ionospheric  plasmas)  experience 
less  spatial  attenuation  than  the  quasi-electrostatic  whistler 
wave  modes  (quasi-modes  of  ionospheric  plasmas)  during 
propagation  from  the  source  regions  to  the  satellite  loca¬ 
tion.  Two  source  mechanisms  are  proposed,  whereby  the 
transmitted  VLF  waves  can  potentially  produce  the  lower 
hybrid  waves.  One  mechanism  comparable  to  that  of  Bell  et 
al.  [1983],  termed  “nonlinear  scattering  process,”  occurs  in 
the  presence  of  ionospheric  density  irregularities.  The  other 
one,  by  contrast,  does  not  require  the  presence  of  ionospheric 
irregularities.  It  is  a  parametric  instability  that,  excited  by 
intense  enough  VLF  waves,  can  generate  lower  hybrid  waves 
and  ionospheric  irregularities  concomitantly  [Lee  and  Kuo, 
1984], 

The  first  mechanism,  presented  in  section  2,  differs  from 
that  of  Bell  et  al.  The  newly  proposed  nonlinear  process  is  a 
scattering  of  the  whistler  mode  VLF  signals  by  preexisting 
density  fluctuations  that  render  a  mode  conversion  to  lower 
hybrid  waves,  instead  of  a  linear  process  that  couples 
energy  into  the  quasi-electrostatic  branch  of  the -whistler 
mode.  We  formulate  a. theory  of  VLF  wave  scattering  off 
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ionospheric,  irregularities,  showing  how  this  process  depends 
upon  the  local  fractional  density  nucrua'-or.- .  ‘.h>-  r?t;r. 
of  their  scale  lengths  to  the  VLF  wavelength,  etc.  The 
scattering  of  VLF  waves  by  ionospheric  density  fluctuations 
generally  causes  elliptically  polarized  modes.  The  induced 
elliptically  polarized  modes,  however,  may  be  predominantly 
electrostatic  under  certain  ionospheric  conditions. 

The  second  mechanism,  described  in  section  3,  produces 
a  spectrally  broadened  signal  via  a  parametric  instability 
[Lee  and  Kuo ,  1984].  The  injected  VLF  wave,  if  intense 
enough,  parametrically  excites  a  magnetic  field-aligned 
purely  growing  mode  together  with  two  sidebands  of  lower 
hybrid  waves.  The  parametric  excitation  process  requires 
the  frequency  and  wave  vector  matching  relations  to  be 
satisfied.  The  two  lower  hybrid  wave  sidebands  are 
characterized  by  their  wave  vectors  which  are  the  sum  and 
difference  of  the  wave  vectors  of  the  whistler  pump  and  the 
purely  growing  mode,  respectively. 

The  lower  hybrid  waves  generated  by  either  of  the  two 
proposed  mechanisms  can  have  wavelengths  much  shorter 
than  that  of  the  incident  VLF  waves.  In  other  words,  much 
greater  wave  numbers  are  introduced  by  the  suggested  source 
mechanisms  for  the  injected  VLF  waves  that,  essentially, 
become  electrostatic  modes.  These  electrostatic  plasma 
modes  with  large  wave  numbers,  detected  by  satellites 
moving  across  the  geomagnetic  field,  give  rise  to  Doppler 
shifts  that  were  observed  as  the  broadened  spectra  (i.e., 
off-carrier  components)  of  the  injected  VLF  waves.  The 
preliminary  results  were  presented  in  Lee  et  a'l.  (1986). 

2.  Nonlinear  Wave  Scattering 


Theory 

A  monochromatic  VLF  wave  transmitted  from  a  ground- 
based  station  into  space  has  been  observed  to  change  from 
linear  into  circular  polarization  f i.e.,  whistler  mode).  If 
a  ducted  whistler  wave  mode,  propagating  exactly  along 
the  geomagnetic  fieid,  is  considered  for  simplicity,  the  wave 
electric  field  may  be  represented  as 

E  =  Eo{x  +  iy)exp(5(fc02  -  w0t)j  (1) 

where  the  z  axis  has  been  taken  along  the  geomagnetic  field 
lines;  u>„  is  the  transmitted  wave  frequency,  and  ko  the 
associated  wave  vector,  assumed  to  be  along  the  z  axis  as 
illustrated  in  Figure  la.  Propagating  into  the  ionosphere 
with  no  irregular  structure,  the  whistler  wave  satisfies  the 
electromagnetic  wave  equation 
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Fig.  1.  Geometry  of  wave  propagation  for  (a)  uniform 
ionosphere  and  (6)  ionosphere  with  irregularities. 
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where  j0  is  the  uniform  oscillatory  current  driven  by  the 
incident  wave  field  in  the  ionospheric  plasma, 


•  _  ie3NeE0 
Jo  ~  m(u,0  -  n,) 

and  e,  Na,  m,  and  0»  are  the  electron  charge,  uniform 
background  plasma  density,  electron  mass,  and  the  unsigned 
electron  cyclotron  frequency,  respectively.  In  the  presence  of 
field-aligned  ionospheric  density  irregularities,  the  scattered 
wave,  allowing  for  electrostatic  modes,  can  be  described  by 

V(V  ■*.)-*>*.  + (3) 

where  6j,  is  the  nonuniform  current  induced  by  the 
interaction  of  the  whistler  wave  field  with  the  density 
irregularities.  We  take  the  x  axis  along  the  direction  of  the 
dens.ty  irregularities  which  are  assumed  to  be  of  the  form 


Sn  =  £nexp[iA:xj 


<4) 


depicted  graphically  in  Figure  16.  The  nonuniform  current 
6j,  can  then  be  expressed  as 

<5j,  =  —e(N06v  +  ~j0)  (5) 

where  6v  is  the  induced  velocity  perturbation.  The  scattered 
wave  field  has  the  general  form  of 

E,  =  [xBt  +  iyjF,]exp[i(fc0z  -  u0t))  (6) 

Solving  the  wave  equation  together  with  the  electron 
momentum  equation,  the  scattered  wave  field  is  found  to 
be  elliptically  polarized, 

E,  =  E0(6n/Na)[x  -f-  (x  +  iy)a]exp[i(k0z  -  ioat) 1  (7) 


where  a  is  proportional  to  the  amplitude  of  the  circularly 
polarized  component  of  the  wave  defined  by 


2upCui0 


ksc3(a>0  -  ne) 


(8) 


where  and  c  are  the  electron  plasma  frequency  and 
the  speed  of  light  in  vacuum,  respectively,  and  k  is  the 
wave  number  of  the  density  perturbation  denned  by  (4). 
Rewriting  the  scattered  field  in  terms  of  the  ionospheric 
irregularity  scale  length,  A  =  2rr/k,  we  find 

E,  =  Eo  exp\ikx)(6h/ No)  exp[:(£0z  -  A>„t)j 

,\2 

x  iy)^|J  (9) 


LP 


CP 


where 

A.  =  \/27r(c/svp«)(fl,/ai0)1/:!  s: 4.4/fc0  (10) 

The  scattered  field  is  composed  of  a  linearly  polarized  (LP) 
part  and  a  circularly  polarized  (CP)  part  as  indicated  in  (9). 
If  the  ratio  (A2/A2)  «  1  (i.e.,  |k|  »  |ko|),  the  scattered  wave 
is  dominated  by  the  linearly  polarized  component  whicl}, 
corresponding  to  a  lower  hybrid  mode,  oscillates  in  the 
direction  of  the  irregularity  wave  vector  k;  hence,  the  wave 
is  electrostatic  in  nature  and  exhibits  a  broadened,  enlarged 
wave  vector  spectrum  relative  to  the  incident  wave.  The 
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Fig.  2.  Schematic  of  Doppler  shift  phenomena. 

Doppier  shift  frequency  measured  by  a  moving  satellite  will 
be  enlarged  as  elaborated  in  the  next  section. 

Analysis 

The  condition  found  above  for  effective  nonlinear  scatter¬ 
ing  into  electrostatic  modes  requires  A  -C  Ac.  Typical  plasma 
parameters  for  the  upper  ionosphere/lower  magnetosphere 
region  of  interest  are  uP,/2r  —  0.7  MHz  and  Qe/2~  —  0.6 
MHz.  Assuming  an  incident  frequency  of-70/2r  =  13.6  kHz, 
we  calculate  Ae  =;  2  km.  The  magnitude  of  the  bandwidth 
expansion  of  the  scattered  VLF  wave  measured  by  a  moving 
satellite  due  to  the  Doppler  shift  effect  is  given  by 

w  k*  v»  .  v>c°*e  „  v- 

-*T 

where  v,  is  the  velocity  of  a  satellite  moving  across  the 
Earth’s  magnetic  field  as  shown  in  Figure  2,  and  k,  = 
kx-rkvi.  Given  v,  >.  8  kni/s  and  A/  ~  100  IIz,  we  find  A  ~ 
80  m  <£  Ac  ~  2  km,  so  the  condition  for  nonlinear  scattering 
is  well  satisfied.  Thus  we  have  found  that  ionospheric 
irregularities  with  scale  lengths  of  several  tens  of  meters 
can  produce  Doppler  shifts  of  the  observed  magnitude. 
Since  k( A)  »  (<)fc<,(A0),  the  simplification  made  in  the 
formulation  of  the  theory,  namely,  the  assumption  that 
the  ducted  whistlers  propagate  exactly  along  the  Earth’s 
magnetic  field,  is  justified.  Under  this  condition,  even  if 
the  ducted  whistlers  have  nonzero  wave  normal  angles,  the 
formulation  of  nonUnear  scattering  mechanism  will  not  be 
significantly  modified. 

Ionospheric  irregularities  can  occur  naturally.  There 
are  many  mechanisms  that  can  generate  field-aligned  iono 
spheric  irregularities  [e.g.,  Fejer  and  Kelley,  1980].  Particle 
precipitation  invoked  by  Bell  et  al  [1983/.  can  produce 
irregularities  with  a  broad  range  of  scale  lengths  in  a  region 
near  the  F  layer  peak,  about  300  km  above  the  Earth’s 
surface.  The  F  region  irregularities,  as  Bell  et  al.  suggested, 
provide  a  favorable  ionospheric  condition  for  the  linear 
scattering  of  the  VLF  waves  and  the  subsequent  coupling 
into  quasi-electrostatic  whistler  wave, modes.  However,  the 
amplitude  of  such  a  quasi-mode  wave  decreases  with  distance 
as  it  propagates  upward  tc  the  satellite  altitudes  of  600- 
380(1  km.  This  arises  from  the  fact  that  the  electrostatic 


whistlers  are  not  eigenmodes  but  quasi-modes  of  ionospheric 
plasmas.  Therefore,  during  the  ionospheric  propagation 
of  quasi-electrostatic  whistlers,  the  ionospheric  plasmas  do 
not  respond  in  phase  to  these  quasi-modes.  This  phase 
mixing  (or  phase  mismatching)  can  attenuate  the  waves 
significantly  during  their  propagation  from  the  source  region 
to  the  satellite  location. 

By  contrast,  our  pioposed  nonlinear  scattering  mechanism 
can  produce  lower  hybrid  waves  which  are  eigenmodes  of 
ionospheric  plasmas,  experiencing  less  spatial  attenuation 
than  the  quasi-electrostatic  whistler  modes  during  propaga¬ 
tion.  This  mechanism  is  effective  provided  that  field-aligned 
ionospheric  irregularities  exist  and  that  the  VLF  wave  fre¬ 
quency  is  greater  than  the  local  LHR  (lower  hybrid  reso¬ 
nance)  frequency.  The  characteristics  of  the  nonlinear  scat¬ 
tering  mechanism,  then,  explain  the  observed  facts  that  the 
spectral  broadening  phenomenon  occurs  only  in  the  pres¬ 
ence  of  impulsive  VLF  hiss  and/or  a  LHR  noise  band  with 
an  irregular  cutoff  frequency,  and  only  for  signals  whose  f-e 
quencies  exceed  the  LHR  frequency  at  the  satellite  loca’.  ■.-■n. 

3.  Parametric  Excitation  of  Lower 
Hybrid  Waves 

In  this  section  we  discuss  the  proposed  second  mechanism 
responsible  for  the  observed  VLF  spectral  b-oadening. 
While  the  nonlinear  scattering  mechanism  does  .lot  work  in 
the  absence  of  ionospheric  irregularities,  lower  hybrid  wave 
modes  can  be  excited  parametrically  by  the  VLF  waves  in 
the  ionosphere.  The  parametric  excitation  of  lower  hybrid 
waves  by  the  injected  VLF  waves  has  been  suggested  by  us 
[Lee  and  Kuo,  1984']  with  the  following  scenario. 

Propagating  from  the  neutral  atmosphere  into  the  iono¬ 
sphere,  the  incident  VLF  wave  can  be  considered  to  be  a 
ducted  whistler  mode  with  k„  =  k0i  along  the  geomagnetic 
field  lines,  taken  to  be  the  z  axis  of  a  Cartesian  system  of  co¬ 
ordinates.  Assuming  the  space-time  dependence  of  the  per¬ 
turbations  to  be  exp[i(k  •  r  -  sat)],  the  proposed  parametric 
instability  can  be  described  by  the  following  wave  frequency 
and  wave  vector  matching  relations: 

L’l-g.  —  ui,  —  v  (12a) 

ki-j  -  xk,  =  zk0  =  ki_  4-  xk,  (126) 

where  the  subscripts  minus/plus  and  s  represent  the 
Stokes/anli-Stokes  components  of  the  high-frequency  lower 
hybrid  sidebands  and  the  field-aligned  zero-frequency  mode, 
respectively.  By  the  definition  of  (12a),  the  Stokes  and 
anti-Stokes  components  refer  to  the  lower  hybrid  sidebands 
whose  wave  frequencies  are  downshifted  and  upshifted, 
respectively,  by  w,  with  respect  to  the  whistler  wave 
frequency-**.),,).  For  a  purely  growing  mode  (i.e.,  Re(w,)=0), 
the  Stokes  and  anti-Stokes  components  of  lower  hybrid  waves 
are  characterized  by  their  wave  vectors,  shown  in  (126), 
which  are  the  difference  and  sum  of  the  wave  vectors  of 
the  whistler  and  the  purely  growing  mode,  respectively. 
Both  the  Stokes  and  the  anti-Stokes  components  of  lower 
hybrid  waves  have  to  be  considered  in  this  two-dimensional 
instability  process,  since  ]k,i  3>  jkd  and  the  concomitantly 
excited  low-frequency  mode  is  a  purely  growing  mode.  Thus, 
jklT|  ~  jkj- j  ~  jk,|;  it  is  found  that  the  lower  hybrid  waves 
have  wave  vectors  large  enough  to  yield  the  observed  spectral 
broadening  due  to  the  Doppler  effect. 

It  should  be  pointed  out  that  a  four-wave  rather  than  a 
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simpler  three-wave  interaction  process  [Beiyer  and  Perkins, 
1976;  Biggin  and  Kelley,  1982)  is  proposed  by  Lee  and 
Kuo  [1984].  Wnile  the  proposed  parametric  instability  has 
a  broad  range  of  VLF  (whistler)  wave  frequency  oi0,  the 
three-wave  interaction  process  [Berger  and  Perkins,  1976}  is 
restricted  to  a  very  narrow  frequency  band,  namely,  ioLH  < 
where  u>LH  is  the  lower  hybrid  resonance 
frequency.  Further,  the  excited  low-frequency  mode  in  the 
work  of  Lee  and  Kuo  ( 1 984 j  is  a  zero-frequency  mode.  This 
four- wave  interaction  mechanism  is  found  to  require  much 
lower  thresholds  than  the  three-wave  interaction  mechanism 
to  excite  lower  hybrid  waves  in  the  ionosphere.  Based  on  this 
proposed  four-wave  interaction  process,  we  can  interpret  the 
observations  of  the  VLF  wave-excited  lower  hybrid  waves  in 
the  ionosphere  indicated  in  the  experiments  of  the  Franco- 
Soviet  ARCAD  3  satellite  [ Berthelier  et  al.,  1982]  and  others 
[e.g.,  Chmyrev  et  al.,  1976]. 

At  this  point,  we  should  comment  on  the  effectiveness 
of  the  parametric  instability  in  causing  the  observed  VLF 
spectral  broadening  phenomenon.  Although  the  production 
of  lower  hybrid  waves  by  the  parametric  instability  does 
not  require  the  presence  of  naturally  occurring  ionospheric 
irregularities,  this  mechanism  needs  a  threshold  to  operate. 
The  thresholds  in  terms  of  whistler  (VLF)  wave  electric  fields 
can  be  calculated  from  equation  (226)  of  Lee  and  Kuo  i  1984): 

Em  =  0M(k2,vl/Qt)(m/e)M-112 


•io^VV'2 


where  k,  v,<,  i',,  and  O,  are  the  wave  number  of  field- 
aligned  ionospheric  irregularities,  the  electron  thermal 
velocity,  the  electroti-ion  collision  frequency,  and  the  electron 
cyclotron  frequency;  r;  is  defined  by  [  1  -t- ( M /m)(k0/k)2]/{l  - 
(M/m)(k0/k)a({lt/urt)2]  where  M/m,  k„,  and  u>pc  are 
the  ratio  of  oxygen  ion  (0'r)  mass  to  electron  mass, 
the  wave  number  of  the  incident  whistler  wave,  and  the 
electron  plasma  frequency,  respectively.  It  should  be  noted 
that  the  parametric  wave  number  k  is  determined  by 
the  dispersion  relation  for  the  lower  hybrid  waves,  ua  = 


This  likely  generates  a  two-humped  spectrum  suppressed 
at  the  carrier  frequency  (/<,)  relative  to  the  sidebands 
peaking  at  approximately  f0±  100  K2,  whereas  the  nonlinear 
scattering  mechanism  simply  broadens  the  bandwidth  of  the 
injected  VLF  waves,  producing  single-peaked  spectra.  This 
can  be  clearly  seen  from  equation  (9),  showing  that  the 
scattered  field  intensity  E,  is  the  incident  field  intensity 
E0  multiplied  by  the  fractional  plasma  density  fluctuation 
(6n/N0  <  1).  This  naturally  leads  to  the  broadened 
VLF  spectrum  peaking  at  the  carrier  frequency.  These 
two  types  of  broadened  VLF  spectra  have  indeed  been 
observed.  Nevertheless,  the  single-peaked  spectra  were 
much  more  frequently  seen  than  the  double-humped  spectra 
(H.  G.  James,  personal  communications,  1988).  The 
characteristics  of  our  proposed  two  source  mechanisms  agree 
well  with  the  observations.  Hence,  the  observations  of  VLF 
spectral  broadening  offer  additional  evidence  supporting  the 
excitation  of  lower  hybrid  waves  by  VLF  waves  in  the 
ionosphere  as  suggested  by  Lee  and  Kuo  [1984]. 

4.  Summary  and  Conclusion 

In  summary,  we  have  investigated  two  possible  source 
mechanisms  that  may  be  responsible  for  the  observed 
spectral  broadening  of  injected  VLF  waves.  In  t.he  presence 
of  ionospheric  irregularities  with  various  scale  lengths  as 
short  as  several  tens  of  meters,  the  nonlinear  scattering  of 
the  VLF  waves  off  these  density  irregularities  can  produce 
electrostatic  modes  with  larger  wave  vectors  which  give  rise 
to  the  apparent  spectral  broadening  through  the  Doppler 
shift  observed  by  the  moving  satellites.  The  amplitude  of 
the  electrostatic  modes  depends  linearly  on  the  amplitude 
of  the  ionospheric  irregularities.  The  broadening  produced 
through  this  mechanism  is  of  the  same  magnitude  as 
the  observed  values.  This  mechanism  produces  single- 
peaked  broadened  spectra  of  VLF  waves,  which  were  most 
frequently  observed  in  the  experiments.  There  are  many 
natural  sources  that  can  generate  ionospheric  irregularities 
with  a  broad  range  of  scale  lengths,  providing  favorable 
ionospheric  conditions  for  the  operation  of  the  nonlinear 
scattering  process. 


-j-  [M/m)(k0/k,  )2y12 ,  where  u>LH  is  the  lower  hybrid 
resonance  frequency  defined  by  uLfJ  =  wfl-/(l  4-  u>£t/f23)1/2. 
If  we  choose  the  following  model  parameters:  v,e  =  1.3  x  106 
(m/s),  ijP'/ 2x  =  0.7  MHz  (i.e.,  No  =  6  x  109  m-3),  vc  =  10 
Hz,  =  0.6  MHz,  the  threshold  field  Em  is  found  to  be 

18  pV/m.  According  to  Inan  and  Bell  [1985],  the  measured 
power  flux  of  transmitted  VLF  signals  is  2  x  10-13  W /m3  at 
the  satellite  location;  namely,  the  measured  field  intensity 
of  whistlers  is  12  pV/m,  which  is  of  the  same  order  of 
magnitude  as  the  threshold  of  the  parametric  instability. 
The  measured  power  flux  (field  intensity)  of  whistier  (VLF) 
waves  in  the  lower  ionosphere  is  2  x  10“9  W/m2  (1.2 
mV/m),  which  is  also  of  the  same  order  of  magnitude  as 
the  instability  thresholds  calculated  by  Lee  and  Kuo  [1984J. 
Therefore,  it  is  possible  for  the  instability  to  be  excited 
though  the  VLF  trar.:  Veters  are  operated  in  pulsed  mode 
with  a  typical  duration  of  seconds. 

It  is  interesting  to  note  that  the  proposed  two  source 
mechanisms  are  expected  to  produce  broadened  VLF  spectra 
with  different  shapes.  The  parametric  instability  will 
preferentially  excite  lower  hybrid  waves  with  the  scale 
length  (.X  =  2~jk)  of  ~  90  m  at  the  satellite  location. 


In  the  absence  of  ionospheric  irregularities,  the  spectral 
broadening  of  the  incident  wave  packet  may  also  stem 
from  the  second  mechanism,  which  involves  the  parametric 
excitation  of  lower  hybrid  waves  [Lee  and  Kuo,  1984].  This 
mechanism  requires  a  threshold  that  can  be  exceeded  by 
the  employed  navigation  transmitters.  These  excited  lower 
hybrid  waves  have  much  shorter  (larger)  wavelengths  (wave 
numbers)  than  those  of  the  injected  VLF  waves.  The 
broadening  introduced  by  this  mechan.sm  is  characterized 
by  the  fact  that  it  produce  wave  frequency  spectra  tiiat 
exhibit  suppressed  field  intensity  at  the  carrier  frequency 
and  enhanced  intensity  at  a  discrete  set  of  frequencies  on 
either  side  of  the  carrier  frequency.  Such  spectra  are  less 
frequently  observed  than  those  produced  by  the  nonlinear 
scattering  mechanism.  The  lower  hybrid  waves  excited  by 
this  mechanism  have  enlarged  wave  vectors;  the  apparent 
broadening  due  to  the  Doppler  effect  is  comparable  to  that 
calculated  for  the  first  mechanism. 

We  conclude  that  the  two  suggested  mechanisms  can 
contribute  additively  to  the  observed  spectral  broadening  of 
injected  VLF  waves  reported  by  Bell  et  al  [1983]  and  others. 
This  spectral  broadening  phenomenon  provides  additional 
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evidence  to  support  the  mechanism  suggested  by  Lee  and 
Kuo  [l984j  for  the  excitation  of  lower  hybrid  waves  by 
injected  VLF  (whistler)  wave?. 
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A  convective  instability  that  leads  to  the  filamentation  of  large-amplitude  Alfven  waves  and 
generates  nonoscillatory  spatially  growing  density  fluctuations  is  investigated  by  using  a  two- 
fluid  plasma  model.  Based  on  linear  perturbation  analysis  for  a  four-wave  parametric  coupling 
process,  the  dispersion  relations  of  the  instability  for  both  a  right-  and  left-hand  circularly 
polarized  Alfven  pump  are  derived.  They  are  then  solved  numerically  for  the  dependences  of 
the  threshold  fields  and  growth  rates  on  the  wavenumber  k  of  the  nonoscillatory  mode,  the 
beta  of  the  plasma,  and  the  frequency  of  the  pumps,  for  both  right-  and  left-hand  circularly 
polarized  Alfven  pumps.  The  relevance  of  the  proposed  instability  with  some  observations  in 
the  space  plasmas  is  discussed. 


I.  INTRODUCTION 

Large-amplitude  Alfven  waves  are  frequently  observed 
in  the  high-speed  streams  of  the  solar  wind,1  in  the  upstream 
Jovian  bow  shocks,2  and  in  the  interplanetary  shocks  and  the 
terrestrial  foreshock.3-4  Moreover,  the  presence  of  such 
waves  in  front  of  supernova  shocks  is  believed  to  play  an 
important  role  in  cosmic  ray  acceleration.5-7  Therefore  the 
stability  of  finite  amplitude  Alfven  waves  is  of  great  interest, 
not  only  in  basic  plasma  physics,  but  also  for  their  applica¬ 
tions  to  space  plasma  physics*-"  and  astrophysics. 1  ^ 1 3 

Various  processes  of  parametric  instabilities  excited  by 
the  AlfVcn  pumps  have  been  investigated.  In  general,  circu¬ 
larly  polarized  pumps  propagating  along  the  background 
magnetic  field  are  considered.  Within  the  one-dimensional 
regime,  two  instability  processes  have  been  studied  exten¬ 
sively.  One  leads  to  the  parametric  decay  instability 
(k>  £0),u-15  and  the  other  one  gives  rise  to  the  modulation 
instability  where  fc^  and  £  are  the  wavenum¬ 

bers  of  the  pump  and  decay  mode,  respectively.  In  both  pro¬ 
cesses,  the  decay  modes  and  sidebands  all  propagate  in,  par¬ 
allel  to  the  Alfven  pump. 

The  modulation  instability  has  also  been  analyzed  in  the 
region  that  describes  the  nonlinear  evolution  of  the  Alfven 
waves  propagating  along  a  static  magnetic  field-  It  is  shown 
that  this  evolution  can  be  governed  by  a  single  “derivative 
nonlinear  Schrddinger  equation,”?0-22  or  more  generally,  by 
a  set  of  three  coupled  equations,23  which  in  turn  are  related 
to  the  nonlinear  Schrddinger  equation  with  known  soliton 
solutions.  The  possible  applications  of  Alfven  sohtons  to  so¬ 
lar  and  astrophysics!  plasmas  have  been  discussed  by  Oven- 
den  et  alP  and  Spangler  and  Sheerin.22  c* 

When  the  decay  modes  and  sidebands  are  allowed  to 
propagate  in  multidimensional  space,  new  decay  channels 
for  the  parametric  instabffities  are  also  pwacnt.3*^  Though 
the  effect  of  multidimensional  space  on  the  noolmear  prop¬ 
erty  of  the  Alfven  waves  is  generally  thought  to  be  impor¬ 


tant,  its  understanding  is  still  far  from  satisfactory,  mainly 
because  of  the  complexity  involved  in  the  analysis.  Recently, 
an  absolute  filamentation  instability  process  has  been  con¬ 
sidered  by  Kuo  et  al. 1  *  Based  on  the  two-fluid  plasma  model, 
a  four-wave  parametric  coupling  process  is  analyzed  for  the 
temporal  growing  instabilities,  which  are  excited,  via  the 
scattering  of  the  Alfven  pumps  into  obliquely  propagating 
Alfven  sidebands  by  the  density  perturbations  that  are  asso¬ 
ciated  with  the  simultaneously  excited  purely  growing  and 
magnetic  field-aligned  nonoscillatory  magnetostatic  modes. 
The  excited  Alfven  sidebands  are  then  added  to  the  pumps 
and  lead  to  the  filamentation  of  the  AlfVcn  pumps.  The  non¬ 
linear  source  for  the  Alfven  sidebands  is  the  beating  current 
driven  by  the  pump  wave  fields  on  the  density  perturbation 
of  the  purely  growing  mode,  whereas  the  driving  forces  for 
the  purely  growing  nonoscillatory  density  perturbations  arc 
the  nonlinear  Loren tz  forces,  which  reduce  to  the  pondero- 
motive  forces  in  the  unmagnetized  case  and  are  perpendicu¬ 
lar  to  the  background  magnetic  field 

In  this  paper,  a  convective  filamentation  instability  pro¬ 
cess,27  in  contrast  to  the  previous  work,"  is  studied.  Since 
the  excited  instabilities  grow  in  space  ( not  in  time)  along  the 
propagation  direction  of  the  pumps,  the  nonoscillatory  den¬ 
sity  perturbations  are  not  magnetically  field  aligned  any¬ 
more.  Consequently,  the  nonlinear  Loren  tz  forces  intro¬ 
duced  by  the  spatial  gradient  of  the  resultant  Alfv£n  wave 
intensity  also  become  not  field  aligned;  The  axial  compo¬ 
nents  of  the  forces  can  drive  density  perturbations  more  ef¬ 
fectively;  than  the  transverse,  components,  which  haw  to 
move  charged  particles  across  the  magnetic  fidd  lines.  From 
this  aspect,  the  convective  instability  process  considered  in 
this  work  is  expected  to  be  characteristically  different  from 
the  absolute  filamentation  usability  process  presented  in 
the  previous  work."  ..  -  A 

In  Sec.  II  the  coupling  equations  for  Alfven  sidebands 
and  the  puidy  growing  inode  arc  derived  A  dispenioQ  rela¬ 
tion  is  obtained  in  Sec;  III.lt  is  analyzed  for  the  threshokl 


fields  and  the  growth  rates  for  various  cases.  The  numerical 
results  are  also  presented  in  Sec.  III.  Finally,  presented  in 
Sec.  IV  are  the  summary  and  discussion. 

II.  COUPLING  EQUATIONS 

Monochromatic  circularly  polarized  Alfven  waves 
propagating  alon»  a  constant  magnetic  field  Bu  in  z  direc¬ 
tions  are  represented  by 


Bp±  =(x±  ,y)Bp ±  e*** ' +  c.c. ,  ( 1 ) 

where  Bp  ±  are  the  unperturbed  wave  magnetic  field  intensi¬ 
ty,  assumed  to  be  constant  and  real  quantities  for  simplicity; 
the  ±  signs  refer  to  the  right-  and  left-hand  circular  polari¬ 
zations,  respectively;  kg  and  co0  are  the  wavenumber  and  an¬ 
gular  wave  frequency,  satisfying  the  dispersion  relation  a>l 
=  (1  ±£>0/ft(  )£o*d>  wherein  oA  and  ft,  are  the  Alfven 
speed  and  the  ion  cyclotron  frequency,  respectively. 

The  quiver  velocities  of  electrons  and  ions  introduced  by 
such  wave  fields  can  be  expressed  as 


r±  ~(1±o’)Vd 


=  —  (2) 
*o  B0 

where  the  subscripts  e  and  /  refer  to  electrons  and  ions. 

In  the  following  analysis  of  nonlinear  wave-plasma  in¬ 
teraction,  the  plasma  is  assumed  to  be  collisionless  and  spa¬ 
tially  uniform.  The  process  under  investigation  is  the  con¬ 
vective  filamentation  instability  of  Alfven  pumps.  The 
pumps  are  scattered  into  obliquely  propagating  sidebands  by 
the  simultaneously  induced  nonoscillatoiy  density  perturba¬ 
tions  responding  to  the  presence  of  the  nonuniform  radiation 
pressure  of  total  Alfven  wave  field.  It  is  a  linearized  version 
of  the  self-focusing  process  for  a  large-amplitude  and  Urge 
cross-sectional  wave  beam  propagating  in  a  nonlinear  medi¬ 
um. 

The  Alfven  sidebands  are  excited  through  the  beating 
current  density  driven  by  the  Alfven  pump  wave  fields  in  the 
density  perturbation  /»„  induced  by  the  self-consistent  radi¬ 
ation  pressure.  The  coupled  mode  equations  for  these  side¬ 
bands  can  be  derived  from  the  following  fluid  equations  for 
electrons  and  ions  together  with  the  Maxwell  equations: 

J^«±  +V-(n05rri  +n,rpt±) 

=  0~4^n±  +  («o 5v.-*  +«*Vt  )•  <3> 

at 

4-SB,  =B0Vx(8rri  xz>,  (4) 

at 

—  5v,±  =  ~c?V-^-  +  ft,(5vt*  -5vr<  )Xz,  (5) 
dt  «o 


Vx5B  t  =  (4an0e/c)  [  (5v, ,  -  SvTt  ) 

+  («y»o)<Vi  >]»  (6) 

where  c,  =  1(7^  -I-  T, )/m,  J 1/2  is  the k»  acoustic  speed. 


spatial  growth  rate  along  thezaxis,  and  (k,a^uu )  represent 
the  lower  and  upper  sidebands.  In  general,  the  upper  side¬ 
band  (k2  =  lq,  +  k  —  ha,  co2  =  co0)  and  lower  sideband 
(k,  =  k,,  —  k  —  bcz,  co,  =  to0)  propagate  together  along  the 
magnetic  field  and  form  a  standing  wave  pattern  across  the 
magnetic  field.  Consequently,  the  self-consistent  radiation 
pressure  and  the  induced  nonoscillatory  density  perturba¬ 
tion  should  also  retain  a  similar  standing  wave  pattern  across 
the  magnetic  field.  Without  losing  the  generality,  we  may 
first  prescribe  the  functional  form  of  the  nonoscillatory  den¬ 
sity  perturbation  to  be 

n,  —  h,e**coskx ,  (7) 

where  h,  is  a  real  amplitude;  quasineutrality  has  been  as¬ 
sumed.  Following  such  a  functional  dependence  for  the  non¬ 
oscillatory  density  perturbation,  the  physical  quantities  of 
the  sidebands  can  be  expressed  accordingly  as 

=  [(x8vri±I  ±iy8V'lty)coskx 


+  iz  <5D„  ± ,  sin  kx]eTe 


<(*b*  -  "o'! 


+  C.C.  , 


Sn  t  —  i8h,  sin  lexeme <tr*  +  c.c. , 

5B.  =  [(i  ±  ia  ^  y)cos  kx  —  iz(k /kg)sin  kx] 


X8B  t  eTe 


+  C.C. , 


where  a.  are  unknown  parameters  defining  the  polariza¬ 
tion  of  the  Alfven  sidebands  and  will  be  determined  shortly; 
the  notation  k„  =  k(l~  he  is  used  to  simplify  the  expression. 

Substituting  (7)  and  (8)  into  (3)-(6),  and  eliminating 
5h  x  from  (3)  and  (5),  we  first  obtain 

Svrtx  —  -  )SB t  /Bq  , 

-  a*  <,a0/kv)8B t  /B0  , 

Shi t,  =  [(kn&/coa)(hJna)(Bp±  /B0) 

—  —  1) 

XOTe,  (o0/SlMSB^  /B0 )  ]//, 

Si>i±r  —  +  (O,/«0 ){(^o^/t?0 )(n,/n0 ) {Bp±  /B0) 

+  (cu0/ko)[l  +  (ko/k)\Sli/a0'i 
X{{ol/klct%  —  \)(±a±  -Oo/ft,)] 
X(SB^/B0)}/f, 

Svi±t  =  -  (ko/kHQSa b)*[(I 

X  {koV^/ao ) (n,/«o )(Br±  /B0 ) 

+  (ao/koHi  o0/alHSB±  /B^\/fv 
&or±I  =  Svi± ,  ± a±  , 

where  /=  1  +  (ko/kj1(£li/a0)1{  I  -  ruJ/QJ) 

-D-  -■  -*  • 


equatioos-Incorporating  the  relations  of  (?)  tn.thase  two 


Let  k = £&  be.  the  filamesutum  wave  vector,  *  be  the  nuneter 
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«±  -{(1  ± [ 1  ~  (ko/kfial/klc1;  -  1)] 

+  (/±Sl//«0)(ftS»> a /^?) 

X  ( 1+  Ar  V*g  >}/{(  l,r±  ) 

X[1  ~  (k0/k)2{a>l/k^  ~  1)]  +  (1  -Stf/colf) 

X  [  fklv\/a *  -  (ko/kiHri/ktf  -  1)]}  (10) 

and  that  of  a  single  mode  equation  for  the  combined  Alfven 
sidebands 


{1  -  (%0/k)2(tl,/a0)2(ly:a±a0/(li)(oZ/kZc*-  1)/ 
+  «±  (k^vi/il,<v0)}(SB ±  /B0) 

=  -  (Vf±a0/(li)ik0k0olA/col) 
X(h,/n0)(Bp±/B0).  (11) 

Equation  (11)  shows  that  the  sideband  fields  are  pro¬ 
duced  by  scattering  the  pump  off  the  density  perturbation 
induced  by  the  radiation  pressure  force.  From  the  z  compo¬ 
nent  of  the  one-fluid  momentum  transfer  equation,  the  den¬ 
sity  perturbation  induced  by  the  radiation  pressure  force  is 
governed  by 


^  d_  nj__  _ 

‘  dz  n0  M, 


(12) 


where  F„  and  Fit  are  the  radiation  pressure  forces  on  elec¬ 
trons  and  ions,  respectively,  and  are  given  by 


d 


X(z- v,±  X<5v*  ±  )  -f  c.c. , 


(13) 


This  is  a  general  expression  of  the  dispersion  relation  for 
the  convective  filamentation  instability  of  circularly  polar¬ 
ized  Alfven  waves.  It  is  analyzed  in  the  following  to  deter¬ 
mine  the  threshold  conditions  and  the  growth  rate  of  the 
instability. 

We  first  determine  the  threshold  conditions.  Let  *  =  0, 
i.e.,  k0  =  kc.  Equation  (15)  is  then  solved  to  obtain  the 
threshold  intensity  to  be 

(Bp±/B0)2h  =  ~(P/2)A±/(B±  +2A±),  (16) 

where  f3  —  c\/v\  is  the  ratio  of  the  plasma  kinetic  pressure 
to  the  background  magnetic  pressure; 

A±=  1  -  (V*)2(n,A>0)2(l  :fa0±a>0/ft,) 

XUol/klc*  -  1  )/f0  :f  a0±  (klv\/O.,(o0)-, 

Bk  -  d//o  ±*>o/toi)(klvi/(t>l)(a0±  +  i+k2/k2); 

a0±  =«±k-»b  an£l  fo  —fllo-ko  ■ 

From  (16),  A±  /(B±  +  2A±  )<0  is  the  additional  crite¬ 
rion  of  the  instability.  This  criterion  limits  the  regions  of  k  of 
the  instability.  It  can  be  manifested  by  a  simple  example. 
Considering  a  special  case  that  /?=  (1  +  Oi0/i\l ),  it  sets 
co\/k\c\  —  1=0.  Thus,^  =  1, 

«*b±  =  ± i(ft,/*>0 )[1  ±  («o/fl, ) (2  +  k*/kl)\ , 

A±  =  (l/2£)[l  +  (o>0 m,)(.k2/kl)} , 
and 

B±  =  (2  +  3k2/2kl)  ±  J(ft,/<»0 ) . 

Now,  we  have  to  analyze  separately  the  criterion  for  each 
pump  polarization. 


where  etJ  =  +  1. 

Substituting  (2),  (7),  (8),  and  (9)  into  (13),  yields 


«o 


+  1  + 


k2  \8B\ 

kfJBT 


[c^/ui,  +  4(Bp±/B0)2]  ' 


(14) 


This  equation  gives  the  nonoscillatory  density  perturba¬ 
tions  as  a  function  of  the  pump  and  sidebands.  Thus  Eqs. 
(11)  and  ( 14)  form  a  complete  coupled  set  of  equations 
describing  the  convective  filamentation  instability  of  circu¬ 
larly  polarized  Alfven  waves. 


HI.  DISPERSION  RELATION 

Substituting  (14)  into  (If),  one  equation  relating  SB  ± 
to  its  complex  conjugate  SB*  is  obtained.  This  equation, 
together  with  its  complex  conjugate,  leads  to  a  dispersion 
relation  as 

|{1  -  (VA)a«V«b)2(l 

X(a>l/k2<*  -l)/f^a± 

+  (W±  (klo2A/a^)(a  ±  +  l+k2/k%) 

X  {3,±  /B0 )2/[c}/ £  +  4(B,±  /B0 )2]}j2 

=  Kp//±c«o/ft/)(ifcolM/®o)(al  +  1  +  AtV*?1) 

X(Br±  /Bof/tf/v2  +4(5,±/^)2]}j2.  (15) 


A.  Right-hand  (RH)  circular  polarization  (+) 

Here  B+  +  2A+  =  (2  +  \/p  +  jft,/ft>o)  +  -  ( 1/ 

B)(,o)0/(li)]{k2/kly,  since  0>l  and  co o/ft,<l, 
B+  +2d+>0  for  all  k.  Hence  the  criterion  reduced  to 
A + <,0,  which  leads  to  the  condition  k  2/k  l  for  the 

instability. 

B.  Left-hand  (LH)  circular  polarization  (— ) 

Since  A  _  >  0,  it  requires  that 

+2A.  =  [2-i(fi,-/<y0)  +  l/iff] 

+  [1+  (l/0)(fi>o/ni)](A'/*3)<O. 
The  region  of  k  of  the  instability  is  then  derived  to  be 

k2/k20  <  [j(ft,/a>0>  -2-  , 

providing  that  Oq/SI,  <  2/(7  +  V33 )  — 0. 157. 

For  the  general  case,  i.e.,  arbitrary  values  of  ft,  the  un¬ 
stable  region  of  k  imposed  by  the  criterion  can  be  determined 
numerically. 

In  the  RH  case,  A + (kJ},co0)  has  a  zero  at  k  =  £*(/?>&<>), 
and  d+<0  for  k>k2,  and  A + >0  for  k<k2.  It  is  shown  that 
B+  -f  2^+>0  for  all  k  in  the  region  p>  1  +  Oq/SI,.  How¬ 
ever,  B++2A+  also  has  a  zero  at  k  =  kx(p,o0)  for 
1  +  coo/ilt  and  thus,  B+  +2A+  becomes  negative  in  the 
region  k<kx.  The  dependences  of  A,  and  ^  on#  for  two 
cases  ofS>0  — 0.3  and  0.01  are  presented  in  Fig.  i,  where 
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FIG.  1.  Boundary  curves  A,(/9,5„)  and  k,((3,a>„)  divide  k-f)  space  into 
three  legions,  where  k  =  k  /k„.  Instability  is  forbidden  in  the  region  between 
two  curves.  Two  cases  for  5>„  =  0.3  and  0.01  are  presented  for  the  RH  circu¬ 
larly  polarized  pump. 


=  ky/k^,  k2  —  k^/ko,  and  a0  =  oyTI,  are  normalized 
quantities.  These  two  curves  k,  and  k2  divide  the  k-0  space 
into  three  regions  and  A+/(B+  +  2A  +  )  >0  in  the  central 
region.  Therefore  there  are  two  instability  regions  Q<k  <  kx 
and  k>k2  for  0<  1  +  co0,  and  only  one  instability  region 
k>k2  for  1  -f  co0.  In  other  words,  both  large  and  small 
scale  filamentation  instability  can  be  excited  by  the  RH  Alf- 
ven  pump  in  low  beta  plasma,  while  only  small  scale  instabil¬ 
ity  exists  in  very  high  beta  plasma. 

In  the  LH  case,  A  _  is  always  positive.  However,  the  sign 
of  B_  +  2A_  depends  on  k,  0,  and  <50.  For  large  co0  (e.g., 


it 


FIG.  3.  The  dependence  of  the  threshold  field  on  wavenumber  k  and  plasma 
P  for  RH  pump  having  frequency  cD0  —  0.3. 


50  =  0. 3 ) ,  B  _  +  2A  _  has  only  a  single  zero  at  k  —  k ,  (/?,to0 ) 
for 0<0X  and  +  2A_  <0  for  0<k<kx;  asw0  becomes 
very  small  (e.g.,  tyo  =  0.01),  two  zeros  of  B_  +  2/4  _  are 
present  at  k  =  k3  and  A:„  respectively,  where  k3  <  kx  and  k3 
exists  only  for  0  larger  than  a  constant  value  03.  Again, 
B_  +  2 A_  <0  for  (0,k3)  <k<kx.  These  boundary  curves 
defined  by  k,  and  k3  are  shown  in  Fig.  2  for  two  cases  of 
&>0  =  0.01  and  0.3.  It  shows  that,  for  a>0  =  0.01,  kx  exists  in 
the  entire 0_ region  and  k3  exists  only  for0>  0.99  =  03\  how¬ 
ever,  only  ky  exists  in  the  region  0  <  0.7  for  o)Q  =  0.3.  In  any 
event,  instability  exists  only  in  the  region  below  the  kx  curve 
and  above  the  k3  curve  or  the  zero  line. 

Presented  in  Figs.  3-6  are  the  dependence  of  the  normal¬ 
ized  threshold  field  r/th  —  (&P±  /B0)lh,  defined  by  (16)  on 


FIG.  2  For  the  LH  pump  case,  die  instability  region  in  A:-/? space  is  confined . 
between  either*:, 09.S,,,. and  *,(/?, £>q)  curvesorthe*:,  curve  and  A  ==  0 
line.  ,  -  _  - 


FIG.  4.  Threshold  field  77,,,  vs*  and/?  for  i  very  low  frequency,  =  0.01, 
RH  pump.  The  dependence  in  very  large  *(>  8)  region  is  not  presented. 
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FIG.  5.  Threshold  versus  k  in  four  different  beta  plasmas  for  the  case 
•f  LH  pump  with  frequency  5„  =  0.3. 


FIG.  7.  The  dependence  of  the  spatial  growth  rate  k  on  *  and  0  for  the 
5„  =  0.3  RH  pump  having  intensity  j f  —  0.05,  where  k  = 


k  —  k/kv,  where  k  is  limited  to  the  range  of  2.5  so  that  the 
finite  Larmor  radius  effect  is  still  small  and  thus  does  not 
invalidate  the  fluid  model  of  the  present  work.  Figures  3  and 
4  are  the  results  of  RH  circularly  polarized  pump  for 
<y0  =  0.3  and  0.01  case,  respectively.  In  each  case,  three  val¬ 
ues  of/?  =  0.5, 1.0,  and  1.5  are  considered.  The  results  show 
that  the  threshold  field  increases  with  0  and  k  in  each  insta¬ 
bility  region.  Moreover,  the  threshold  field  increases  very 
fast  with Lin  the  lower  k  region  (i.e.,  0  <*<*, )  and  becomes 
infinite  at  k  = 

Figures  5  and  6  correspond  to  LH  pump  wave.  Again, 
a0  —  0.3  and  0.01  are  used,  respectively,  as  the  representa¬ 
tive  parametric  values.  For  the  high  frequency  pump  case, 
ci0  —  0.3,  instability  can  only  be  excited  in  the  relatively  low 


0  plasma,  0  <0.7.  As  S0  reduces  to  0.0 1 ,  instability  exists  in 
all  0  ranges  of  interest.  The  results  show  that  the  threshold 
field  also  increases  with  0  and  k  and  in  general,  decreases 
with  <y0.  It  becomes  infinite  aik  —  k.  and  k  —  k^  (Fig.  6  for 
0=1  and  1.5  cases). 

We  next  present  the  growth  rate  of  the  instability.  As¬ 
suming  that  the  pump  intensity  tj2  =  (&,±  /B0)2  =  0.05, 
( 15)  is  then  solved  for  the  spatial  growth  rate  k.  Shown  in 
Figs.  7-10  are  the  functional  dependences  of  the  normalized 
growth  rate  k  =  K/k0  on  k  and  0  for  «0  -  0.3  and  0.01  and 
for  RH  and  LH  pump  cases,  respectively.  For  the  RH  pump 
(Figs.  7  and  8),  the  growth  rate  of  small-scale  instability 
{k>  1)  appears  to  increase  with  0,  while  the  growth  rate  of 
large-scale  instability  (k  <  1 )  decreases  with  0.  In  general, 


FIG.  8,  Growth  rate  versus  A:  and  0  for  the  5„ = 0.01  and  tj2  =  0.05  RH 
pump. 
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than  that  of  large-scale  instability.  Figures  9  and  10  corre¬ 
spond  to  the  LH  pump  and  50  =  0.3  and  0.01,  respectively. 


a  peak  in  the  unstable  band.  The  peak  growth  rate  tends  to 
decrease  with 0,  except  in  the  region  near  0= 0.5,  the  peak 
growth  rate  increases  again  and  then  decreases  to  zero  at 
0—0. 7.  For  the  low-frequency  case,  the  dependence  be¬ 
comes  quite  irregular.  This  is  because  there  is  an  extra  for¬ 
bidden  region  0 <k<k3  for  0>03  =  0.99  (for  the 
S0  ss  0.01  case).  In  general,  the  growth  rate  decreases  with 0 
for  0<  0.99  and  increases  with  £  for  0>  0.99.  The  growth 
rate  also  has  a  finite  jump  near  k  =  k3  for  the/?  >  0.99  case. 
This  is  because  the  threshold  approaches  quickly  to  infinite 
at  *  =  k3. 

IV.  SUMMARY  AND  DISCUSSION 

The  convective  filamentation  instabilities  of  large-am¬ 
plitude,  circularly  polarized  Alfven  waves  propagating 
along  the  background  magnetic  field  have  been  studied.  Sim¬ 
ilar  to  the  absolute  filamentation  instability  process,"  the 
convective  instabilities  arc  also  excited  via  the  scattering  of 
the  unperturbed  Alfven  pumps  into  sidebands  bv  the  simul¬ 
taneously  excited  purely  growing  density  perturbations.  The 
nonlinear  source  for  the  Alfven  sidebands  is  also  the  same.  It 
is  the  beating  current  driven  by  the  pump  wave  field  on  the 
density  perturbation  of  the  purely  growing  mode.  However, 
the  characteristic  features  of  the  two  instability  processes  arc 
quite  different.  The  absolute  instability  grows  in  time,  while 
the  convective  instability  grows  spatially  along  the  path  of 
the  pump.  The  purely  growing  density  perturbation  is  asso¬ 
ciated  with  the  magnetostatic  mode  for  the  absolute  instabil¬ 
ity,  it  becomes  an  electrostatic  mode  in  the  convective  case. 
The  nonlinear  forces  generated  to  drive  the  purely  growing 
modes  are  also  different.  The  magnetostatic  mode  is  driven 
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FIG,  9.  Growth  rate  versus  k  and  0  for  the  Sy,  =  0.3  and  if  =»  0.05  LH 
pump.  .  _  1  /  .  r 
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FIG.  10.  Growth  rate  versus  k  and  0  for  the  5„  =  0.01  and  r/2  =  0.05  LH 
pump. 


by.a  transverse  nonlinear  Lorcntz  force,  while  the  electro¬ 
static  mode  is  driven  by  an  axially  directed  nonlinear  Lor- 
entz  force.  Consequently,  the  dispersion  relations  mid  the 
results  of  threshold  field  and  growth  rate  are  all  different  for 
each  instability  process.  Solving  the  dispersion  relation 
(IS),  we  have  determined  the  threshold  fields  and  growth 
rates  of  the  convective  filamentation  instabilities  for  the 
cases  of  the  RH  and  LH  circularly  polarized  pump  waves. 
We  have  also  determined  the  functional  dependences  of  the 
threshold  field  i/th  on  the  wavenumber  k  of  the  nonoscilla- 
tory  mode,  the)?  of  the  plasma,  and  the  frequency  <u0  of  the 
pump  (see  Figs.  3-6).  Similar  dependences  for  the  spatial 
growth  rate  *  =  #r/7co  for  pump  intensity  iy2  =  0.05  have 
been  examined  and  presented  in  Figs.  7-10. 

We  now  apply  the  results  of  the  proposed  instability 
process  to  explain  some  of  the  observations  associated  with 
large-amplitude,  low-frequency  Alfvenic  fluctuations  ap¬ 
pearing  in  the  sclar  wind,  in  the  upstream  regions  of  plan¬ 
etary  bow  shock,  and  in  the  interplanetary  shock.  For  the 
high-frequency  pump  case,  e,g.,  a>0  ~  0.3,  instability  can  be 
excited  by  both  the  RH  and  LH  pump  in  low/?(  <  0.7)  plas¬ 
mas,  however,  only  the  RH  pump  is  unstable  in  the  high 
0{>l)  plasmas.  The. finite  amplitude  AlfVen  waves  that 
were  observed  in  the  high-speed  streams  of  the  solar  wind 
and  discussed  by  Abraham-Shrauner  and  Feldman  are  the 
left-hand,  circularly  polarized  waves  with  5o  —  0.3  and 
p2  =  0.05.  The  results  of  Fig.  9  indicate  that  these  waves  are 
stable  to  the  proposed  instability  for  5  >  0.7  and  become  un¬ 


stable  for/?<0.7.  The  high-frequency,  right-hand  polarized 
waves  have  been  observed  in  the  interplanetary  shocks  and 


are  unstable  in  general  and  the  excited  instability  has  larger 
growth  rate  in  high  beta  plasma  and  larger  *(  >  1)  region. 
Low-frequency,  left-hand  circularly  polarized  Alfven 


bow  shock.  As  shown  in  Fig.  10,  these  waves  are  unstable  in 
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instability  is  also  very  broad.  Though  the  observed  frequen¬ 
cy  is  about  5„  =  0.02,  the  results  of  Fig.  10  for  <y0  =  0.01  are 
qualitatively  the  same  and  hence,  no  additional  figure  is  pre¬ 
sented.  The  proposed  instability  is  also  found  unstable 
against  the  low-frequency  RH  pump  in  the  very  small  k  re¬ 
gion.  As  shown  in  Fig.  8,  the  instability  prefers  to  be  excited 
in  the  low  beta  plasma.  The  existence  of  RH  waves  can  be 
speculated,  because  in  such  a  low-frequency  region  the  dis¬ 
persion  relations  of  right-  and  left-hand  circularly  polarized 
Alfven  waves  are  almost  identical. 

In  general,  the  instability  conditions  for  the  RH  pump 
and  the  LH  pump  are  different,  even  in  the  low-frequency 
region  where  Alfven  modes  become  approximately  degener¬ 
ate.  This  is  because  the  sidebands  in  both  cases  are  elliptical- 
ly  polarized  and  shear  mode  in  nature,  while  the  pumps  are 
circularly  polarized  in  the  opposite  sense.  Thus  the  induced 
nonlinear  Lorentz  forces  on  electrons  and  ions  are  different 
for  differently  polarized  pumps.  Consequently,  different  in¬ 
stability  conditions  are  required  by  differently  polarized 
pumps. 
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Abstract.  Intense  magnetosonic  waves,  originally 
propagating  at  the  right  angle  with  the  interplanetary 
magnetic  field,  can  excite  a  purely  growing  mode  along 
the  interplanetary  magnetic  field  together  with  two 
symmetric  magnetosonic  sidebands  propagating,  obliquely 
across  the  magnetic  field.  This  instability  process  leads 
to  the  filamentation  of  the  magnetosonic  piimp  waves. 
These  two  excited  magnetosonic  sideband  modes  propagate 
together  perpendicularly  across  the  magnetic  field  and, 
meanwhile,  form  a  standing  wave  pattern  along  the 
magnetic  field.  The  thresholds  of  this  filamentation 
instability  can  be  exceeded  in  the  solar  wind  environment. 
It  is  predicted  that  the  density  fluctuations  produced  by 
the  filamentation  instability  along  the  interplane’ -.i-y 
magnetic  field  have  wavelengths  greater  than,  at  least,  a 
few  Earth  radii.  The  polarization  of  the  obliquely 
propagating  magnetosonic  waves  excited  by  the 
filamentation  instability  is  determined  by  the 
characteristics  of  the  magnetosonic  pump  waves  and  the 
environmental  plasmas. 


Introduction 

Recently,  there  has  been  a  considerable  effort  in  the 
study  of  the  stability  of  Alfven  waves  in  the  space 
plasmas.  It  is  attributed  to  the  frequent  appearance,  as 
detected  by  the  space  crafts,  of  large-amplitude  Alfven 
waves  in  the  high-speed  streams  of  the  solar  wind 
[Abraham-Shranner  and  Feldman,-  1977],  in  the  upstream 
Jovian  bow  shocks  [Goldstein  et  al.,  1985],  and  in  the 
interplanetary  shocks  [Vinas  et  al.,  1984],  In  general, 
these  waves  appear  to  be  circular  polarization  and 
propagating  along  the  background  magnetic  field. 
Moreover,  the  ISEE  1  and  2  spacecrafts  have  detected 
both  magnetosonic  waves  and  Alfven  waves  propagating 
upstream  into  the  solar  wind  with  phase  velocities  much 
less  than  the  solar  wind  velocity.  Hence,  these  waves  are 
blown  backward  over  the  spacecraft  and,  consequently,  an 
apparent  -change  in  polarization  of  waves  occurs. 
Intrinsically  rigljt  (left)  hand  polarized  waves  are  seen  to 
have  an  apparent  left  (right)  handed  polarization  in 
spacecraft  records  [Hoppe- and  Russell,  1983]! 

The  parametric  decay  instability  and  modulation 
instability  of  large-amplitude,  circularly  polarized  Alfven 
waves  propagating  in  the  solar  wind  along  the  magnetic 
field  have'  been  analyzed  extensi vely  by  Lashmore-bavies 
[1976],  Goldstein  [1978],  Derby  [1978],  Wong. and  Goldstein 
[1986],  Longtin  and  Sonnenip  [1986],  and  Terasawa  et  al 
[1986].  The  filamentation  instability,  competing  with  these 
instabilities,  can  break  up  the  large-amplitude  Alfven 
waves  and  induce  not  - only  plasma  density  fluctuations  but 
magnetostatic  fluctuations  [Kuo  et  al.,  1988]..  Starting 
from  small  perturbations  in  plasma  density,  the 
filamentation  instability gives;  rise  to  a  modulation  of  the 
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plasma  dielectric  constant -and  then  wave  distribution  that, 
in  turn;. enhances  the  density  perturbations. 

The  modulation  instability  has  also  been  analyzed  in 
the  region  that  describes  the  nonlinear  evolution  of  the 
Alfven  wave  propagating  along  a  background  magnetic 
field  It  is  shown  that  this  evolution  can  be  governed  by 
a  single  "derivative  nonlinear  Schrodinger  equation" 
[Mjolhtis,  1976;  Spangler,  1986;  Hada  et  al,  1989].  The 
possible  applications  of  Alfven  solitons  to  solar  and 
astrophysical  plasmas  have  been  discussed  by  Ovenden  et 
al  [1983] .and  Spangler  [1985]. 

While,  significant  research  has  been  directed  to  the 
Alfven -waves,  the  magnetosonic  waves  have  received  much 
less  attention.  The  magnetosonic  waves  are  hybrid  and 
elliptically  polarized  modes  that  can  propagate  obliquely 
across  the  magnetic  field.  They  become  decoupled  from 
the  Alfven  waves  when  the  direction  of  their  propagation 
is  exactly  perpendicular  to  the  ambient  magnetic  field. 
The  purpose  of  the  present  work  is  to  investigate  the 
filamentation  instability  of  the  magnetosonic  waves  that 
propagate  perpendicularly  across  the  magnetic  field.  This 
instability-  excited  obliquely  propagating  magnetosonic 
sideband  modes.  The  concomitantly  excited  nonoscillatory 
(purely  growing)  modes  have  their  wave  normal  parallel  to 
the  magnetic  field.  By  contrast,  the  purely  growing 
modes  excited  by  the  filamentation  instability  of  Alfven 
waves  are  magnetic  field-aligned  [Kuo  et  al.,  1988,  1989]. 
Hence,  the  characteristics  of  the  filamen ration  instability 
under  consideration  are  distinctively  different  from  those 
of  the  instability  discussed  in  our  previous  work  [Kuo  et- 
aL,  1988,  1989], 

This  paper  is  organized  as  follows.  The  coupled 
mode  equations  are  derived  in  section  2.  They  show  how 
the  magnetosonic  sidebands  and  the  purely  growing  modes 
are  parametrically  coupled  through  the  magnetosonic  pump 
wave.  In  section  3  the  dispersion.relation  is  obtained  and 
analyzed  to  derive  the  threshold  conditions  and  the 
growth:  .rates  of  the  filamentation  instability.  The 
relevance  of  the  present  work  to  some  observations  is 
discussed.  Finally  presented  in  section  4  are  summary  and 
conclusions'. 


Filamentation  Instability  of  Magnetosonic  Waves 


Consider  a  magnetosonic  wave  propagating 
perpendicularly  across,  the  interplanetary  magnetic  field 

I?0«zB0.  Its  wave  magnetic  field  is  represented  by 


Bp.?  zBpCxplKkoX-Wot)]  +  c.c.  (1) 

where  the  unperturbed  wave-  field  intensity  is 
assumedftobeconstant;  and.a::rral -quahtityfor.rimplimty. 
The  fwave,  number .  kj; 

are  related'  by  the  dispersion  felatipn 
and  the  ion  acoustic  speed,  respecti  vejy.  v  The^  yelpcity 
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responses  of  electrons  and  ions  to  the  magnetosonic  pump 
wave  can  .be  written  as 

Vp.  =  [x+i(  k^V^iA)  yK  1Vko)  ( B^B0>  exp{i(  hpX-ept)  3 
+  c.c. 

'Vpi  =  xfWo/yfBp/Sofcxrtit^tJi  +  c.c.  (2) 


In  the  following  we  will  investigate  a  filamentation 
instability  which  c as  break  up  a  large-amplitude 
magnetosonic  wave  into  filaments.  The  magnetosonic  wave 
is  assumed  to  be  initially  uniform  as  described  by  ( 1).  for 
simplicity.  The  filamentation  instability  results  from  small 
fluctuations  in  plasma  density  which  lead  to  modulation  Of 
plasma  dielectric  constant  and  spatial  wave  distribution, 
which,  in  turn,  enhance  the  plasma  density  fluctuations. 
Such  a  positive  feedback  process  involves  a  four-wave 
coupling  process:  through  the  magnetosonic  pump  wave 

(w0,  kc),  two  side-band  perturbations  (iu±)  are 
coupled  with  a  purely  growing  density  perturbation 
(w,kz).  The  wave  frequency  and  wave  vector 

matching  conditions,  are  u+-w=w0ew.+«*  and  1c  ±=Tcp±k,  are 
satisfied  where  the  asterisk  denotes  the  complex 
conjugate.  Illustrated  in  Figure  1  is  the  wavevector 
matching  condition,  showing  the  relative  orientation  of 
these  wave  vectors.  It  will  be  sown  that  these 
perturbations  can  be  excited  simultaneously.  In  other 
words,  they  grow  exponentially  with  the  same  growth 
rate,  q=-iw  at  the  expense  of  magnetosonic  pump  wave 
energy.  This  instability  requires  a  threshold  to  occur, 
which  is  going  to  be  determined  later.  In  the  analysis  of 
filamentation  instability,  pump  field  intensity  is  assumed 
to  be  a  constant  maintained  by  a  source.  Pump  depletion 
effect  which  may  contribute  to  the  saturation  of 
filamentation  instability  is,  therefore.ignored  within  the 
present  framework  of  instability  analysis. 

Our  proposed  theory  is  developed  on  the  basis  of  the 
two-fluid  plasma  model.  The  fluid  model  is,  in  general, 
valid  provided  that  the  condition  k2<^/rf«l  is 
satisfied  where  k,  c,,  and  are  the  wave  number 
of  the  purely  growing  mode,  the  ion  acoustic  speed,  and 
the  ion  gyrofrequency,  respectively.  Therefore,  the 
analysis  of  filamentation  instability  is  limited  to  the 
region  of  k/kp<1.2.  To  analyze  the  instability. 


Fig.  1.  The  relative  orientation  and  wave  vector  matching 
relation  for  the  magnetosonic  pump  (too  jco),  two  sidebands 
(©i.y  and  the  purely  growing  mode  (co  jc)  as  required  in  the 
proposed  four-wavecduplihg  process.  . 


the  complete-  set  of  fluid  equations  will  be  linearized  with 
respect,  to  the  first-order  perturbations  in  physical 
quantities  such  as  density  (n) ,  velocity  (v),  wave  electric 

field  (E)  and  wave  magnetic  field  (if),  These 
first-order  vphysical- quantities,  associated  with  the  sideband 

modes  are  denoted  by  <5n±,  5V,j+,5Ii±,  and  S~§+  those 
associated  with  the  purely  growing  mode  n,, 

"v,,*,  acc*  18,.  Thus,  the  total  density  is  given 
by  n=np+ns+(  Sn++5n_+c.c.) ;  similarly,  "v«j="v«4 
+c^W5'  ^ej++£v,j_TC.c:)‘  "2 *'!;,+(  EpfSE++£E_+c. c.) , 

and  B=B0+{  B^+5^+*5B-.+c.c.)  for  the  concerned  four-wave 
interaction  process. 

The  purely  growing  mode  is  associated  with 
electrostatic  disturbances,  varying  spatially  along  the 

magnetic  field,  iie.,  k=zk.  The  fluid  equations 
used  to-  derive  the  coupled  mode  equation  for  the  purely 
growing  mode  include  the  linearized  continuity  equations 
for  electrons-  and  ions: 

IT  +  no'i|  v»«»  =  0  *  Jf  n.+no  ^  V»U  (3) 

and  the  combined  linearized  equations  of  motion  for 
electrons  and  ions: 

m.  I  v.«  +  |  v»ir.  *  (F„ 

+  Fu)  -  -  m,  (cj/np)  |  n. 

In  deriving  the  above  equations,  quasi -neutrality  has  been 
assumed.  F^  and  Fif  represent  the  z  components  of  the 
nonlinear  Lorentz  forces  experienced  by  electrons  and 
ions,  respectively,  namely,  F,^*<-mt^( "v •  V~vej) 

+  q#  ("ve^/ c)xB>,  where  angle  brackets  denote  the 
time  average  over  the  pump  wave  period.  The  explicit 
expressions  of  are  given  by 


*Vm.l  =  |  (Vp,.  .. 

+i<*>  (-4?  )  £  fKlxViy  '  vp 'IxKh 


Sv%iy)*c.c 


where  £Vei=fv^++f?ej_  iS  tlle  velocity  responses  of 
electrons  and  ions  to  the  combined  sideband  mode; 
<^(«j)-+l(-l).  The  expression,  (4)  reduces  to  the 
form- of  .ponderoraotive  force  in  an  unraagnelized-pfasma. 
.Without  losing  the  generality,  the  density -fluctuations  of 
the  purely'growing'mode  car.  be  expressed  as 

n,  =  n,  exp(Tt)cos  kz 

where  n,  and  7  ate  the  .  realF  amplitude-  and  -the 
growth  rate  of  the  density  fltKtuatjbns^respwtiwJy; 

Combining  the  equations  in  (3)  yields  the  coupled 
mode  equation  for  the  purely  growing'mode  as  ;  1 

(M’.KS./iV  =  (k/miKf^+Ffc)  (5) 
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where  Fe;-I=Fe^exp(7t)'sinkz  is  assumed.  One  can 
see  from  (.4)  and  (5)  that  the  density  fluctuations 
associated  with  the  purely  growing  mode  stem  from  the 
differential  interaction  of  the  magnetosonic  pump  wave 
and  sideband  modes.  Hence,  expressions  for  the  physical 
quantities  of  the  combined  sideband  mode  have  to  be 
determined  in  a  '  self-consistent  manner.  The  purely 
growing  mode  is  realized  to  lie  a  nonlinearly  driven  mode. 
However,  in  the  .‘absence  of  the  nonlinear  Lorentz  force, 
(5)  reduces  to  the  linear  dispersion  relation  of  the  ion 
acoustic  wave. 

The  magnetosonic  sideband  modes  are  excited  through 
the  beating  current  driven  by  the  magnetosonic  pump 
wave  field  in  the  density  fluctuations  of  the  purely 
growing  mode.  It  should  be  noted  that,  since  the 
magnetosonic  pump  wave  propagates  perpendicularly  across 
the  magnetic  field  and  the  purely  growing  mode  oscillates 
spatially  along  the  magnetic  field,  two  magnetosonic 
sidebands  are  excited  symmetrically  around  the  pump  wave 
and  they  propagate  obliquely  across  the  magnetic  field 
(see  Figure  1).  Thus,  these  two  sideband  modes  propagate 
together  in  the  propagation  direction  of  the  pump  wave 
and  form  a  standing  wave  pattern  along  the  magnetic 
field.  In  the  following  derivation  of  the  coupled  mode 
equations  for  the  two  sideband  modes,  we  can,  therefore, 
combine  them  into  a  single  coupled  mode  equation  for  the 
combined  sideband  mode. 

This  combined  sideband  mode  is  characterized  by  the 
following  physical  quantities  with  self-consistent 
expressions: 

fin  ■  Sn  exp(Tt)cos  kz  exp£i(  l^,x-t*r0t)  ]  +  c.c. 

=  K&v.jx  +  iy«vejy)cos  kz  +  izSveis  sin  kz]exp(Tt) 

•  exp[i(k^x-t^t)]  +  c.c. 

SB  =  [(k/koK-ix+ayJsin  kz  +  2  cos  kz]SB  expfqt) 

•  exp[i(k„x-w0t)J  +  c.c.  (6) 

where  a,  an  unknown  parameter  defining  the  polarizations 
of  the  sideband  modes,  will  be  determined  later  on.  The 
fluid  equations,  used  to  derive  the  coupled  mode  equation 
for  the  combined  sideband  mode,  include  the  linearized 
continuity  equations  for  electrons  and  ions: 

|  (fci/no)  +  V  .  [£ve  +  (n,/n0)’vpt]  = 

0  *  Jf  (®“/"o)  +  v  •  t^i  +  (n»/n0)"vpj]  (7) 

the  linearized  Maxwell  equation  for  induction  law: 

|  («B/B0)  -  V  x  (f?t  x  2)  (8) 

the  combined  equations  of  motion  for  electrons  and 
ions: 

fv \  =  -  c2,  V  (Sn/iig)  +  (\(6vrlTvt)  x  2  (9) 

and  the  linearized  Maxwell  equation  for  Amperes’  law: 

V  x  (5B/Bo)  =  (fyv*A)[(*?r*V#) 

.  +,  (V^Ktpj  - V,,,)]  (10) 


Substituting  (6)  into  equations  (7) -(10),  we  obtain 

=  (<V*o)(«B/B0) 

Mey  =  Q<VkoK«B/B0) 

=  -  («o/k0)[(k2oC2./n2i)(®/no)- 

( l-aWo/^HiB/Bo)]^  1-uytfj) 

^iy  =  (VtyUkoC2, /<■><))(  no)  +  Mfyw0)-1] 

(w^koHJg/Bo)}^  l-wytf;) 

S?u  =  (kcVfc’oKffi/no) 

+  °<kvyfy(«g/B0) 

®/no  =  {[( l-ooytyA  i-^o/n2;)  K«8/b0) 

+  (V^)(Bp/Bo)>/[l-k*cy^o  +  kW^r^O )]  01) 

Using  the  results  in  (11)  together  with  the  x  and  y 
components  of  ( 10) ,  we  can  derive  the  following  coupled 
mode  equation  for  the  combined  sideband  mode: 

[“V^+k^MvVc2,)  +  k2(k2+k20)vJAc2,/“2o 
+  *  W[kVVwJo(  l'l^,/^0)])(ffl/B0) 

=  -  ^oI^aO-^./w1,,)  -  c2,](n./n0)(Bp/B0)  (12) 

as  well  as  the  expression  for  the  parameter  « 

“  =  (Wo/^)[l  ^(k^kVcy^ol/ 

[  1-c2,/vja-(  kVA/^0)  ( l-kV./^o)  1  ( 13) 

In  the  derivation  of  (12)  and  (13),  « 

1  has  been  assumed.  Equation  ( 12)  shows  that  the 
magnetosonic  sideband  fields  are  produced  by  the 
magnetosonic  pump  field  in  the  density  fluctuations  of  the 
purely  growing  mode.  Equations  (5)  and  (12)  thus  form 
the  complete  set  of  coupled  mode  equations,  which  can  be 
based  on  to  analyze  the  proposed  filamentation  instability 
of  magnetosonic  waves  in  the  solar  wind  environment. 

Aaalysis  of  the  Proposed  Instability 

The  dispersion  relation  of  the  concerned  instability  can 
be  derived  from  the  two  coupled  mode  equations,  (5)  and 
(12).  Its  derivation  is  briefly  described  as  follows.  With 
the  aid  of  (2),  (4),  (6),  and  (11),  the  coupling  term  on  the 
right-hand  side  of  (5)  can  be  expressed  as 

k(Fe+P.)/mi  =  k’lvVc^I-c*^)/ 

(l-k’cV^oKBp/BoKSB/Bo  +  «g‘/Bo)  . 

-  2{l?cl/(l-kJci./«*0)](»./n0)(Bp/By*  (14) 

Substituting^  14)  into  (5)  and  combining  it  wth  ,(  12)  yield 
the  dispersion  relation  :j--  ■- 
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K/K0 

Fig.  2.  a/5  versus  k/k0  for  several  selected  values  of  P  in  the  range  of  0.25-2.5;  it  shows  how  the 

polarization  of  the  excited  magnetosonic  sidebands  is  determined  by  the  characteristics  of  the  magnetosonic 
pump  waves  and  the  environmental  plasmas. 
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f+kV,  =  2k20vJA{(i-«l+/l)/(/M-^2/ir0)] 
ll-PtU-P)(l-a5)/(M-tf/y?0)] 
(^l-^/^ol/K/SUlMl+^/l^o)  *  aSUb  Mltg/lt-fi] 

(15) 

where  P=<^t/^\  is  the  ratio  of  the  plasma 
kinetic  pressure  to  the  background  magnetic  pressure,  and 
£=“0/ In¬ 
setting  t=0  in  (15)  determines  the  threshold  field  of 
the  filamentatioB  instability 

(Bpth/B0)2  =  (k2/2k10)^(/M)J-/E(l+k2/kJ0) 

+  o£(/hl)(kyk2-0]/{I-/W)/ 

(/W-/3t2/ks0l[l-«^l)(  ] 

(pH-pf/fy?  -  ffl+k2/^) 

+  Q$^l)(^o/k?-/S)]/(^l-^/^0)}  (16) 

In  terras  of  (16),  the  growth  rate  of  the  instability 
derived  from  (15)  is  found  to  be 

1  =  [(Bp/Bpth)J-l]1/Jkc,  (17) 

The  threshold  condition  defined  by  ( 16)  is 
conveniently  expressed-  for  the  normalized  threshold  field 
intensity  (Bpth/B0)  as  a  function  of  the 

normalized  wave  number  (k/kj,).  The  parameters- 
P(  "C?,  /  v?A) ,  £(«<%/ 1^):,  and:  a  (  defined 
by  (I3)j  characterize  the  solar  wind  plasinas  and  the 


propagating  magnetosonic  waves.  In  fact,  the  parameter 
a  can  also  be  expressed  in  terms  of  (k/k^>),  6,  and  p 
as  follows: 

o  =  5{I-2{l-KVko)V/(l+^/{l-^[Oc/kb)2/ 

(l+«][l-(k/k0)2/?/(l+«])  (18) 


Displayed  in  Figure  2  is  ( a/S)  as  a  function  of  (k/kg) 


for  several  values  of  p  chosen  from  the  range  of  0.25- 
Z5.  The  result  shows  that  a  is,  in  general,  larger  than 
one  for  5<0.3  which  is  the  region  of  practical 
interest.  Especially,  a»l  for  P<  1  and  £«1.  In  this 
parameter  range,  the  decay  sidebands  become  linear 
polarization  on  the  y-z  plane  and  are  dominated  in  the 
nature  of  shear  mode. 

The  threshold  fields  (BpUj/B0)  as  a  function 
of  ( k/kg)  are  shown  in  Figures  3-6.,  For  illustrative 
purposes,  the  values-  of  P  are  again  selected  from  the 
range  of  0.25-2.5  and  5=0.01,  0.03,  0.06,  or  0.3  is  chosen. 
The  corresponding  growth  rates,  for  (Bp/B0)J=0.05, 
as  a  function  of  k  are  given  in  Figures  7-10.  Several 


outstanding  features  of  the  proposed  instability  can  be 
seen  from  these  figures.  For  instance,  the'  threshold  field 
generally  increases  with  k,  the  wave  number  of  the 
excited  purely  growing  mode,  and  it  finally  reaches 
infinity  at  certain  k<kQ.  While  the  threshold  field 
increases  with  6  in'  the  region  0<jS<l,  it  decreases  . with 
p  for  fi>  1.  The  growth  rate  has  a  peak  value  (denoted 
by  7„„)  at  certain  k  (denoted  by  k^).  We 
note  that  both  7mtx  and  k,,,  decrease  with  P 
for  0<6<1,'  but  they  increase  with  P  for„.£»l.  By 
contrast,  while  T^..  decreases-,  with  £,  k^ 
increases  with  £.  The  thresholds  and  growth  rates  of 
the  instability  vary; with  the  magnetosonic  wave  frequency. 
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Fig.  3.  The  normalized  threshold  field  intensity  (Bp/Bo)  versus  (k/ko)  with  several  selected  values  of  P  in  the 
range  of  0.25-2.5  for  8=0.01. 


It  is  found  that  the  threshold  fields  (growth  rates) 
increase  (decrease)  with  the  magnetosonic  wave 
frequency. 

Let  us  base  upon  the  observed  characteristics  of 
magnetosonic  waves  (Hoppe  and  Russell,  1983;  and 


references  therein]  to  examine  the  excitation  of  the 
proposed  filamentation  instability.  The  magnetosonic 
waves  observed  by  the  ISEE  1  and  2  typically  have  low 
frequencies  of  order  of  0.01  Hz,  which  is  about  0.1  times 
the  local  proton  gyrofrequency  (i.e.,  S=w0/f^=0. 1). 
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The  local  interplanetary  magnetic  field  (B0),  therefore, 
is  of  order  of  100  nT.  The  magnetosonic  waves  typically 
have  amplitudes  (Bp)  of  about  10  nT,  namely, 
Bp/B0  ~  0.1  with  wavelengths  ~1  RE.  The  parameters 
adopted  in  the  analysis  of  our  proposed  instability  are 
generally  in  agreement  with  these  observations.  Under 


the  condition  of  Bp/B0  ~  0.1,  the  thresholds  of 
the  proposed  instability  can  be  exceeded  for  k/k„ 
(=Aq/A)  less  than  0.1  or  so,  where  and  A  are  the 
wavelength  of  the  magnetosonic  wave  and  the  scale  length 
of  the  excited  purely  growing  mode,  respectively.  An 
important  prediction  of  the  proposed  process  is  that 
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Fig.  7.  The  growth  rate  (Y),  assuming  (Bp/Bo)2=0.05,  as  a  function  of  k  for  5=0.02. 


plasma  density  fluctuations  with  scale  lengths  longer  than 
at  least,  a  few  RE  ( Earth  radius)  can  be  generated  along 
the  magnetic  field  by  magnetosonic  waves.  The 
polarization  of  the  obliquely  propagating  magnetosonic 
waves  is  determined  by  the  value  of  a  which  in  turn,  is 
a  function  of  the  characteristic  parameters  of  the 
environment:  8,  p,  and  ko  (see  equation  (18)). 


However,  we  note  that  there  seems  to  exist  a  discrepancy 
between  the  observed  and  the  predicted  propagation 
directions  of  magnetosonic  waves.  The  waves  upstream  of 
Earth’s  bow  shock  with  frequencies  of  0.01  Hz  were 
observed  to  propagate  nearly  parallel  to  the  magnetic 
Held,  while  the  predicted  magnetosonic  waves  propagate 
obliquely  across  the  field.  We  suspect  that  this 


Fig.  8.  Yvetsos  (k/ko)  for  5=0.03. 


‘  734 

Phys.  Fluids  B1  (4),  April  1369 

-  0699-8221 / 89/040734-07301 .90 

©1969  American  Institute  of  Physics  '  734 

— f Jr_-  __ 

- 

Kuo  and  Lee  Filamentation  of'Magnetosonic  Wave 


for  S=0.06 


0-025 


Fig.  9.  y  versus  (k/ko)  for  5=0.06. 


discrqaiicy  may  arise  from  the  intrinsic  difficulty  in-thc 
in  situ  measurements  of  upstream  wtves  by  st-ellites. 

U*  speed  of  solar  wind  is  far  greater  than  those  of 
satellites,  all  measures  waves  appear  to  be  nearly  parallel 
propagating.  The  predicted  oblique  propagation  of 
magnetosonic  waves  excitedly  niamentatioc  instability 
awaits  the  corroboration  of  future  observations. 


■4.  Summary  and  Conclusions 

The  excitation  of  a-  parametric  instability  -by 
magnetosqnic  waves  in  the  solar  wind  environment  is 
proposed.  The  magnetosonic  pump  wave  originally 
propagates  at  the  right  angle  with  the  magnetic  field. 
This  instability  can  excite  a  purely  growing  mode  along 
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tbs  magnetic  field-  together  with-  two  symmetric 
magnetosonic  sidebands  propagating  obliquely  across  the 
magnetic  field.  These  two  sideband  modes  thus  propagate 
together  perpendicularly  across  the  magnetic  field  aadr 
meanwhile,  form  a  standing  wave  pattern  along  the 
magnetic  field.  The  physical  process  responsible  for  the 
instability  can  be  described  a*  the  scattering  of  the 
magnetosonic  pump  wave  into  -obliquely  propagating 
sidebands  by  the  simultaneously  excited  plasma  density 
fluctuations,  which  are  associated  with  a  purely  growing 
mode.  This  instability  leads  to  the  fibmentation  of  the 
magnetosonic  pump  wave.  The  thresholds  of  this 
fi lamentation  instability  can,  indeed,  be  exceeded  in  the 
soiar  wind  environment.  The  .density  fluctuations 
produced  by  the  fi lamentation  instability  along  tbb 
magnetic  field  are  predicted  to  have  wavelengths : greater 
than,  at  least,  a  few  Earth  radii.  As  shows  in  (18),  tbe 
polarization  of  the  obliquely  propagating  magnetosonic 
waves  excited  by  the  filamentation  instability  is 
determined  by  the  environmental  characteristic  parameters: 

P(  =  e\/^A)»  and  kj.  The 

predicted  oblique  propagation  of  magnetosonic  waves, 
however,  has  hot  been  confirmed  by  the  available 
observations  yet 
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Abstract — Possible  generation  of  whistler  waves  by  Tromso  HF  heater  is  investigated.  It  is  shown  that  the 
HF  heater  wave  can  parametrically  decay  into  a  whistler  wave  and  a  Langmuir  wave.  Since  whistler  waves 
may  have  a  broad  range  of  frequency,  the  simultaneously  excited  Langmuir  waves  can  have  a  much  broader 
frequency  bandwidth  than  those  excited  by  the  parametric  decay  instability. 


1.  INTRODUCTION 

There  is  considerable  interest  in  investigating  gener¬ 
ation  and  amplification  of  VLF/ELF/ULF  waves  in 
the  ionosphere  by  high  power  HF  waves  (see.  e.g., 
Stubbe  el  al. ,  1982:  Lunnen  el  al..  1985).  Several 
processes  have  been  suggested  as  potential  mech¬ 
anisms  causing  them  such  as  the  modulation  of  elec¬ 
trojet  current  (Stubbe  and  Kopka.  1981 ;  Fejer  and 
Krenzien.  1982).  the  thermal  filamentation  instability 
of  HF  heater  (Kuo  and  Lee.  1983)  and  parametric 
excitation  of  Alfvcn  (ELF)  waves  (Papadopoulos  et 
al.,  1982).  In  this  paper,  we  discuss  a  possible  mech¬ 
anism  for  VLF/LF  wave  generation  in  the  ionosphere 
by  HF  heater. 

As  described  in  Sections  2  and  3.  the  proposed 
mechanism  is  the  parametric  excitation  of  whistler 
waves  by  HF  heater.  The  concomitantly  excited  elec¬ 
trostatic  waves  arc  Langmuir  waves.  Thus,  the  con¬ 
sidered  process  is  the  parametric  decay  of  HF  heater 
wave  into  a  Langmuir  wave  and  a  whistler  wave. 
This  process  is  different  from  the  so-called  parametric 
decay  instability  of  HF  heater  waves  on  two  aspects. 
One  is  that  the  parametric  decay  instability  generates 
a  Langmuir  wave  and  an  ion  acoustic  wave ;  both  of 
the  excited  waves  are  electrostatic  mode  waves.  The 
other  one  is  that  the  ion  acoustic  wave  has  a  frequency 
of.  typically,  a  few  kHz  in  the  ionosphere,  while  the 
whistler  wave  may  have  a  broad  range  of  frequency 
covering  VLF,  LF,  and  even  MF  bands.  Therefore. 
Langmuir  waves  excited  by  the  proposed  mechanism 
can  have  a  broader  frequency  bandwidth  than  those 
pt  'duced  by  the  parametric  decay  instability. 

Observations  of  the  HF  heater-produced  VLF / 


ELF.LiLF  waves  were  constantly  made  at  Tromso. 
Norway  with  the  E1SCAT  heating  facilities.  Stimu¬ 
lated  electromagnetic  emissions  with  a  broad  frequency 
band  of  several  tens  of  kHz  were  seen  in  the  Tromso 
heating  experiments  (Stubbe  et  al..  1984  and  refer¬ 
ences  therein).  Hence,  our  proposed  mechanism  is 
specifically  analyzed  for  high-latitude  heating  experi¬ 
ments  with  the  intention  to  compare  the  theoretical 
results  with  observations  made  at  Tromso.  Norway. 
A  brief  discussion  is  finally  given  in  Section  4. 


2.  THEORY 

As_rncntioned  earlier,  the  EISCAT  heating  facilities 
at  Trdmso.  Norway  were  observed  to  produce  VLF/ 
ELF/ULF  waves  in  the  ionosphere.  The  Tromso 
HF  transmitter  is  steerable  to  launch  heater  waves 
propagating  along  the  earth's  magnetic  field.  In  view 
of  these  facts,  we  will  present  our  theory  as  follows  in  a 
simple  geometry  of  propagation,  which  is  appropriate 
for  describing  the  parametric  excitation  of  whistler 
waves  by  HF  heater  at  Tromso.  Norway. 

Consider  a  left-handed  (L-H)  circularly  polarized 
heater  wave  propagating  upward  along  a  downward 
geomagnetic  field  in  an  overdense  ionosphere.  The 
ionospheric  density  («)  is  assumed  to  have  a  linear 
profile,  increasing  with  altitudes  (r)  as  n  =  «0(1  +:/L) 
where  L  is  the  scale  size  of  the  background  iono¬ 
spheric  inhomogeneity.  The  heater  wave  is  reflected 
at  a  critical  height  before  reaching  the  ionospheric  F 
peak  and  becomes  a  standing  wave,  whose  energy 
follows  an  Airy  function  distribution  illustrated  in 
Fig.  1. 
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Fig.  1  Standing  wave  pattern  ol' heater  wave  near  the  reflec¬ 
tion  height. 


The  first  peak  of  the  Airy  function  is  located  at 
z  =  0  and  its  width  is  denoted  by  2A.  The  heater  wave 
field  in  this  region  is  then  modeled  by  Er  =  (.v— n‘)er 
4-r.e..  where  r.p  =  [e„,  ( 1  —  /riA)]  exp  ( —  uont)  whose 
phasor.  «...  (l-/r  A),  is  proportional  to  a  complex 
Lorcntzian  function.  The  velocity  response  of  elec¬ 
trons  to  the  heater  field  is  found  to  be 

V,  =  —  /(.v  -  i)  )eF.pim[(u) „  ~fi,)  +  h\]  +  c.c.  ( 1 ) 

where  v,  is  the  electron  collision  frequency. 

We  will  show  that  the  parametric  decay  of  the 
heater  wave  (<o„.k„)  into  a  downgoing  (or  upgotng) 
whistler  wave  (w.k)  and  an  upgoing  (or  downgoing) 
Langmuir  wave  (<y,.k,)  can  occur  in  the  aforemen¬ 
tioned  peak  region  of  heater  wave  field.  The  following 
wave  frequency  and  wave  vector  relations.  <u„  =  <!>+<•>, 
andk(O)  =  —  k,(0)  =  ~:k.  arc  satisfied.  First  derived 
is  the  coupled  mode  equation  for  the  whistler  wave, 
which  is  excited  by  the  beating  current  driven  by  ths; 
pump  field  in  the  density  perturbation  associated  with 
the  excited  Langmuir  wave. 

The  density  of  beating  current  as  the  driving  source 
of  the  whistler  wave  is  given  by  =  —e$n?Vr 
where  <5 «,(  =  <5«,cxp(— hi,t)-t-c.c.)  is  the  density  per¬ 
turbation  associated  with  the  Langmuir  wave:  Vr 
shown  in  (I),  is  the  velocity  response  of  electrons 
to  the  heater  wave  field.  The  density  perturbation 
can  be  obtained  from  the  Poisson's  equation,  on,  = 
~(l/4iH?)(r£//cr),  where  s,  is  the  amplitude  of 
Langmuir  wave  field  <5E,=  re,  =  re,  exp  (-/w,/).  In 
addition  to  the  beating  current,  a  linear  current  is 
caused  by  the  whistler  wave  with  the  background 
plasma,  <5J„  =  —end\u  where  SV„  is  the  velocity 
response  of  electrons  to  the  whistler  field  and  n  is  the 
background  electron  density ;  <)V(,  has  a  similar  form 
as  SVP ,  namely, 

'5V„  =  -  / (.v  -iy )e&Jm[{w  ~  Cle)  +  ivr] + c.c.  (2) 


where  e,,  is  the  amplitude  of  the  whistler  wave  field. 
(5E,„  =  (.v— iy)e.„+c.c.  =  (x—iy)e„  exp  ( — ia>,) + c.c. 
Substituting  the  respective  expressions  for  and 
<5J/;  into  the  whistler  wave  equation  yields  the  coupF'd 
mode  equation 

i  I-,  +  /[o»;(c)/(to — +  iv, )] -  -r—^fe,, 

(<  /  ct  iz  j 


=  i[elm(w„ -Q,+  />,)]  ~  (cp ^ e,*j  (3) 


where  ior(z)  =  Ane'n(:)jm. 

We  next  derive  the  coupled  mode  equation  for  the 
Langmuir  wave.  The  driving  force  of  the  plasma  oscil¬ 
lation  is  the  Lorentz  force  experienced  by  the  oscil¬ 
lating  electrons  in  the  wave  magnetic  fields.  Respond¬ 
ing  to  the  electric  fields  of  the  heater  wave  and  whistler 
wave,  electrons  oscillate,  respectively,  with  velocities 
V,  and  (W.,  transverse  to  the  background  magnetic 
field.  These  cross-field  motions  interact  with  the  wave 
magnetic  fields.  <)B, ,  and  dB,,.  and  result  in  a  longitudinal 
Lorentz  force.  -(r/f)(iiV„xBf+ V^xPB,*)  which  is 
responsible  for  the  longitudinal  plasma  oscillation. 
Including  this  force  in  the  electron  momentum  equa¬ 
tion  and  combining  this  equation  with  the  electron 
continuity  equation,  we  obtain  the  following  coupled 
mode  equation  for  the  Langmuir  wave 

r- 

=  -  2 (cw;//»)  j^e,  r*Joj*(wn  -Qc  +  tv,) 

+  e*  1  riu)„(w*-Qf-  ivr)  (4) 
cz 

where  V,  is  the  thermal  velocity  of  electrons. 

It  is  clear  from  the  two  coupled  mode  equations. 
(3)  and  (4),  that  the  whistler  wave  and  Langmuir  wave 
are  coupled  through  the  HF  heater  wave,  which  acts 
as  the  pump  of  the  concerned  parametric  instability. 
To  proceed  the  instability  analysis  further,  it  is  desir¬ 
able  to  reduce  these  two  equations  into  first-order 
differential  equations.  This  can  be  done  by  intro¬ 
ducing  the  transformations  (Rosenbluth,  1972: 
White  era/.,  1973), 


f.,~  a,  exp 


j/  A,;-W|/+J^ 


AA-,(.-)dc 


y  *  * 


! 
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=  a:  exp 


i  A,r--w-/4 
_  '  Jo 


AA:(r')dr 


(6) 


where  A,  =  A,(0)  =  A  =  — A:(0),  tof  =  c/,;  =  to^'(O) 
+k]V;  and  co,  =  co  =  A ic2Qelioj;(0),  AA,  =  [a)_;(0) 
— a>;(;)]/2A, KJ  and  AA,  =  (oj/2A:Qrr;)[w'(r)  —  w; 
(0)].  Note  that  ©i+w,- «i0  =  <5  for  a  mismatch  fre¬ 
quency  and  AA  =  AA ,  4  AA,  =  A'r,  where  A'  =  —  (oi„(0)/ 
2kL)((o/Qec2+  1/A',). 

Substituting  (5)  and  (6)  into  (3)  and  (4)  and 
assuming  that  |ra,/cfj  «  wja,!,  |«/,/ct|  «  o>|a,j. 
\cajc:\  «  |A|||a,|.  and  |ca;/or|  «  |A2||«,1,  we  derive 
the  following  first-order  differential  equations  for  o, 
and  a. 


{eim)[(o-(:)klw ,  w(«o  -  Q,)j 


•«,«?exp[-/A'eV2]  (7) 
and 


(Jr  -  ]«;=  (i''/M)[Q,(w;-<j)A/ft>;(r)tw(1-fl,)] 

exp  [  -  ik'z2j2)  (8) 

where  K,  =  AF;  w,  and  l':  =  2Ac:Qr«;(0)  are  the 
group  velocities  of  Langmuir  wave  and  whistler  wave, 
respectively:  v,  =  v,/2  and  v,  =  v/o/Cl,  The  simpli¬ 
fied  coupled  mode  equations.  (7)  and  (8),  will  be 
ana! yzed  in  the  next  section  for  the  parametric  exci¬ 
tation  of  Langmuir  wave  and  whistler  wave  by  a  non- 
uniform  HF  hea'er  wave  in  a  nonunifoitr.  ionospheric 
plasma. 


3.  ANALYSIS 

After  being  Laplace-transformed  in  time,  equations 
(7)  and  (8)  can  be  written  as 

^  A  +  v  i  4  K,  ~J 4,  =  (e/7n)[w;(r)A/to:a»j(w0-Q,)] 

•  arA  ?  exp  (  -  /A'"/2) 4 a  10(*)  (9) 

\^+Vi+iS..vi0At 

-  (elm)[nr(a>2  - <5)A/(u;(;)(cu0  -  Q,)] 

•  apAf  exp  ( —  ik':'/2)  -F  a20(r)  ( 10) 

wherep  is  the  Laplace  transform  variable.  A ,  =  ^[at] 
and  flj  =  Se[a2],  where  if  [  ]-  J*  [  ]exp(-pt)d/  rep¬ 
resents  the  Laplace  transform;  al0(z)  and  a20(z)  are 
the  initial  values  of  at  and  a2,  respectively.  In  the 


following  eigenmode  analysis,  a  ,„(;:)  and  a20{z)  can 
be  set  to  be  zeros. 

Combining  (9)  and  (10)  leads  to  a  second-order 
differential  equation  for  A , : 


+/)(:)—+/:(:) 


A,  =0 


(H) 


where 


/;(-)  =  (7>+v,)/F,  ~(p4v2-/<5)/F, 

4  [ik'z  - 1  l(L+:)~  i/(A  -  /.')]  ( 1 2) 


and 


j:(:)  =  (eim)z[k2nr((o2~S)/u)f(02 

x  (ft>o— K,(!  4c:.A3) 

+  [(/>+ v ,)/  V,\  [/A'r  -  1  l(L  4  r)  -  //(A  -  ij)J 

—  [(/?+vi)/F|)[(p4v:  — /t))/K  j],  (13) 

Introducing  that  f,(r)  -=  4'(r)exp[- rj;",y , (-') dr'], 
equation  (II)  can  be  further  transformed  into 


+m 


'F(r)  =  0 


(14) 


where /  (_•)  =./;(r)-(l  2)[(L2)/;(r)  +  (d,d.-)/,(r)]. 

Both  the  effects  of  the  background  density  inhomo¬ 
geneity  (i.e„  f^x)  and  finite  inrtability  zone  (i.e.. 
A  ^  oo )  are  taken  into  account  and  included  in  f(z). 
However,  since  L  »  A.  the  effect  of  background  den¬ 
sity  inhomogeneity  is  negligible  in  comparison  with 
that  of  the  finite  instability  zone.  Therefore.  /'(-)  can  be 
approximately  expressed  as 

m  -x-'/(14--7A’)-(1  /4)(( p 4 v , )/ K, 

4(p4v:  — /V5)/K.  — /A'r4//(A-/r)]:  — ( 1/2) 

x[/A'4l/(A-/r)!]  (15) 

where 

x:  =  (e/m)i[A2fl,.(e>j— d)/o>|Cu,(co0— Q,):]eo/F,  F2. 

The  eigenvalue  of  (14)  needs  to  be  determined  for 
finding  the  threshold  fields  of  the  heater  wave  and  the 
growth  rates  of  the  instability.  The  requirement  of 
threshold  fields  is,  in  general,  imposed  by  two  additive 
effects.  They  are  the  coiiisional  loss  of  excited  modes 
and  the  convective  loss  across  the  boundaries  of  finite 
instability  zone.  These  two  effects  are  evaluated  sep¬ 
arately  as  follows. 


(a)  Coiiisional  loss 

We  first  consider  the  damping  of  excited  modes  due 
to  collisions  by  letting  A  -+  oo.  It  is  equivalent  to  the 
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study  of  instability  in  a  uniform,  unbounded  iono¬ 
spheric  plasma.  Then.  A ,  and  /,(z)  become  spatially 
independent  functions:  equation  (13)  reduces  to 

f:  =  xl~l(p  +  v,)iyt][(p+v3)IV,A.  (16- 

Note  that  6  —  0.  k'  =  0  because  of  the  neglect  of  bac' 
ground  density  inhomogeneity  (i.e..  L  -*  co). 

The  dispersion  relation  of  the  instability,  deter¬ 
mined  from  equation!  II),  is  obtained  by  setting/;  =  0 
which  leads  to  the  growth  rates  ( P )  of  the  instability 
to  be 

/»  =  (l/2){-(vl+v:)  +  {(vl— v,)2+4aiKl  K,]1’}. 

(17) 

Taking  p  -  0  in  (17),  we  derive  the  threshold  con¬ 
dition  for  exciting  the  instability 


eeB 


[nUoj^-Q,) 


=  (v,v:w,/fl,A-) 


(18) 


(/j)  Convective  loss 

We  next  examine  the  effect  of  the  finite  instability  zone, 
namely,  the  convective  loss  of  excited  modes  which  is 
assumed  to  dominate  over  the  coliisional  damping 
loss  in  determining  the  thresholds  of  the  instability. 
In  this  case,  equation  ( 1 5)  is  simplified  to  be 

/(:)  =  x:.(l  +c;/A:)  — (1/4) 

x  f/»(v,  H-v,)/K,  ^;]2-(I/2AJ).  (19) 


in  obtaining  (19).  we  have  also  assumed  that  0  »  KyA. 
The  zeros  of f(z)  are  then  found  to  be 

±&{4*-V\v:;{pHvi  +  v1y- 

+  2(V yV ,/A) '1  —  I } ' : 

=  ±-V  (20) 

These  are  the  positions  of  two  real  turning  points  for 
equation  (14).  which  has  the  following  solutions: 


*f"(z)  J(r)exp 

TO  =  r-(,)exPr 


for  r  >  z, 


for  z  <  r. 


well  behaved  at  z  ±  co.  Thus,  spatially  localized 
but  temporally  growing  modes  exist  provided  that  a 
positive  eigenvalue  pcan  be  found  from  the  eigenvalue 
equation  defined  by  the  Bohr-Sommerfeld  quan¬ 
tization  condition  (Landau  and  Lifshitz.  1965) 

(»i+|)ir=  \  (21) 


w1  •:»  =0.1,2 — stands  for  the  number  of  modes 

•  eigenfunction  *P(z)  between  the  two  turning 
z,  and  r,. 

•  integral  in  (21)  can  be  expressed  in  terms  of 
complete  elliptic  integral  for  numerical  analysis.  How¬ 
ever.  this  integral  can  be  much  simplified  for  z,/A  «  1. 
and  a  simple  approximate  form  of  (21)  results  as 

z,/A  ^  [(2m  +  l)/xA]' \  (22) 

Substituting  (20)  into  (22)  to  eliminate  z,  and  con¬ 
sidering  the  m  =  0  case  for  the  requirement  of  lowest 
thresholds,  we  obtain  the  growth  rate  of  the  insta¬ 
bility.  It  is 

/>  =  [2K,K,/A(k',  +  P':)]((aA)5/(l-i-aA)-]/2],:. 

(23) 

Setting  p  =  0  in  (23)  leads  to  the  following  expression 
for  the  threshold  fields  of  heater  wave 

«... 

jn((o0-Sle) 
where  u)p0  =  (4ire:/i0/m) 1 :. 


V,c 
‘Or,  0A 


(24) 


4.  DISCUSSION 

We  have  formulated  a  mechanism  whereby  whistler 
wave  and  Langmuir  wave  can  be  parametrically  ex¬ 
cited  by  an  HF  heater  in  the  overdensc  ionospheric 
heating  experiments  at  Tromso.  Norway.  The  required 
threshold  fields  of  heater  wave  are  derived  and 
examined  separately  for  two  cases,  wherein  either 
the  coliisional  loss  or  convective  loss  is  considered 
to  be  the  dominant  damping  process  of  the  pro¬ 
posed  instability.  A  quantitative  analysis  of  the 
instability  is  carried  out  here  w'ith  the  following 
parameters:  ojt,:'2n  =  4.04 MHz.  Q.J2n  =  1.4 MHz. 
V,  =  1.3  x  105ms' A  ~  150m.  v,  ~  500Hz.  and 

<0,0  Si  ft>„. 

It  is  found  that  both  calculations  of  the  threshold 
fields  from  (18)  and  (24)  have  nearly  the  same  mag¬ 
nitude,  i.e.  lvm'  ’.  Tnis  indicates  that  the  coliisional 
loss  and  convective  loss  play  equally  important  roles 
as  the  linear  damping  mechanisms  of  the  proposed 
instability.  The  effective  threshold  field  imposed  by 
both  damping  mechanisms  is,  therefore.  v-'2vm'!. 
The  growth  rate  of  the  instability  with  the  effective 
threshold  field  can  be  calculated  from  (23)  with  modi¬ 
fied  width  (A')  of  the  first  peak  of  the  heater  wave 
energy  distribution.  The  modified  width  is  given  by 
A'  =  A((l  -*)/(!  +  w)]1'-3  where  ,v  =  £'*/«o  and  c,h  is 
the  threshold  field  (given  in  (18))  imposed  by  the  coi- 
lisionai  loss  as  the  dominant  damping  mechanism. 
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Conceptually,  the  width  modification  results  from  the 
reduction  of  heater  wave  energy  by  collisions. 

The  effective  threshold  fields  (~  i.4vm  ')  of  the 
instability  can  be  exceeded  by  the  peak  field  intensity 
of  pump  waves  transmitted  from  the  Tromso  HF 
heater,  especially  after  the  swelling  effect  on  the  wave 
field  has  been  taken  into  account.  If  the  effective  field 
intensity  of  heater  wave  near  the  reflection  height  is 
assumed  to  be  3vm~'.  whistler  waves  with  fre¬ 
quencies  ranging  from  a  few  tens  of  kHz  up  to.  say. 
several  hundred  kHz  can  be  excited  in  a  few  seconds 
according  to  (23).  The  theory  predicts  that  Langmuir 


w'aves  excited  concomitantly  with  whistler  waves  at 
Tromso  may  have  a  much  broader  frequency  band¬ 
width  than  those  excited  by  the  parametric  decay 
instability.  The  frequency  bandwidth  of  excited  Lang¬ 
muir  waves  is  determined  by  the  wave  frequency  of 
excited  whistler  waves.  Stimulated  electromagnetic 
emissions  with  a  broad  frequency  band  were  indeed 
observed  in  the  Tromso  heating  experiments. 
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Abstract  The  temporal  evolution  of  HF-enhanced  plasma 
lines  (HFPLs)  is  investigated  theoretically.  The  HFPLs  refer 
to  the  radar  echoes  at  frequencies  near  the  sum  and  difference 
of  the  radar  frequency  and  the  HF  heater  wave  frequency. 
These  enhanced  spectral  lines  are  caused  by  backscatter  of 
radar  signals  from  plasma  waves  having  a  wavenumber  of  18 
m-L  A  nonlinear  theory  has  been  developed  to  explain  both 
the  intensity  overshoot  and  altitude  expansion  of  HFPLs  ob¬ 
served  at  Arecibo,  Puerto  Rico  [Djuth  and  Sulzer,  1989]  with 
good  agreement 

Introduction 

Parametric  instability  plays  an  important  role  in  ionospheric 
HF  heating  experiments  to  produce  various  nonlinear  plasma 
phenomena.  One  of  the  most  pronounced  signature  of  this 
instability  is  the.  constant  observation  of  HF-enhanced  plasma 
lines  (HFPLs)  at  Arecibo,  Puerto  Rico,  which  refer  to  the  en¬ 
hanced  backscatter  spectrum  of  the  Arecibo  430  MHz  radar 
signals  at  frequencies  near  430  MHz  ±  Thf  [e.g.  Carlson  et 
al.,  1972).  Here  fiff  is  the  frequency  of  the  HF  heater  wave. 
The  enhanced  radar  echoes  at  these  two  sidebands  provide 
useful  information  on  the  characteristic  features  of  those 
parametrically  excited  upgoing  and  clowngoing  plasma  waves. 
One  of  the  interesting  characteristics  of  the  observed  HFPLs  at 
Arecibo  is  the  plasma  line  overshoot  phenomenon  [Showen- 
and  Behnke,  1978;  Showen  and  Kim,  1978], 

This  overshoot  phenomenon  was  further  investigated  in  re¬ 
cent  experiments  on  the  temporal  evolution  of  HFPLs  by 
Djuth  et  al.  [1986].  The  experimental  results  at  high  HF 
power  (80  MW  ERP)  showed  that  the  HFPLs  exhibited  an 
initial  growth  for  a  few  tens  of  milliseconds  before  reaching  a 
maximum  intensity,  and  then  drastically  reduced  their  strength 
in  a  time  period  of  a  second.  An  interpretation  of  this  phe¬ 
nomenon  has  been  offered  in  Kuo  et  al.  [1987],  suggesting 
that  anomalous  damping  introduced  by  incoherent  scattering 
of  electron  orbits  by  the  excited  plasma  waves  results  in  a 
mode  competition  process.  It  is  found  that  anomalous  damp¬ 
ing  is  generally  larger  for  plasma  waves  with  greater  propaga¬ 
tion  angles  with  respect  to  the  geomagnetic  field.  Hence,  the 
plasma  waves  picked  up  by  the  Arecibo  430  MHz  radar  at  the 
propagation  angle  of,  approximately,  40°  can  be  suppressed 
by  those  pfopagatingat  smaller  angles  after  the  heater  is  on  for 
about  a  second.  Similar  mode  competition  process  leading  to 
pi asma  line  overshoot  has  also  been  considered  in  the  early 
work  by  Perkins  et  al.  [1974].  Jri  their-  formulation  of  the 
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nonlinear  damping  rate;  the  effect  of  geomagnetic  field  is  ne¬ 
glected.  Their  nonlinear  damping  mechanism  is  based  oh  res¬ 
onant  velocity  diffusion  of  charged  particles  in  the  beat  wave 
fields  (Le.  nonlinear  Landau  damping),  which  becomes  much 
less  effective  than  that-caused  by  orbit  (spatial)  diffusion 
mechanism  [Dum  and  -Dupree,  1970]  in  demagnetized 
plasma  if  there  is  no  specific  cyclotron  resonance  involved. 

In  the  heating  experiments  conducted  by  Djuth  and  Sulzer  at 
Arecibo  in  April;  1988,  the  HF  heater  was  operated  in  a  pulse 
mode  of  1  minute  on  and  9  minutes  off.  During  these  obser¬ 
vations,  the  HF  heater  was  operated  at  5.1  MHz  and  approxi¬ 
mately  80  MW  effective  radiated  power  was  transmitted. 
Radar  measurements  with  improved  altitude  and  temporal 
resolution  were  made  in  the  experiments  to  investigate  the  time 
evolution  of  both  the  intensity  and  location  of  HFPLs.  The 
observations  showed  that,  after  the  HF  heater  was  turned  on, 
the.  intensity  of  HFPLs  exhibited  overshoots.  During  the 
same  time  period,  the  altitude  interval  of  HFPLs  also  ex¬ 
panded  to  an  extent  greater  than  that  predicted  by  a  linear  the¬ 
ory.  The  upper  cutoff  altitude  (~  275.8  Km)  was  very  close 
to  the  reflection  point  of  HF  waves  in  the  ionospheric  plasma. 
While  the  location  of  HFPLs  never  exceeded  the  HF  reflection 
height,  it  could  generally  extend  downward  over  1  tb  2  kilo¬ 
meters  within  30  -  50  milliseconds  after  the  HF  heater  was 
turned  on  [see  Figure  1  of  Djuth  and  Sulzer,  1989], 

A  theory  is  offered  in  the  present  paper  to  explain  this  newly 
observed  downward  extension  of  altitude  interval  of  the 
HFPLs.  The  proposed  theory  is  based  on  the  "fake"  heating 
of  the  bulk  plasma  by  excited  Langmuir  waves  [e.g.  Segdeev 
and  Galeev,  1969,  pp  63-64],  In  addition  to  parametric  in¬ 
stabilities,  strong  turbulence  processes  involving  the  produc¬ 
tion  of  cavitons/solitons  have  been  suggested  as  possible 
sources  of  HF-enhanced  plasma  waves  in  the  Arecibo  iono¬ 
sphere  [e.g.  DuBois  et  al.,  1990;  Payne  et  al.,  1984].  The 
"fake  heating"  should  occur  regardless  of  the  source  of  in¬ 
duced  Langmuir  waves.  To  illustrate  the  effect,  we  will  use 
parametric  instability  theory  so  that  Langmuir  wave  intensities 
can  be  more  readily  calculated  and  the  modal  equations 
derived  in  our  previous  work  addressing  the  plasma  line  over¬ 
shoot  phenomenon  [Kuo  et  al.,  1987]  can  be  directly  incorpo¬ 
rated  in  the  present  study.  An  "apparent"  increase  in  the 
plasma  temperature,  which  is  proportional  to  tire  energy.den- 
sityof  excited  Langmuir  waves,  can  cause  the  lowering  of  the 
resonant  location  of  HFPLs  picked  up  by  tire  radar.  Thus, 
following  the  growth  ajidsaurrationof  Langmuirwaves,  an 
expansion  of  the  originatihg  altitude  intefval  of  HFPLs  is  de- 
tected.  :  V  "  - 

.  ,  .  ■  The^i^  NunjCTi^Jl^^jB:  . 

As  Langmuir  waves  are  excited  by  HF  heater  via  parametric 
decayinsrabtiityrtiieirwavefieldscanaffect  electronmotioh 
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and  cause  diffusion  in  the  velocity  space.  Such  a  diffusion 
process  can  be  described  by  a  quasilinear  equation  wherein  the 
diffusion  coefficient  can  be  separated  into  two  parts.  One  part 
arises  from  the  resonant  interaction  of  waves  with  electrons 
and  only  very  few  electrons  contribute  to  this  part  of  the  dif¬ 
fusion  coefficient.  By  contrast,  the  nonresonant  (or  adiabatic) 
wave-electron  interaction  involves  the  main  body  of  the 
plasma  distribution  in  the  plasma  oscillations. 

The  oscillatory  kinetic  energy  associated  with  the  excited 
Langmuir  waves  increases  with  the  field  intensity  of  the 
waves  and  the  main  body  of  the  plasma  appears  to  be  heated. 
This  effective  heating  of  electrons  by  Langmuir  waves  leads  to 
the  apparent  increment  in  the  electron  temperature  by  x  = 
2kEk2/47tne,  where  Ej^  is  the  spectral  intensity  of  Langmuir 
waves  and  ne  is  the  background  electron  density  [e.  g. 
Sagdeev  and  Galeev,  1969].  Consequently,  the  dispersion 
relation  of  Langmuir  waves  has  the  following  modified  form: 
©Z2  =  cope2+f2e2sin20+3k/2(Te+t)/ni<.  where  ©/,  (Ope,  k/,  fL, 
me  and  Te  have  the  standard  meaning  and  6  is  the  propagation 
angle  of  Langmuir  wave  with  respect  to  the  geomagnetic  field. 
Note  that  E/±  and  (0/  are  fixed  by  the  relations  t/±  = 
+(2Er±Jchf)  [Showen,  1979]  for  the  radar  detected  HFPLs, 
where  Er  is  the  wavevector  of  the  radar  signal  and  <0/  =  ©hf 
from  the  frequency  matching  condition  of  the  parametric  in¬ 
stability.  In  order  to  maintain  the  maximum  growth  rate  of  the 
parametric  decay  instability,  it  requires  that  ©pe^k^x/me 
remain  a  constant  so  that  the  sideband  of  the  instability  pro¬ 
cess  can  always  satisfy  the  dispersion  relation.  Consequently, 
the  matching  altitude  of  the  radar  detected  HFPLs  has  to  move 
downward  for  a  smaller  ©pe  because  plasma  heating  occurs  at 
altitudes  below  the  F  peak.  The  instability  of  Langmuir  waves 
associated  with  HFPLs  in  the  original  height  will,  however, 
continue  to  grow  as  long  as  the  instability  thresholds  are  still 
exceeded  by  the  heater  wave  field.  Then,  the  lowering  of  the 
HFPLs'  matching  altitude  actually  represents  the  expansion  of 
the  originating  altitude  interval  of  the  detected  HFPLs.  This 
process  of  lowering  the  matching  altitude  of  HFPLs  can  be 
described  by  d©pe2/dt  =  -3(k/2/m«)  dx/dt. 

A  linear  density  profile  is  assumed  for  the  heated  iono¬ 
sphere,  i.e.,  (0^2= ©pco^l+x/L),  where  L  is  the  linear  scale 
height  of  the  ionospheric  density  and  ©pc2  is  the  electron 
plasma  frequency  at  the  initial  location  of  the  concerned 
HFPLs.  Then,  the  above  equation  can  be  written  as 

dx/dt  =  -  (3k/2L/me©pco2)dx/dt  (1) 

showing  the  temporal  change  of  the  height  of  the  detected 
HFPLs.  Let  le  =  4jtsin8jdkk2Eke2  represent  the  spectral  in¬ 
tensity  of  Langmuir  waves  propagating  at  the  same  angle  0 
with  respect  to  the  geomagnetic  field.  Thus,  x  is  defined  to  be 
Jo  *  IedS^jnte,  where  0O  is  the  averaged  maximum  spectral 
angle  of  the  excited  Langmuir  waves  [Kuo  et  al.,  1987]. 
Substituting  this  expression  for  x  into  (1)  yields 

dx/dt  =  -(3k^L/4jtnemetOpeo2)fI8.d0o/dt+Jo^lWdt  d0]  (2) 

where  0O  and  le  are  governed  by  the  rate  equations  derived  in 


our  earlier  work  [Kuo  et  al.,  1987]  on  the  study  of  HFPLs 
overshoot  phenomenon.  They  arc 

dlf/dt  =  [fiecosQ  —  pe  —  2de  —  Ttie1/2]l0  (3) 

and 

d0o/dt  =  -(bsinSQoCos^Oo/tasinSOo  +  ble„sin50ocos50o 

+  (2  -  5sin20o/cos20o)(aeocos0o  -  pe.  -  tTe„Ie01/2]} 

Jo0°[aecds0  -  Pe  -  2de  -  tfle’^Iesin^cosZOdO  (4) 

where  b  =  a6c2/32pi5B02©02,  a  -  a2/Pi,  ae  =  2aVk0,  pe  = 
2(Po+Piko),%  =  (Ak/ko)3/2ekoi/2/(47tsin0mcmi©oc1)1/2,  and 
de  =  [hac2ko2sin20/2Bo2©o2]io8'l0’sin20’cos0’d0’ ;  0O  is 

determined  via  ae.cos0o-  Pe.  -  2de.  -  fjeje.  =  0  in  equation 
(3)  because  plasma  waves  with  this  propagation  angle  have 
zero  nonlinear  growth  rate;  a  =  6EJ1  (metni©^!/2,  ko  =  ct2 
cos20/4Pi2  and  po = Vej/4 ;  2Eo,  ©o,  c,,  and  vc;  are  the  ampli¬ 
tude  of  the  HF  heater  wave,  the  heater  wave  frequency,  the 
ion  acoustic  speed,  and  the  electron-ion  collision  frequency, 
respectively.  Pi  defined  by 

-3/2]  is  proportional  to  the  ion  Landau  damping  rate ;  Ak/k0 
given  by  (^©oPo/koV^2  is  the  bandwidth  of  spectral  lines 
contributing  to  the  anomalous  damping  rate  of  the  Langmuir 
cascading  process  ;  h  =  2ko/Akflco  and  vle  is  the  electron 
thermal  speed  defined  by.fiyrae)1/2. 

The  three  coupled  equations  (2)-(4)  describe  the  temporal 
evolution  of  x,  le  and  0O  which  correspond  to  the  altitude  of 
HFPLs,  the  spectral  intensity  of  excited  plasma  waves,  and 
the  maximum  propagation  angle  of  plasma  waves,  respec¬ 
tively.  For  a  numerical  analysis,  we  rewrite  these  equations  in 
terms  of  the  dimensionless  variables  x  =  x/L,  %  =  2pot,  and  Ig 
=  laAtaieTe  as  follows: 

dx/'d|  =  -Ape.(d0o/dt)+/o0‘,(dIa/dt)d0  (5) 

dle/d^  =  [ai2cos29  -  1  -BcosSOsin^/oMesin2©’ 

cos0’d0’  -  Ccos-^sin-wOle  1/2)10  (6) 

d0o/d^  =  -  [  Bsin38ocos40o/[2ai2sin20o  +  Ble, 
sin^oCos^o  +  (2  -  5sintyo/cos2Q0) 
(aftosty  - 1)  -  C(5/2  -  7510200/005200) 

cos-20osin-W20oIe.W2]  ^’(die/d^sm^osOdO.  (7). 

where  ai  =  a/a*,  a*  =  2(P0Pj)»/2,  A  =  Sk^v^/©^2,  b  = 
(^aisPoSv^mpeoa/p^Qe2©^,  and  C  =  Pi^Wp^ZOai2 
Pb2V(e2.  Using  the  parameters  relevant  to  the  Aretibo  heating 
experirnerittuv*  =  13xl05rn/sec,c,  =  O^ 

=  1100?K,  ©o/2x  =  5.1  MHz,  Xjt  =  2u/kjr=  35  cm  (for  the 
plasma  mode  detected  by  the  430  MHz  radar),  Vd  =  500  Hi. 
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and  aj  =  2.0,  we  then  have  0O  =  125  H2.  p1  -  31.8  m/sec,  A 
=  1/50,  B  =  120  and  G  =  0.325.  .Thus,  the  time  dependence 
ofithe  matching  altitude  and  intensity  of  HFPLs  can  be  deter¬ 
mined  by  integrating  equations  (5H7). 

Shown  in  Figure  1  is  the  numerical  result  of  temporal  evo¬ 
lution  of  the  spectral  intensity  I400  of  the  Langmuir  waves 
simulating  that  detected  by  the  Arecibo  430  MHz  radar  at  the 
angle  of  40°  with  respect  to  the  geomagnetic  field.  The  ob¬ 
served  overshoot  phenomenon  of  HFPLs  is  successfully  re¬ 
produced  by  our  model.  Further,  the  calculated  evolution  of 
the  HFPLs’  height  is  presented  in  Figure  2,  showing  that  the 
matching  altitude  of  HFPLs  moves  downward  as  the  instabil¬ 
ity  grows.  The  evolution  of  the  “fake”  temperature  increment 
x  is  also  shown  in  the  same  figure.  As  mentioned  before,  this 
represents  the  expansionof  the  originating,  region  of  the 
HFPLs  because  the  Langmuir  waves  responsible  for  HFPLs 
are  excited  not  only  from  their  original  location  but  also  ex¬ 
tending  to  a  lower  altitude  due  to  "fake  heating"  of  electrons 
by  the  Langmuir  waves.  The  calculated  distance  of  lowering 
is  0.05  times  of  the  ionospheric  scale  height  (L)  after  the  HF 
heater  is  turned  on  for  50  milliseconds.  The  ionospheric  scale 
height,  L,  deduced  from  the  experiments  is  about  22  kilome¬ 
ters,  leading  to  the  expansion  of  originating  altitude  interval  of 


Fig.  1  Temporal  evolution  of  the  intensity  Lo°  of  HFPLs 
detected  by  Arecibo’s  430  MHz  radar. 


Fig.  2  Normalized  matching  altitude  x-.and  the  fake 
temperature  x  vs.  time  L 


HFPLs  to  be  1.1  kilometers.  These  theoretical  results  agree 
well  both  qualitatively  and  quantitatively  with  the  observations 
[Djuth  and  Sulzer,  1990]. 

In  the  formulation  of  present  theory,  mode  Competition 
mechanism  is  included  for  overshoot  while  cascading  mecha¬ 
nism  is  included  for  the  saturation  of  the  plasma  lines.  Only 
when  the  first  mechanism  dominates  over  the  second  mecha¬ 
nism,  overshoot  can  occur.  Therefore,  it  is  useful  to  find  out 
the  relationship  between  overshoot  phenomenon  and  the  HF 
intensity.  The.'dependence  of  die  overshoot  intensity  and 
overshoot  time  on  the  pump  in  tensity  are  concerned. 

Let  ti,  t2  and  t3  represent  the  times  for  the  plasma  line  in¬ 
tensity  to  reach  its  peak,-then  drop  to  one-half  and  then  one- 
tenth  of  its  peak  value,  respectively.  Note  that  tj  approaches 
infinity,  if  this  plasma  line  does  not  overshoot;  The  depen¬ 
dence  of  these  characteristic  times  on  the  heater  intensity  is 
displayed  in  Figures  3(a)  and  3(b).  Shown  in  Figure  4  is  the 
dependence  of  the  peak  intensity  of  HFPLs  on  the  heater  in¬ 
tensity.  The  results  show  that  both  overshoot  times  and  peak 
intensity  of  HFPLs  decrease  generally  with  the  heater  inten¬ 
sity.  It  is  also  shown  that  the  overshoot  phenomenon  occurs 
oniy  when  the  heater  intensity  exceeds  the  threshold  of  the 
parametric  decay  instability  by  a  value  ct]  ~  1.42.  These  pre¬ 
dictions  await  to  be  corroborated  in  the  future  experiments. 


fig.  3  Funtiooal  depertdence  of  the  times  (t^tiXfor  HFPLs’ 
intensity  to  reach  its  peak  (tj),  to  drop  to  one-half  (12)  and 
one-tenth  O3)  of  its  peak  value  on  the  normalized  pump 
intensity  aj. 

(a)  tivsctj. 

(b) t2andt3vsct{. 
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Summary 

Our  theoretical  model  [Kuo  et  al.,  1987]  developed  origi¬ 
nally  for  explaining  the  intensity  overshoot  of  the  HFPLs  has 
been  extended.  The  extension  of  the  model  is  based  on  the 
fact  that  rionresonant  heating  of  plasma  electrons  by  tire  para¬ 
metrically  excited  Langmuir  waves  can  yield  an  apparent  in¬ 
crement  in  electron  temperature.  This  fake  heating  of  elec¬ 
trons  introduces  an  extra  term  into  the  dispersion  relation  of 
Langmuir  waves.  The  modification  of  the  dispersion  relation 
moves  the  matching  location  of  the  HFPLs  to  a  lower  altitude 
where  the  plasma  density  is  less  than  that  at  their  initial  loca¬ 
tion.  This  model  has  successfully  reproduced  both  phenom¬ 
ena  of  the  intensity  overshoot  and  the  downward  expansion  of 
the  originating  altitude  interval  of  HFPLs  in  good  agreement 
with  the  observations.  The  numerical  results  also  show  that 
both  the  overshoot  time  and  peak  intensity  decrease  generally 
with  the  pump  intensity  and  overshoot  occurs  only  when  the 
pump  intensity  exceeds  the  threshold  intensity  of  the  paramet¬ 
ric  decay  instability  by  a  factor  aj  -  i.42. 

It  is  noted  that  thermal  filamentation  instability  can  generate 
large  scale  density  perturbation  (>  1  Km  scale  length)  which 
may  also  cause  the  shift  of  the  plasma  line  altitude..  However, 
two  facts  stand  on  the  way  to  rule  out  the  possibility  cf  such  a 
cause  mechanism  for  the  observed  phenomena.  One  is  its 
growth  time.  It  takes  at  least  a  few  seconds  for  the  develop¬ 
ment  of  large  scale  density  irregularities  through  thermal  fila- 
mentatior.  instability  [Kuo  and  Schmidt,  1983}.  The  other  one 
is  the  filamentation  direction  which  is  perpendicular  to  rather 
than  within  the  magnetic  meridian  plane  for  an  O-mode  pump 
wave  [Kuo  and  Schmidt,  1983]. 
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The  coupling  of  hydromagnetic  Alfve'n  waves  is  studied  numerically  in  a  dipole- 
field  model  of  the  magnetosphere.  The  two  coupled  hydromagnetic  equations 
derived  by  Radoski  are  solved  as  an  implicit  boundary-value  problem,  namely 
the  boundary  conditions  at  the  magnetopause  are  determined  self-consistently. 
Thus  the  calculated  wave-field  distribution  inside  the  magnetosphere  can 
match  all  known  linear  characteristic  features  of  the  stormtiine  Pc5  waves 
observed  on  14/15  November  1979  from  satellites  A  set  of  proper  boundary 
oonditions  is  found,  excellent  agreement  between  the  numerical  results  and 
observations  is  demonstrated.  Based  on  the  very  limited  spatial  coverage 
\L  %  6-6  and  within  a  latitudinal  region  (—10°.  10°)),  of  the  data  provided 
by  the  satellites,  the  theoretical  model  can  successfully  reconstruct  the  global 
micropulsations  in  the  magnetosphere  and  identify  the  source  regions  of 
hydromagnetic  waves. 


1.  Introduction 

Recent  observations  of  a  stormtime  Pc5  evimt  during  14/ lo  November  1979. 
near  geosynchronous  orbit  ,'Rigbie  etui.  1982:  Takahashi.  Higbie  -k  Baker  1985) 
provide  comprehensive  data  on  the  characteristic  features  of  the  stormtime  Pc5 
waves,  where  Pc5  waves  appear  as  geomagnetic  fluctuations  having  their 
oscillation  periods  in  the  range  from  !50  to  600  s.  This  event  was  characterized 
by  long-duration  (50  hi  highly  compressional  magnetic  perturbation,  and 
modulation  in  energetic  particle  fluxes.  The  azimuthal  phase  velocity  Jl^i  = 
4-14  km  s_1  and  azimuthal  wavenumber  m  =  20-120  were  measured.  However, 
the-  field-aligned  structure  of  the  stormtime  Pc5  waves  was  not  reported  until 
recently  by  Takahashi  el  al.  (1987).  .Magnetic  field  data  from  the  four  satellites 
SCATHA  (P78-2).  GOES2.  GOES3  and  GEOS2  have  been  analysed  to  examine 
the  wave  structure  along  the  ambient-  magnetic  field  line  (L  %  6-6)  within  a 
latitudinal  region  (—10°,  10°).  Five  distinct  features  have  been  identified.  If  has 
been  observed  that  the  compressional  component  of  the  wave  has  a  node  near 
the  equator,  offering  evidence  of  an  antisymmetric  standing-wave  structure, 
while  the  transverse  component  has  a  symmetric  standing-wave  structure  with 
an  antinode  near  the  equator.  Moreover,  nodes  for  the  transverse  components 
at  several  degrees  off  the  dipole  equator  have  also  been  observed.  They  indicate 


160 


31. 11.  Whang.  8.  P.  Kuo  and  31.  C.  Lee. 

that  these  components  of  the  magnetic  perturbation  have  a  very  short  parallel 
wavelength  about  3 Re.  which  is  much  shorter  than  30 lir  for  the  fundamental 
mode  of  the  shear  wave  at  L  %  6-6.  For  the  wave  frequency  of  0002  Hz.  the 
field-aligned  phase  speed  of  the  Pc5  wave  is  estimated  to  be  about  40  km  s'1,  an 
order  of  magnitude  lower  than  the  estimated  Alfven  sjieed  (Takahashi  et  nl. 
1985).  It  has  been  shown  that  the  two  transverse  components  of  the  wave  have 
different  wavelengths,  as  evidenced  by  the  result  that  their  nodes  are  located  at 
different  magnetic  latitudes.  One  possible  cause  leading  to  this  result,  suggested 
by  Takahashi  et  al.  (1987),  is  the  inhomogeneity  of  the  magnetosphere  along  the 
field  line.  The  relative  phase  among  the  magnetic  field  components  is  either  0°. 
2: 90°  or  180°.  and  changes  by  180°  as  the  observing  satellites  move  across  the 
nodes  of  the  wave.  The  compressions’  component  oscillates  90°  out  of  phase 
with  the  azimuthal  component.  It  oscillates,  however,  either  in  phase  or  180° 
out  of  phase  with  the  radial  component.  An  interpretation  of  these  relationships 
based  on  wave  propagation  in  an  inhomogeneous  medium  is  also  provided.  The 
last  of  the  identified  features  is  the  presence  of  harmonies.  This  is  manifested  by 
the  wave  form  of  the  compressional  magnetic*  field  component,  which  oscillates 
in  time  with  different  heights  for  two  consecutive  peaks.  This  feature  is 
associated  with  the  nonlinearity  of  the  piastna. 

Many  mechanisms  have  been  suggested  to  explain  the  excitation  of  the 
stormtime  Pco  waves.  These  include  the  drift  mirror  instability  (Hasegawa 
I960),  bounce  drift  resonance  excitation  of  L’LF  waves  (Southwood  1976).  drift 
compressional  instability  (Hasegawa  1971 :  Xg.  Patel  &  Chan  1984)  tnd 
coupling  between  a  drift  mirror  wave  and  a  shear  Alfven  wave  (Walker  ai. 
1982).  Recently.  Cheng  &  Lin  (1987)  performed  a  comprehensive  eigenmode 
analysis  for  both  the  drift  mirror  and  drill  compressional  instabilities.  The 
results  show  that  two  of  the  five  features  deduced  from  the  observations,  the 
antisymmetric  standing-wave-  structure  and  the  phase  relations  among  the  field 
components,  can  be  explained  by  the  eigenmode  structures  of  the  drift  mirror 
instability  in  the  dipole  magnetic  field  configuration.  However,  the  difficulty  in 
using  the  drift  mirror  instability  to  interpret  stormtime  Pc5  waves  comes  from 
the  instability  condition.  I'sing  the  “BCATHA  ion  data  (17  eV-300  keV). 
Takahashi  *7  al.  {1987)  have  shown  that  the  instability  condition  cannot  be 
satisfied.  They  have  also  shown  that  all  the  existing  theories  have  limited 
success  in  explaining  the  results  of  the  event.  It  is  therefore  recognized  that  a 
new  theoretical  development  to  explain  the  excitation  mechanism  of  stormtime 
Pco  waves  is  desirable. 

Information  cn  the  spatial  structure  of  the  waves,  both  radial  and  latitudinal 
variations,  seems  to  play  a  key  role  in  guiding  the  theoretical  development. 
However,  only  very  limited  spatial  coverage  of  the  data  provided  by  the 
satellites  is  available.  In  the  present  work- we  study  theoretically  the  coupling, 
and  propagation  of  hvdromagnetie  waves  in  the  dipoie-field  model  of  the 
magnetosphere.  The  set-  of  coupled  hvdromagnetie  wave  equations  (Radoski 
1967)  is  solved  numerically.  -The  solution  satisfies  the  perfect -conductor 
boundary  conditions  at  the  ionosphere,  together  with  the  imposed  wave 
structures  drawn  by  Takahashi  el  al.  (19S7)-from  the  magnetospherie  storm 
-events  during  14/15  November  1979.  The  local  data  obtained  by  satellites  are 
then  used  ;to  reconstruct  the  global  wave  structures  of  the  event  in  terms  of 
hvdromagnetie  waves  in  the  magnetosphere.  The-  boundary  conditions  at-  the 
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magnetopause  are  determined  self-consistentlv  so  that  the  spatial  dependence 
of  the  solutions  along  the  L  =  0  0  line  matches  all  four  linear  features  drawn 
from  the  satellite  data  (Takahashi  el  nl.  1987). 

Section  2  describes  the  dipole-field  model  of  the  magnetosphere  and  the 
coupled  hvdromagnetic  wave  equations  used  in  the  following  numerical 
analysis.  The  numerical  results  and  a  comparison  with  the  satellite  data  are 
presented  in  §3.  A  summary  of  the  present  work  and  conclusions  are  given  in 
§4.  The  method  of  numerical  analysis  is  briefly  described  in  an  Appendix. 


2.  Governing  equations 

It  is  generally  believed  that  geomagnetic  mieropulsations  are  associated  with 
hvdromagnetic  waves  that  can  be  excited  through  different  types  of  plasma 
instabilities  in  the  magnetosphere  or  on  its  boundary,  or  in  the  solar  wind.  In 
general,  the  sources  of  these  plasma  instabilities  can  be  considered  to  be 
localized  in  comparison  with  the  range  of  the  magnetosphere.  Therefore  only  a 
source-free  hvdromagnetic  wave  equation  will  be  used  in  the  following  to  study 
the  propagation  of  hydromagnet ic  waves  in  the  magnetosphere.  The  MHD 
equations  are  used  for  the  analysis  of  hvdromagnetic  perturbations  in  the 
magnetosphere,  which  can  be  considered  to  be  an  ideal  MHD  system.  It  is  then 
straightforward  to  convert  the  MHD  equations  into  the  ideal  MHD  wave 
equation,  which  is  a  single  second-order  time-derivative  equation  for  the  wave 
electric  field  E: 

~  =  Ax{AxlVx(VxE)]}.  (1) 


where  A  =  B0/(47rp)  is  the  Alfven-wave  velocity. 

An  orthogonal  dipole  co-ordinate  system  is  employed  to  describe  the  wave 
propagation  in  the  magnetosphere.  The  unit  vectors  are  along  the  principal 
normal  to  the  field  line  (v).  parallel  to  the  field  line  (ji)  and  in  the  azimuthal 
direction  (<|<)  respectively.  More  specifically,  the  co-ordinates  are  given  by  v  = 
(sin'-0)/r.  which  is  constant  along  a  dipolE-field  line,  ft  =  *(cos<?)/r\  which  is 
constant  along  an  orthogonal  trajectory  of  the  dipole  field  lines,  and  which 
is  the  ordinary  azimuthal  spherical  polar  co-ordinate;  f  =  r/lle  is  the  normalized 
radius  in  spherical  co-ordinates. 

Having  a  time-harmonic  dependence  of  the  form  e~,u .  and  an  oscillating 
longitudinal  variation  of  the  form  e lm*.  the  two  scalar  components  of  (I)  in  the 
dipole  model  of  the  magnetosphere  become  (Radoski  1967) 


H, 


0 •  Wi  .  0 

Here  r  is  related  to  r  and  /i  through  the  relation  vr+fi2?'  =  1, 
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with  E,,  the  toroidal  component  of  the  wav  >  electric  field. 

=  (7*4  E^, 

with  Ep  the  poloidal  component  of  the  wave  electric  field,  and  p  is  the  plasma 
mass  density,  with  the  spatial  dependence  p  =  p0(l/r)n,  where  p0  is  the  plasma 
mass  density  on  the  earth  and  n  is  an  adjustable  positive  density  index.  Since 
a  dipole  magnetic  field 


_  M0(i>+4,Pr*\t 


is  assumed, 


or  _  fwRp,2  ■inp0ri~n 
V  +  4fi27* 


Equations  (2)  and  (3)  describe  two  basic  types  of  hvdromagnetic  waves.  One 
is  called  the  guided  or  toroidal  mode;  it  channels  its  energy  only  along  the 
magnetic  field  lines  of  forces.  The  other,  called  the  isotropic  or  poloidal  mode, 
propagates  equally  in  all  directions  at  the  Alfven  speed  (Radoski  1976).  In 
general,  these  two  modes  are  coupled  to  each  other  through  the  inhomogeneities, 
as  shown  by  the  coupling  terms  on  the  right-hand  sides  of  (2)  and  (3).  However, 
for  the  special  case  of  axial  symmetry,  m  =  0.  (2)  and  (3)  become  decoupled,  and 
thus  these  two  modes  are  separate.  In  this  case  (2)  and  (3)  represent  the 
eigenequations  of  the  toroidal  and  poloidal  modes  respectively.  Along  the 
magnetic  field  line  (£).  the  phase  velocity  of  the  toroidal  mode  has  magnitude 
equal  to  the  local  Alfven  speed.  Hence  the  toroidal  mode  is  sometimes  also 
called  the  Alfven  mode.  On  the  other  hand,  the  poloidal  mode  propagates 
isotropically,  with  total  phase  speed  equal  to  the  Alfven  speed.  While  (B^.E^ 
represent  the  field  components  of  the  toroidal  mode.  (B  ,Bv>Ep)  represent  those 
of  the  poloidal  mode.  The  corresponding  fluid  velocity  perturbations  are  found 
to  be  in  the  <|*  and  v  directions  respectively.  V^,  the  response  of  E,„  is 
perpendicular  to  the  wave  propagation  direction,  which  is  in  the  (p..  4*)  plane 
(for  the  m.  =  0  case),  while  V,..  the  response  of  E ^  has  a  component  along  the 
propagation  direction.  On  the  basis  of  the  shear  compressional  features  of  the 
velocity  responses,  the  shear  (or  slow)  and  the  compressional  (or  fast)  modes  are 
also  called  the  toroidal  and  poloidal  modes  respectively. 

Since  the  toroidal  mode  only  channels  its  energy  along  the  magnetic  field,  it 
resonates  only  with  the  local  field  line,  and  its  oscillating  frequency  is.  m 
general,  characterized  by  the  length  of  the  associated  geomagnetic  field  line. 
Hence  periods  of  geomagnetic  micropulsations  associated  with  such  field-line 
resonant  modes  vary  with  the  L  value  of  the  field  line.  It  is  believed  that  such 
field-line  resonant  modes  are  directly  responsible  for  the  geomagnetic 
micropulsations  (Dungey  1954:  Chen  &  Hasegawa  1974:  Radoski  1974:  Chen  & 
Cowley  1989).  On  the  other  hand,  the  energy  of  the  poloidal  mode  can  flow 
across  the  field  line  and  fill  the  volume  of  the  magnetosphere.  The  characteristics 
of  the  poloidal  mode  depend  on  the  size  of  the  magnetosphere,  which  is  a  feature 
of  the  cavity  mode  (see  e.g.  Kuo.  Lee  &  Wolfe  1987).  Without  coupling 
( m  =  0),  the  field-line  resonant  mode  (toroidal  mode)  will  not  respond 
to  perturbations  starting  outside  the  field-line  resonance  region.  However, 
perturbations  in  the  magnetosphere,  in  general,  stem  from  the  solar  wind, 
which  is  outside  the  magnetosphere.  Therefore  the  study  of  the  general  case 
(m  =#  0),  with  coupling  between  toroidal  and  poloidal  modes,  is  clearly  crucial 


1(33 


Global  micropulsations 


to  the  understanding  of  dynamical  processes  such  as  magnetospheric 
substorms  resulting  from  the  solar-wind-magnetosphere  interaction. 

For  solving  (2)  and  (3).  the  boundary  conditions  must  be  specified.  We 
assume  that  the  magnetosphere  extends  from  L  =  2  to  />  =  10.  the  ionosphere 
is  a  perfect  conductor,  and  the  tangential  components  of  the  electric  field  at  the 
inner  boundary  are  zero.  These  assumptions  lead  to  the  fixed  boundary 
conditions  P(v,p  =  ±0-25)  =  0  =  T(v,/i  =  ±0-25) .  where  P  =  ie^  and  T  =  e,.. 
The  boundary  conditions  at  the  magnetopause.  P(v  =  OT./r)  and  Ttv  =  0*1 . //). 
are  generally  imposed  by  an  external  source  such  as  the  solar  wind.  They  can 
thus  be  assumed  to  be  any  functions  of  p  varying  with  the  physical  situation. 
For  good  agreement  between  numerical  results  and  recently  reported 
observational  data  (Takahashi  et  al.  1987).  the  boundary  conditions  at  i>  =  0-1 
are  found  to  be  of  the  form 


P(v  =  0-1.//)  =  P0(p)  =  . jin 


m. 


(-t) 


T(v  =  0-I.h)  =  O.  (o) 

where  P0(/i)  is  a  spatially  distributed  sinusoidal  perturbation,  q  is  an  odd  integer 
less  than  nine,  and  k  is  the  harmonic  number  of  the  sinusoidal  perturbation.  For 
practical  reasons,  the  toroidal  mode  at  the  outer  boundary  is  cut-off.  i.e.  T  —  0 
in  the  frequency  range  of  the  observed  perturbations.  Using  the  relations 
derived  from  Faraday's  law. 


the  numerical  solutions  for  the  magnetic  field  can  be  determined  from  the 
solutions  for  the  eiectric  field. 


3.  Numerical  results  and  comparison  with  observations 

In  the  appendix  our  numerical  method  is  described  and  the  coupled  equations 
(2)  and  (3)  are  converted  into  two  difference  equations  (A  1)  and  (A  2)  for 
numerical  analysis.  Equations  (A  1)  and  (A  2}  are  then  solved  to  obtain  the 
numerical  solutions  of  the  coupled  hydromagnetic  wave  equations  (2)  and  (3). 
The  local  data  acquired  by  four  satellites  during  the  occurrence  of  the  magnetic 
substorm  event  on  14/15  November  1979  will  be  used  to  prescribe  conditions 
to  be  matched  by  the  numerical  solutions.  These  data  on  magnetic  fluctuations 
were  recorded  with  flux-gate  magnetometers.  The  geostationary  satellites 
GOES2,  GOES3  and  GEOS2  were  located  at  the  geographic  equator,  with 
geographic  longitudes  of  107°  W,  135°  W  and  15°  E.  which  correspond  to 
magnetic  latitudes  of  9°  X.  5°  N  and  1°  X,  by  taking  centred  dipole  graphic  co  • 
ordinates.  GCATHA  had  an  elliptical  orbit  with  an  orbital  period  of  23-6  h.  The 
apogee  and  perigee  were  7-8/?f  and  5-3 Rt  respectively,  and  the  inclination  of  the 
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Fiovre  I.  Xoimalized  electric  fields  e,-,  and  e,  at  the  magnetopause  with  k  =  4  and  q  —  o. 


orbit  was  7-8°.  During  the  course  of  the  14/15  November  event,  its  magnetic 
latitude  varied  between  —7°  and  0°.  The  details  of  the  experiments  were 
described  by  Fennell  (1982)  with  SCATHA.  by  Grubb  (1975)  with  GOES2  and 
G0ES3.  and  by  Knott  (1982)  with  GE082. 

A  normalized  quantity 

<ff> 

<#>  +  </$>  +  </$> 

was  introduced  in  Takahashi  e.t  at.  (1987)  to  describe  the  relative  amplitudes  of 
each  component  of  the  magnetic  field,  which  varies  with  magnetic  latitude. 
Here  a  =  »\  ^  or  //.  and  (B'l)  is  the  mean-square  amplitude  of  the  magnetic  field 
component,  which  was  obtained  by  integrating  its  power-spectral  density  over 
the  frequency  range  of  the  wave.  In  the  same  spectral  analysis  the  relative 
phase  of  the  magnetic  field  components  was  also  calculated.  The  relative  phase 
of  component  x  with  respect  to  component  /?  is  denoted  by  <j> x/j. 

The  best  match  of  numerical  results  with  observational  data  is  found  when 
the  harmonic  numbers  of  the  sinusoidal  perturbation  of  the  poloidal  mode  at 
the  magnetopause  as  given  by  (4)  are  M  =  20.  k  =  4  and  q  =  5.  The  normalized 
wave  electric  field  at  the  magnetopause  with  the  boundary  condition  (4)  has  the 
form  shown  in  figure  1.  Equations  (A  1)  and  (A  2)  are  then  solved  to  obtain  the 
field  distributions  in  the  entire  magnetosphere.  The  consistency  of  the  numerical 
results  with  observations  is  checked  as  follows.  The  numerical  results  for  the 
magnetic  field  around  the  equator  in  the  magnetic  latitudinal  range  —  11°  <  A 
<  1 1°  at  a  radial  distance  L  ss  6-6,  i.e.  v  =  0T515.  are  extracted  for  comparison. 
The  ionospheric  boundary  is  chosen  to  be  about  l-84i?#  in  our  mode).  Following 
the  same  procedure  as  Takahashi  el  al.  (1987)  for  data  analysis,  the  normalized 
amplitude  and  the  relative  phase  are  determined  by  the  numerical  solutions  of 
the  coupled  hvdromagnetic  wave  equations  as  shown  in  figures  2  and  3.  These 
results  can  now  be  compared  directly  with  the  observational  data  summarized 
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FiouRE  3.  Intercomponent  cross-phase  of  the  wave  given  from  numerical  calculations. 

in  figures  5  and  6  of  Takahashi  et  al.  (1987.  hereinafter  referred  to  as  ref.  I). 
Excellent  agreement  is  found  between  our  numerical  results  (i.e.  figures  2  and 
3)  and  the  observations  (i.e.  figures  5  and  6  of  ref.  I).  ’ 

As  stated  in  §1,  five  distinct  features  have  been  identified  from  the 
observational  data.  These  features  are  the  basis  for  the  further  test  of  the 
theoretical  results.  First,  the  property  of Bz  having  short  parallel  wavelength  of 
a  few  earth  radii  (about  3 Re)  as  shown  in  figure  5(c)  of  ref.  1  is  reproduced  by 
the  numerical  result  shown  in  figure  3(c)  for  B/r  As  shown  in  figures  5 (a)  and  ( b ) 
of  ref.  1,  the  nodes  of  Bx  and  By  appear  at  different  locations,  indicating  that 
these  two  wave  components  have  different  parallel  wavelengths.  Such  a 
multiple-wavelength  phenomenon  also  appears  in  the  numerical  results  as  can 
be  seen  from  figures  2(a)  and  (6)  on  comparing  the  locations  of  the  nodes  of  B,. 
and  Bf.  The  antisymmetric  standing-wave  structure  found  in  figure  5(c)  of 
ref.  1  is  also  reproduced  by  the  numerical  results  shown  in  figure  2  (c). 

The  relative  phases  between  each  pair  of  magnetic  field  components  are 
obtained  by  calculating  the  latitude  dependence  of  the  relative  amplitude  of  the 
pair.  The  relative  amplitude  between  the  v  and  5 Ir  components,  for  example,  is 
defined  as 

where  a  =  v  or  xjr.  The  relative  phases  for  the  three  different  components  are 
evaluated  and  displayed  in  figure  3.  They  agree  quite  well  with  the  observational 
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Figure  4.  Normalized  magnetic  field  of  radial  component  (m  =  20.  k  =  4.  7  =  ■'>)■ 


results  shown  in  figure  fi  of  ref.  1.  Both  the  numerical  and  observational  results 
show  that  Ba  and  Br  oscillate  180°  out  of  phase  for  6  >  5°  and  in  phase  for 
0  <  5°:  B/t  and  are  out  of  phase  by  +90°. 

The  last  feature  deduced  from  the  satellite  data  is  the  frequency  doubling 
often  observed  near  the  magnetic  equator.  The  second-harmonic  waves 
observed  by  GEOS2  and  GEOS3  appear  along  the  field-line  direction  near 
magnetic  latitude  A  x  0-5°  and  along  the  azimuthal  direction  at  magnetic 
latitude  A  %  4-7°  respectively.  It  is  generally  believed  that  the  harmonic  is  due 
to  the  nonlinear  effect  of  the  plasma.  The  present  study  only  considers  MHD 
waves  in  a  linear  magnetospheric  plasma,  so  the  nonlinear  harmonic  generation 
process  cannot  be  revealed  in  the  numerical  results.  Since  the  nonlinearity  of 
the  magnetospheric  plasma  is  weak  because  of  the  low  plasma  density,  its  effect 
on  the  results  of  the  presently  extrapolated  global  reconstruction  is  expected  to 
be  small. . 

As  mentioned  in  §  1.  the  satellites'  orbits  cover  only  a  small  and  specific  region 
of  the  magnetosphere.  Hence  the  collection  of  the  data  by  them  is  bounded  in 
a  local  region.  By  contrast,  our  numerical  results  are  applicable  to  the  entire 
magnetosphere.  Therefore  the  present  -work  provides  a  way  to  use  iocal 
information  to  reconstruct  the  global  hvdromagnetic  waves  in  the  mag¬ 
netosphere.  Consequently,  detailed  information  on  the  coupling  between  the 
toroidal  and  poloidal  modes  of  hvdromagnetic  waves  can  be  extracted.  Figures 
4-6  show  the  global  structure  of  the  three  magnetic  field  components.  Their 
local  values  at  L  =  6-6  and  — 11°  <  A  ^  11°  lead  to  the  results  shown  in  figures 
2  and  3.  The  /t  axis  (from  0  to  025)  is  directed  along  the  magnetic  field  fine  and 
covers  the  region  from  the  equator  to  the  north  pole.  The  other  half  of  the  field 
distributions  can-  be  obtained  by  its  symmetrical  counterpart.  The  v  axis 
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Figure  5.  Normalized  magnetic  field  of  east -west  component  (»t  —  20.  k  =  4.  q  =  5). 


Figure  6.  Normalized  magnetic  field  of  north-south  component  (m  =  20.  k  =  4.  q  —  5). 

represents  the  outward  direction  of  the  field  line:  v  =  0-1  (r  =  10/?,)  is  at  the 
magnetopause,  while  v  =  0-5  (r  =  2RC)  is  at  the  ionosphere. 

The  global  structure  of  the  field  components  presented  in  figures  4-6  shows 
that  the  dominant  regions  of  the  wave  field  in  terms  of  its  amplitude  and  spatial 
variation  are  somewhere  near  the  inner  boundaries  of  t  he  region  inside  L  —  6  6. 
This  suggests  that  the  location  of  the  source  of  the  instability  responsible  for  the 
observed  Peo  wave  at  L  =  6  6  may  not  be  at  L  =  6-6.  Therefore  the  use  of  the 
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SCATHA  ion  data  to  deduce  if  the  drift  mirror  instability  is  the  cause  of  the 
observed  Pco  wave  may  not  be  conclusive.  Although  Takahashi  el  ill.  (1987) 
have  shown  that  the  instability'  condition  cannot  be  satisfied  at  L  =  6-6.  it  is. 
in  fact,  possible  that  the  condition  for  drift  mirror  instability  may  be  satisfied 
at  the  actual  location  of  the  source  of  instability,  which  is  near  the  inner 
boundary  of  the  magnetosphere  as  indicated  in  our  theoretical  analysis. 


5.  Summary  and  conclusion 

We  have  investigated  hydromagnetie  oscillations  of  the  geomagnetic  field  in 
the  earth's  magnetosphere  to  explain  low-frcquency  geomagnetic  pulsations.  A 
dipole-field  model  of  the  magnetosphere  is  adopted  for  the  analyses  of  the 
problems,  whereas  the  plasma  density  distribution  is  matched  to  in  silu 
measurements  with  satellites.  Considering  an  ideal  MHD  plasma,  a  hydro- 
magnetic  wave  equation  for  the  vector  wave  field  is  derived.  In  terms  of  the 
dipole  co-ordinates,  this  equation  is  then  decomposed  into  two  mixed-type 
coupled  partial  differential  equations  with  variable  coefficients  for  two  coupled 
hydromagnetie  wave  modes.  A  numerical  code  has  been  developed  for  solving 
these  two  equations  using  a  finite-difference  technique.  The  stability  of  the 
numerical  scheme  is  examined  by  both  the  Fourier  and  matrix  methods,  using 
the  following  procedure:  the  stability  conditions  are  first  determined  by  the 
Fourier  method  and  then  checked  by  the  matrix  method.  Since  two  different 
methods  are  used,  one  can  be  sure  that  the  stability  conditions  are  determined 
consistently.  This  code  is  employed  +o  .->tudy  an  implicit  boundary-value 
problem,  namely  finding  the  self-consistent  boundary  conditions  at  the 
magnetopause  so  that  the  wave  field  distribution  inside  the  magnetosphere 
matches  the  characteristic  features  drawn  from  the  data  information  obtained 
by  the  satellites.  Then,  on  the  basis  of  local  data,  we  can  reconstruct  the  global 
structure  of  the  stormtime  Pco  waves  in  the  magnetosphere. 

The  characteristic  features  of  the  hydromagnetie  waves  derived  from  our 
theoretical  model  are  found  to  be  consistent  with  those  of  the  observed 
stormtime  Pco  waves.  For  instance,  the  compressional  component  of  the  waves 
has  an  antisymmetric  standing-wave  structure  evidenced  by  a  node  near  the 
equator,  while  the  transverse  component  exhibits  a  symmetric  standing-wave 
structure  with  an  antinode  near  the  equator.  That  the  two  wave  components 
have  different  parallel  wavelengths  (referred  as  the  multiple-wavelength 
phenomenon)  is  also  seen  in  our  numerical  results.  The  calculated  relative 
phases  for  the  three  wave  components  agree  quite  well  with  observations.  The 
appearance  of  a  harmonic  wave,  which,  however,  cannot  be  reproduced  by  the 
present  model,  is  believed  to  originate  from  nonlinear  effects  of  the 
magnetoplasma. 

The  reconstructed  global  structure  of  the  micropuisations  can  provide  useful 
information  on  the  mechanism(s)  generating  the  stormtime  Pco  waves.  Our 
numerical  results  show  that  strong  variation  and  large  wave  structure  occur 
near  the  boundary  of  the  magnetoplasmas  spanning  the  magnetosphere  and 
ionosphere.  This  suggests  that  the  source(s)  of  instability  mechanisms 
responsible  for  the  observed  Pco  waves  are  located  in  that  region,  namely  near 
either  the  magnetopause  or  the  ionosphere. 

In  conclusion.  the  global  micropuisations  in  the  magnetosphere  can  be 
reconstructed  by  the  present  model.  The  success  of  this  theoretical  model  has 
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been  demonstrated  by  comparing  the  numerical  results  with  a  107!)  event 
(Takahashi  et  al.  1987)  for  which  the  stormtime  Pco  waves  were  measured  by 
the  four  geosynchronous  satellites  SCHATHA.  GOES2.  GOES3  and  GEOS2. 
This  model  has  also  identified  the  source  regions  of  micropulsations  to  be  at 
either  the  magnetopause  or  the  ionosphere.  Moreover,  it  can  correlate 
geomagnetic  activities  inside  the  magnetosphere  with  solar  wind  activity  in  the 
neighbourhood  of  the  magnetopause.  The  present  model  has  neglected  possible 
wave  dissipation  processes.  We  have  assumed  that  the  dissipation  of  waves  in 
the  magnetosphere  has  been  compensated  by  sources  located  at  the  boundaries. 
This  is  realized  by  the  fact  that  the  observed  waves  were  last  for  a  long  duration 
of  about  50  h. 

This  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research  Grant 
No.  AFOSR-88-0127.  The  numerical  work  was  performed  at  the  Pittsburgh 
Supercomputing  Center,  supported  by  the  National  Science  Foundation.  The 
authors  wish  to  thank  Professor  Jerry  Shmoys  for  critical  comments  on  the 
original  manuscript. 


Appendix.  Numerical  methods 

The  governing  equations  (2)  and  (3)  co!  le  a  set  of  mixed-type  coupled 
partial  differential  equations  (PDEs)  with  \  ..able  coefficients.  Equation  (2)  is 
a  parabolic  type  of  PDE.  while  (3)  is  of  elliptic  type.  The  two  equations  are  first 
converted  into  difference  equations  for  numerical  analysis.  This  is  done  bv  a 
finite-difference  technique  using  a  staggered  grid  with  grid  increments  A//  and 
Ac.  Both  of  the  employed  computational  grids  are  carefully  chosen  to  ensure 
second-order  accuracy  in  the  resulting  system  of  algebraic  equations.  A  fully 
implicit  scheme  is  used  for  this  numerical  algorithm.  The  resulting  finite- 
difference  equations  are 


s,  bi  ~ 2  Ar  aif  b4  -f  jp  a{w  b4  +  4  .  b^j  TP 1 
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Figure  7.  Dependence  of  the  minimum  Ar  on  in.  with  S/i  =  0-004. 
where  1  <  /  ^  /—  1  and  1  <_/<•/—  1. 
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Methods  are  adopted  to  establish  the  stability  conditions  of  the  finite 
difference  scheme  numerically.  The  Fourier  method  (Potter  1972)  is  first  used 
to  examine  parametrically  the  stability  conditions  of  the  designed  numerical 
scheme.  These  conditions  are  then  incorporated  by  applying  the  matrix  method 
(Ames  1977).  This  method  is  used  to  exemplify  the  validity  of  these  stability 
conditions.  Since  the  designed  algorithm  is  effective  for  a  parametric  study  of 
the  stability  conditions,  the  optimum  parameters,  which  can  increase  the 
accuracy  of  the  numerical  solution,  can  be  determined  together  with  the 
stability  analysis.  For  instance,  considering  the  case  of  longitudinal  mode 
number  m  =  100,  we  find  that  the  stability  conditions  of  both  methods  are 
satisfied  for  A/t  =  0-004  and  Ai>  ^  0-01334.  Hence  the  minimum  Av  =  001334  is 
used  in  the  numerical  analysis.  The  numerical  programs  for  the  coupled 
equations  (2)  and  (3)  and  the  stability  analysis  have  been  executed  on  the  Cray 
X-MP/48  of  the  Pittsburgh  Supercomputing  Center. 
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Using  the  same  procedure  to  perform  the  stability  analysis  for  different  in 
values,  the  dependence  of  the  minimum  A v  on  in  can  be  determined.  The  result 
is  shown  in  figure  7.  which  shows  that  the  minimum  Ar  decreases  rapidly  with 
in  for  20  <  m  ^  40  and  approach  a  constant  value  for  m  >  40.  Since  the 
minimum  Ar  is  smaller  for  larger  m.  the  numerical  solution  is  expected  to 
achieve  better  accuracy  for  the  case  of  layer  in. 
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